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Introduction 


The  complete  calculation  and  representation  of  sound  fields 
involves  considerable  difficulties  even  in  the  cases  of  simple 
radiator  arrangements.  ■'"For "this  reason,  no  systematic  treatment 
of  the  subject  is  available  to  the  author's  knowledge . > In  the 
following  treatment,  an  attempt  is  made  to  develop  the  basic 
formulae  and  concepts  and  to  apply  them  practically  in  a series  of 
simple  examples.  Due  to  the  extraordinary  simplification,  it  ap- 
pears justifiable  at  first  to  confine  the  calculation  to  a region 
at  a great  distance  from  the  radiator.  With  this  limitation,  the 
transition  from  the  non-directional  to  the  directional  radiator  is 
made  very  simply  when  the  formulae  for  the  non-directional  sound 
emitter  are  modified  by  an  additional  factor  designated  as  the 
directional  characteristic  or  the  radiation  factor.  Also,  with 
thi3  simplification,  the  calculation  of  the  sound  field  of  membranes 
(which,  in  contrast  to  the  piston  membrane,  do  not  need  to  have  con- 
stant vibration  amplitudes  and  may  possess  nodal  lines)  then  en- 
counters no  particular  difficulties. 

The  second  part  treats  the  representation  of  the  sound  field  in 
the  immediate  neighborhood  of  the  radiator.  The  case  first  investi- 
gated, that  of  two  point  radiators,  shows  that  even  here  substantially 
more  complicated  circumstances  exist.  One  is  therefore  forced  to 
describe  the  sound  field  at  every  point.  This  is  done  most  simply 
when  the  constant  pressure  amplitude  curves  are  drawn.  In  the  case 
of  the  circular  piston  membrane,  which  is  of  considerable  practical 
significance,  the  values  on  the  normal  axis  and  also  in  the  immediate 
neighborhood  of  the  membrane  can  be  very  plainly  stated.  From  the 
fact  that,  with  diameters  increasing  in  relation  to  the  wavelength, 
an  ever  increasing  number  of  null  values  and  maximum  values  result 
on  the  axis,  it  follows  that  the  piston  membrane  with  increasing 
radius  produces  no  sound  field  even  in  the  axis  neighborhood  compa- 
rable with  a plane  sound  wave.  In  the  cases  where  the  diameter  is 
not  too  great  in  relation  to  the  wave  length,  the  calculations  of 
the  sound  field  can  generally  be  carried  out  at  nearby  points.  The 
graphical  representation  of  the  curves  of  constant  pressure  amplitude 
provides  a simple  conspectus  of  such  sound  fields. 

The  spherical  radiators  which  are  treated  in  Part  Three  have  a 
particular  significance.  For  sound  velocity  amplitudes  given  on  the 
sphere,  the  sound  field  can  be  calculated  even  for  nearby  field  points 
It  is  shown  in  a series  of  examples  that  the  fonnulae,  which  at  first 
appear  involved,  are  quite  useful  for  practical  calculation. 
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In  order  to  avoid  a complicated  appearance  of  the  text,  the 
theoretical  derivations  are  given  only  in  so  far  as  is  necessary  for 
comprehension,  and  references  are  given  to  the  existing  literature 
for  the  remainder.  Also,  for  simplicity,  the  calculation  and  repre- 
sentation of  the  sound  velocity  amplitudes  are  avoided  and  only  the 
pressure  amplitude  is  used  to  describe  the  sound  field.  Thereby  the 
use  of  the  velocity  potential  is  completely  avoided  since  the  concept 
possesses  mainly  a mathematical  significance  and  according  to  experi- 
ence provides  no  pleasure  for  the  practical  physicist. 

All  considerations  are  based  on  a given  velocity  amplitude  of  the 
vibrating  membrane;  or,  what  amounts  to  the  same  thing,  a definite 
sound  velocity  amplitude  immediately  before  the  membrane  (imagined 
stationary)  is  prescribed. 


viii 


Part  One 


THE  SOUND  FIELD  AT  A GREAT  DISTANCE  FROM  THE  RADIATOR 
1 . Nondirectional  radiation 


We  start  with  radiators  of  zero  order,  i.e.,  with  radiators 
whose  radiating  surface  moves  cophasally  outwards  and  inwards. 

The  simplest  case  is  a breathing  or  pulsating  sphere  (Fig.  1-a). 
Here  the  velocity  amplitude  is  the  same  for  all  points  of  the  radi- 
ating surface.  In  general,  the  motion  of  the  radiating  surface  is 
described  by  the  behavior  of  the  velocity  amplitude  tt>  — wf'"' 
of  the  membrane,  w will  in  general  be  different  at  different  points 
of  the  radiating  surface. 

It  is  advantageous  to  introduce  the  velocity  amplitude  w in- 
stead of  the  displacement  amplitude  a (for  which  the  relation 
w~'Znn a exists)  3ince  the  latter  achieves  a particular  signifi- 
cance only  in  very  infrequent  cases. 

We  first  consider  the  series  of  radiator  forms  represented  in 
Fig.  1.  I he  radiators  1-a  to  1-d  are  rotationally  symmetric  about 
a vertical  axis  through  the  center.  The  radiators  1-e  to  1-g  are 
likewise  rotationally  symmetric  or  have  a rectangular  cross  section 
ljJce  1-h.  The  rest  position  of  the  membrane  i3  indicated  by  a 
heavy  broken  line,  the  position  of  maximum  displacement  is  dotted, 
and  the  rigid  supporting  wall  is  indicated  by  a solid  line.  The 
maximum  amplitudes  of  the  periodically  (sinusoidally)  moving 
membrane  are  indicated  by  the  greatly  enlarged  arrows.  The 
linearity  of  the  general  wave  equation  depends  on  the  assumption 
that  the  pressure  and  velocity  amplitudes  of  the  sound  field  be 
small.  This  assumption  is  fulfilled  in  practice  except  for 
unusually  strong  pressure  variations  (o.g.,  explosion  waves). 
Furthermore,  the  separation  of  those  membranes  vibrating  sym- 
metrically with  respect  to  the  middle  plane  (e.g.,  in  ^'igs.  l~e 
and  1-f)  is  assumed  so  small  that  the  membrane  zero  position  can 
be  regarded  as  practically  coincident  with  this  middle  plane. 

It  then  turns  out  that  the  sound  field  can  be  very  simply 
determined  if  the  two  following  assumptions  are  fulfilled! 

A.  TTie  dimensions  of  the  radiator  in  every  direction  are 
small  compared  to  the  wave  length. 
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B.  The  field  point  lies  at  a sufficiently  great  distance 
from  the  radiator. 

From  the  general  formulae  (to  be  explained  later)  it  then  follows 
that  the  behavior  of  the  sound  pressure  is  represented  in  the  simple 
form: 


p0  — p0ei('a‘-ir+*i2) 


where  the  sound  pressure  amplitude  p0  is  given  by 

c ■ a • F ■ vi- 


Po  = 


2 Al- 


and the  total  radiated  sound  power  10  by 

L0  = (pll%ca) =4  c • an  ' — -is— 


Aa 


(1) 

(2) 

(3) 


We  see  that  the  pressure  amplitude  is  described  in  a very 
simple  manner  by  the  standard  quantities  and  that  each  of  these  quan- 
tities appears  linearly*  ?he  significance  of  the  various  symbols  in 

as  follows: 


/t 
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Jig.  1 Simple  radiator* 


<3  is  the  density  of  the  medium 
c is  the  velocity  of  sound  in  the 
medium 

F is  the  total  radiating  surfacd^ 
w is  the  mean  velocity  amplitude 
m of  the  radiating  surface 
r is  the  distance  of  the  field 
point  from  the  center  of  the 
radiator 

X is  the  wave  length 
k is  the  wave  number  ( 2n/\  ) and 
co  is  the  circular  frequency  2*n . 

The  product  o • a is  termed  the  specific 
acoustic  resistance.  Fur  plane  waves  its 
magnitude  is  determined  by  the  ratio  of  the 
sound  pressure  amplitude  to  the  sound  velocity 
amplitude.  For  the  two  values  with  which  we 
are  principally  concerned,  those  for  air  and 
water,  we  have  43  and  1,5  x 10*  absolute  C.G.S. 
units  respectively  as  the  value  of  c *0  . 

In  the  general  case  the  velocity  amplitude 
w(x,y,z)  varies  for  the  individual  points  of 


^ For  membranes  which  vibrate  in  a rigid  wall  we  will  denote  the 
8urfaoe  radiating  into  the  halt  space  by  F so  that  the  "2'*  in  the  de- 
nominator of  equation  2 drops  out. 


— 


the  membrane. 


In  this  case  wm  is  given  by 

wrn  — y jw(x,y,z)dF.  ^ 

F 

The  influence  of  the  membrane  is  given  by  tne  product  F ■ u-,„ 

If  we  define  the  volume  swept  by  the  oscillating  surface  during  each 
half  period  as  the  deformation  volume  of  the  radiator,  then  two  radi- 
ators of  equal  frequency  have  equal  deformation  volumes  if  the  product 
F ■ wm  has  the  same  value  for  both.  This  is  so  since  the  displacement 
amplitude  a(x,y,z)  differs  from  the  velocity  amplitude  w(x,y,z)  only 
by  a constant  factor  (w  = 2tm*a)  . 

From  the  fact  that  the  quite  varied  vibration  forms  of  Figs . 1-a 
to  1-h  produce  the  same  sound  field  for  equal  values  of  F ■ w,„  , we 

will  quite  generally  conclude:  For  radiators  which  are  small  compared 
to  the  wave  length  and  whose  membrane  motion  takes  place  in  phase, 
the  sound  field  at  a comparatively  great  distance  i3  determined  by 
the  formulas  (1),  (2),  and  (3)*  This  means  that  under  assumptions.  A 
and  B two  sound  radiators  for  which  the  frequency  and  the  product  F ■ 
are  identical,  produce  the  same  sound  field,  independently  of  the 
separate  values  of  displacement  amplitude  and  radiating  surface.  In 
Fig.  2 three  different  radiators  corresponding  to  Fig.  1-e  and  1-f 
are  represented  with  equal  deformation  volumes,  i.e.,  equal  values 
of  F ■ w,„  , 

The  shaded  surfaces  represent  the  maximum  amplitudes  in  one 
direction  of  the  vibrating  membranes  in  two  mutually  perpendicular 
planes.  The  velocity  amplitude  corresponding  to  that  of  Fig.  2-b  is 
given  by  w = (i  -oa/(?r.)a.  Then  for  the  three  (3)  cases  of  Fig.  2 the 
equal  values 

j qU,  2 nj  (i-e®/eii)Ve 

0 

result  for  Fwm  . 

It  is  advantageous  to  set  up  a simple  example  as  a standard  form 
since  from  it  the  sound  pressure  can  ba  specified  under  quite  general 
conditions  because  of  the  linear  influence  of  all  of  the  independent 
variables . As  such  a standard  fonn  we  consider  a piston  membrane  in 
a rigid  wall  which  produces  a sound  pressure  amplitude  of  1 dyne/cm2 
at  a frequency  of  800  cps  on  the  normal  axis  at  100  cm  distance  (which 
corresponds  to  the  sound  pressure  amplitude  for  normal  speech  immedi- 
ately in  front  of  the  mouoh  of  the  speaker) . <Ve  thus  find  that  with 
a velocity  amplitude  of  10  cm/sec  the  surface  radiating  into  the  half 
space  must  be  10  cm2.  It  is  to  be  particularly  noted  that  formula  2 
is  valid  for  all  plane  beamed  radiator  arrangements  if  the  field  point 
lies  on  the  normal  axis  at  a sufficient  distance  from  the  radiator. 


I 2 

;j  (*U  n l 


a‘n 


(4a) 
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This  will  be  explained  in  the  next  chapter.  Therefore  a radiating 
surface  of  100  cm^  also  produces  a sound  pressure  amplitude  of  1 
dyne/cm?  at  a field  point  10  m distant  on  the  normal  axis.  Finally, 
the  formula  is  likewise  valicf^for  a spherical  radiator  of  zero  order 
for  a sufficient  distance  of  the  field  point. 

The  corresponding  sound  power  per  cm^  for  the  standard  form  is 
(according  to  3) 


4 nr*  2-43 


= 1,16  -10-4  erg/sec  cm8  = 1,16  • 10~»  Watt/cm2 . 


(4b) 


In  water  the  sound  pressure  amplitude  v n * 800  dynea/cm^ 
results  for  the  same  assumptions  (with  the  same  frequency  - therefore 
a different  wave  length)  and  the  sound  power  becomes 


inr 2 


82  • 10* 

2 -Ts . io6  ~ 2,13.erg/8eccm2  - 2,13- 10~7  Watt/cm2. 


(4c) 


One  thus  sees  what  a decisive  influence  the  "sound  resistance" 
(specific  acoustic  resistance)  exerts.  In  order  for  the  membrane  to 
impress  the  same  velocity  amplitude  in  water  as  in  air  at  the  same 
frequency  184-fold  power  must  be  produced.  On  the  other  hand  for  a 
membrane  which  is  to  radiate  at  an  equal  wave  length  the  same  power 
in  air  as  in  water  it  follows  that 


u>„ 


i/ubJW 

\ 43  ~~ 


cv^CO, 


(4d) 


i.e.,  the  velocity  amplitude  w,  of  the  membrane  in  air  is  then  60 
times  as  great  as  the  velocity  amplitude  w.:.  in  water  while  the 
displacement  amplitude  of  the  membrane  in  air  is  about  13  times  as 
grea^  as  that  in  water.* 


.2, 

v'  This  is  so  with  the  provision  of  Note  1. 

* That  the  number  13(60/4.5)  should  be  changed  to  270(60*4.5) 
is  shown  by  the  following; 

Using  the  statement  mRda  at  the  end  of  the  paragraph  following 
eq,  4,  one  has  from  eq,  4dt 

^ 60  -»■  ~~  -v  to 

~ >'ur  _ c 

Sinoe  A is  to  ue  constant  for  the  two  cases,  we  have  ^ / 

, r\  hi  4 

or  *lur~  ~ **-•  — V.  if  . Banco  — - ^ be? 1 V.  ^ — 2 70, 

*4  C.  u •** 
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2.  Beamed  Radiation 
A.  The  directional  characteristic 
(a)  For  a definite  frf  : .ency 


We  now  drop  assumption  (A)  but  i os©  instead  the  condition  that 
the  membrane  vibrates  as  a double  membrane  cophasally  toward  the  two 
sides  - the  middle  plane  at  every  moment  of  the  vibration  being  the 
symmetry  plane.  Examples  of  such  membranes  are  represented  in  Figs. 
l~e  and  1-f.  oince  the  whole  process  above  the  symmetry  plane  repre- 
sents the  process  below  the  symmetry  plane,  a rigid  wall  can  be  in- 
troduced at  the  position  of  the  symmetry 
plane  without  changing  the  sound  process 
in  any  way.  In  place  of  the  double-side dly 
acting  membrane,  one  can  then  just  as  well 
imagine  the  one  half  of  the  sound  field 
produced  by  a simple  membrane  vibrating 
in  a rigid  wall.  This  is  significant 
since  a membrane  vibrating  in  a large 
rigid  wall  appears  in  practice  rather  fre- 
quently. without  particular  emphasis  in 
every  case  it  is  assumed  in  the  following 
examples  that  the  radiating  surface  con- 
sists of  such  double  membranes  vibrating 
symmetrically  to  the  central  plane  or  of 
membranes  vibrating  in  a rigid  wall  for 
which  only  the  one  half  space  is  considered  which  is  cut  off  by  the 
rigid  wall  and  into  which  the  one  half  of  the  membrane  with  the  sur- 
face F radiates.  If  we  imagine  that  .the  elementary  waves  produced 
by  the  individual  elements  of  the  radiator  act  together  (as  point 
sources)  at  the  field  point,  it  is  then  clear  tliat  the  concurrent 
action  in  uniform  phase  at  far  distant  field  points  prevails  for 
all  directions  as  long  as  the  individual  elements  are  at  a distance 
d from  one  another  which  is  small  compared  to  X Otherwise, 
the  elementary  waves  arrive  at  the  field  point  with  a phase  differ- 
ence which  is  no  longer  negligible  (Fig.  3)  with  a resulting  lack 
of  spherical  symmetry.  ..o  can  take  account  of  the  latter  conditions 


Fig. 


3,  Nondirectional 
radiation. 


^ The  qualifying  condition  boro  would  rather  seem  to  be  that 
the  over-all  dimension  of  the  radiating  ensemble  be  small  compared 
to  ~ The  statement  in  the  text  could  be  construed  to  moan  this 
If  ^ wero  ta'<en  to  represent  the  distance  between  any  two  elements 
(including  the  two  most  distant  elements)  instead  of  being  merely 
the  distance  between  two  adjaeont  elements  as  is  shown  in  Fig.  3. 
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by  applying  the  directional  characteristic  SR  . hence  in  place  of 
(2)  we  have  the  formula 


V = Po  • SR » 


(5) 


wherein  p0  has  the  value  given  by  (2)  while  SR  is  defined  as  the 
integral  over  the  surface  F: 

SR  = f w(x ,y)  • dF . ( 


Here  the  rigid  wall  with  the  membrane  (in  the  position  of  rest) 
is  thought  to  be  in  the  plane  of  the  coordinate  system.  Since  the 

XY-plane  is  the  symmetry  plane  for  the 
sound  field,  we  can  confine  our  con- 
siderations to  the  space  lying  above 
the  XY-plane  and  therefore  disregard 
the  symmetrically  downward  oscillat- 
ing part  of  the  double  membrane  (Fig. 
4). 


Let  a , {3  and  y be  the  direc- 

tion angles  for  the  corresponding 
field  point  line,  x and  y the  co- 
ordinates of  the  point  belonging 
to  the  element  of  Integration  aF , 
w(x,  y)  the  corresponding  velocity 
amplitude  and  let  k = 2 w/X  . we 
note  that  if  the  wave  length  X 
is  great  compared  to  x and  y - i.e., 
compared  to  the  dimensions  of  the  membrane  - the  exponent 
the  Integral  approaches  zero  so  that 


Fig.  4.  Beamed  radiation. 


.great 

under 


•N  " »)*r 


wm . F 
V 


(6a) 


and  therefore  by  (4)  becomes  ecpial  to  1 as  was  to  be  expected,  ^s 
has  already  been  indicated,  it  is  essential  to  note  ttiat  sJt  like- 
wise bocomos  1 for  a large  F if  the  field  point  line  coincides 
with  the  4-axis  (since  a ^ 9ir  and  /I  ---  W° ) . 

The  validity  of  (5)  results  from  the  general  basic  relation 
(for  a surface  F on  an  infinite  plane)  given  by  Rayleigh v^s 


The  Theory  of  Sound,  Sec.  273. 
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(7) 


<p  = 


1 C6<p<i 
2 nj  dn 


Sty 


-dF. 


When  we  introduce  the  pressure 


P = 


(*) 


in  place  of  the  velocity  potential  <p  and  replace  the  normal  velocity 
$2  by 

on 


to  f=  w(x, 

we  immediately  obtain 

p = w{x. y) e'  f dF. 

>• 

V/e  will  now  consider  the  sound  field  at  a great  distance  from 
the  radiating  surface  and  accordingly  assume  that  the  field  point  P 
is  at  a sufficient  distance  R from  0 30  that  the  connecting  lines 
from  it  to  the  individual  radiator  elements  dF  can  be  regarded  as 
parallel.  (A  more  accurate  formulation  will  be  given  later  (p.  59)») 
If  we  then  drop  a perpendicular  from  dF  which  is  determined  by 
the  coordinates  (x,  y)  to  OP  intersecting  OP  at  Q then 

OQ  — x cos  a + y cos  / i and  r = R — (.»•  cos  a + y cos  /}). 

Since  u > x cos  -x  y cos  (i  we  can  replace  r by  R in  the 
denominator  and  obtains 


(9) 

(10) 


\)  ~ e»<wt-  iH)  I' w(Xi  y)eU(tco»«  + vcmfi) (if 
y 


(11) 


and  from  this,  considering  (2),  (4),  and  (6),  there  results  the 
formula  (5). 


v Ho re  the  symbol  r designates  the  distance  from  the  field 
point  to  the  surface  element  whereas  the  same  symbol,  as  employed 
in  equations  (l)  and  U),  signifies  the  distance  from  tho  field 
point  to  tho  oenter  of  tae  radiator  (seo  tho  list  of  definitione 
on  P*  2). 
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We  shall  assume  for  what  follows  that  the  radiating  surface  or 
the  radiating  system  possesses  a point  of  syneaetiy,  which  coincides 
with  the  coordinate  origin,  such  that  w(x,  y)  - w(-x,  -y) . 

(That  is,  tw  membrane  points  situated  in  image  fashion  with  respect 
to  the  membrane  point  of  symmetry  have  the  same  amplitude  and  phase.) 
The  integral  (ll)  may  then  be  written  in  the  real  form;  * 


SR- 


F 


w(x,  y)c os  [k(xcostx  + ycosft)]dF . 


(Ha) 


We  will  furthermore  confine  the  field  point  line  to  a definite 
plane  (the  measuring  plane  (die  Peilebene))  which  is  to  coincide 
with  the  ZY-plane . Then  since  a = 90°  and  p 4-  y = 90°  (Fig.  4) 


SR 


w{ x.  y)cos[kys\ny]dF . 


(12) 


In  many  cases,  interest  in  the  behavior  of  sr  does  not  extend 
to  the  whole  measuring  plane  but,  particularly  with  sharply  beamed 
radiators,  only  to  the  immediate  neighborhood  of  the  principal 
maximum  at  y = « 


For  this  range  cos[£ysinyj  can  be  replaced  by  l - 4 Wy*  sin2y  and 
it  is  found  that 


SR  = 1 — 


ia8insy  /' 

T wj?J 


w(x}  y)yHF . 


F 


(12a) 


If  we  imagine  the  radiating  surface  element  dF  at  the  point 
(x,  y)  of  the  membrane  replaced  by  an  element  of  mass 

dvi^~U'^dF 

tit  Jr 


T,  = "~;rl  w(x,y)yhlF  - J t/dm 

F F 

denotes  the  moment  of  inertia  vdth  respect  to  the  X-axis  of  the 
stationary  membrane  thus  covered  vdth  mass  (total  mass  1),  and 
for  small  values  of  y we  can  write 


(12b) 

(12c) 


Apparently  this  is  the  real  form  of  eq. 


6 instead  of  eq.  11. 


-a- 


fe-'OSv 


9?  = 1 — £ A2sm2y  • Tx . \ \ JJ 

This  formula  is  often  suitable  for  a quick  determination  of  the 
directivity  (die  Peilschirfe)  of  radiator  arrangements. 

As  an  example  we  will  investigate  four  simple  radiator  arrange- 
ments (Fig.  5): 

1.  Two  point  sources 

2.  The  circular  line  densely  filled 
with  radiators 

3.  The  rectangular  piston  membrane 

4.  The  circular  piston  membrane 

Z 

tie  find  the  moment  of  inertia  T and  the 
corresponding  approximation  formulae  without 
7 -s  difficulty! 


1.  Tx  = r2,  SR,  •—  1 — -g  k2sin2y  • r2, 

2.  7\=  y,  91,= 


(13a) 


3.  T3=  3 


4.  T4  = 


3RS  = 1 — fcJ8inay  • T 


SRt  = 1 - -g-  &2ainay  • i- 


For  these  four  cases  the  directional 
characteristics  may  be  generally  calculated 
by  carrying  out  the  integration  in  Eq,  (6)^ 
There  results  (see  Fig.  5^ 


e<*r»tai-+e-<*r»tar  2nr 

SRt  — — cos  sin  y 


Fig.  5.  Calculation  of  3ig  = .,2  fe<*,JrC0l,!‘-,  »™*0>cfs 
the  directional  Ljnr 


characteristic 


-=  i r^vdf  = sin  y) , 


Cleotr.  Nachr.  Techn. , Vol.  4 (1927),  pp.  239-253. 


IS- 


v From  equations  (l3a),the  following  meanings  of  r are  to  be 
inferred!  r = one-half  distance  between  elements,  for  two  point 
sources , 

X = radius,  for  ring  and  piston t 
r - one  half  the  side  of  the  square  piston  membrano. 
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•nyt .<  vetyfys  v-'r  <yr'. 


+ 7 

S' 


% 


irl**'*»dy=  2nr 


fixer  . \ 

h-8my) 


(16) 


-ainy 


. f f ,f  ■'.(-»-“>•) 

SR4  = -fn  / gdo  I ei^nnyoo^dep  = -J  gdpJ0(kgainy)  =2--^-— 


0 


- am  y 


(17) 


'tilth  the  abbreviation  -y~  sin  y — sin  y — x , it  is  found  that 


2 </,(x) 


iRi  = coax,  8i2  = /0W,  <K3  = ™,  = ^ 


(17a) 


where  j m or  j^x)  signify  Bessel's  functions  of  the  zero  or 
first  oraer  respectively.  By  the  vjell-knovm  series  development  the 
correctness  of  the  approximation  formulae  may  be  immediately  con- 
firmed. One  thus  obtains 


cos  x ~ 1 — 2 + 


J0M  ■--- 1 - v + 


sin  x _ x4 
~x~~  ~ 1 ~ 6 + 


i (17b) 


2 J,(x) 


- 1 


In  order  to  obtain  a general  measure  for  the  directivity  we 
will  inquire  as  to  the  angle  <p  in  the  measuring  plane  (ZY)  which 
the  field  point  line  forms  vdth  the  Z-axis  when  the  sound  energy 
is  reduced  to  half  the  maximum  value,  (i.e.,  the  amplitude  is  re- 
duced from  1 to  1 / y 2)  (Fig.  6). 

We  win  denote  this  angle  as  the  half-value  beam  width  in  this 
measuring  planed.  A simple  calculation  then  yields  vdth  sufficient 
accuracy  for  practical  purposes  the  fonowing  simple  relations: 


^ It  is  here  assumed  that  there  is  some  directivity  so  that  the 
half-value  width  amounts  at  most  to  d0°. 


*><k  V'T^t  *••  • 


<P1  = 

15°  A/d 

for 

fz  = 

20°  i/d 

for 

tyz  = 

25°  i/d 

for 

<p*  = 

30°  i/d 

for 

The  error  here 

is  ; 

portionally  for  smaller  values  of  X 


(18) 


d and  decreases  pro- 


The  angles  15°,  20 °,  25°,  and  30°  are  then  characteristic  of 
the  directivity,  "e  will  emphasize  this  by  a special  term  - the 
directivity  coefficient  (das  Peilmass).  The  directivity  coefficient 
is  in  general  different  for  every  measuring  plane.  In  order  to 
obtain  the  directional  characteristic  for  a measuring  plane  com- 
pletely, the  characteristic  functions  will  be  plotted  as  a function 

of  x and  a variable  scale  will  be  placed 
below  the  x-axis,  which  permits  reading 
off  the  directional  characteristic  for 
each  value  of  d/X  and  each  angle  v . 

Thus  the  complete  behavior  of  the  direc- 
tional characteristic  as  a function  of  y 
in  the  above  four  examples  and  for  all 
values  of  d/X  which  lie  between  0 and  g 
can  be  immediately  ascertained  from  Pig.  7* 
If,  for  example,  the  directional  character- 
istic for  the  densely  covered  circular  line 

8?-y0(ysiny)  is  desired  for  d/X  = 3, 

and  y — 40°  , one  follows  the  line  d/X  = 3 
to  its  intersection  with  the  40°  line  and 
reads  off  the  corresponding  ordinate  on 
curve  2.  It  is  thus  found  that  9?  = +0,15 
(see  the  broken  line  in  ^'ig.  ?)• 

As  an  additional  example,  we  mention  the  straight  line  group 
which  consists  of  a number  of  uniform  nondirectional  individual 
radiators  arranged  at  equal  intervals  on  a straight  line,  due 
to  its  special  significance  in  practice,  we  will  investigate  this 
arrangement  very  closely.  Let  the  radiators  lie  on  the  Y-axis 
at  a distance  d from  each  other  and  let  them  have  the  coordinates: 


Pig.  6.  Definition  of 
the  beam  sharpness. 


yp  > Up  - 1 > • • • y - - 0 * y - v » 


(18a) 


^ere  the  coordinate  origin  is  to  coincide  with  the  middle  radi- 
ator for  an.  odd  number  of  radiators  {»  - 2</  ■ 1)  and  is  to  lie  midway 
between  the  middlemost  radiators  for  an  even  number  cf  radiators 
(n  = 2 q)  • the  first  case  one  would  then  set  d q ■ i,  p i = q 2 
etc.,  while  in  the  second  case  one  would  set  p - q j,  „ . i „ .i  Qtc. 
(Pig.  8), 
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Tig.  7.  General  representation  of 
the  directional  characteristic. 

1.  Directi*  ;al  charaeterictic  of  the  system  consist- 
of  tvo  radiator*  (Distance  d).  2,  Directional  char- 

eaterisMc  of  the  densely  oovered  circular  line  (Di- 
ameter dj.  3.  Directional  characteristic  of  the 
densely  covered  rectangular  eurfaee  (Length  d)  in  the 
symmetry  plane.  4.  Directional  characteristic  of  the 
densely  covered  circular  surface  (Plameter  d). 


In  general,  then,  for  an  arbitrary  index  l , the  relation 
yl  = ld  is  valid.  If,  for  abbreviation,  vne  set 

„iJ:d silly  _ , 

(18b) 

vre  then  find  the  directional  characteristic  from  (6)  when  we  replace 
the  integral  by  a summation  sign 

m = + p 

^-~Sz'"  a&) 

m — — p 


The  geometric  series  on  the  right  may  be  easily  summed  and  one 
obtains: 


3J  = 


1 


1 2"+  i _ 2- <"+  *’ 

n 


2—1  n z\  _ z-  \ 


from  which  follows  (when  eiil,8in>'  is  again  introduced  for  z) 


5R  = 


sin  ^>  + 

yjfcrfain^ 

n sin 

\kd~~Y~" 

Tsmyj 

(18d) 


(18e) 


Since  for  odd  numbers,  <p  + -4;  is  replaced  by  9~i  - % and  for 

even  numbers  <j>  + -1  is  replaced  by  cf  =.  Jl  , the  general  result 
* * ^ z 


is  that 


SR 


nnd  . 

x s,1‘  ^ 

7t  d . 
•smy 


(19) 


- where  k=y  is  substituted  for  k. 


In  order  to  obtain  a clear  understanding  of  the  general  be- 
havior of  the  directional  characteristic  for  different  values  of 
d/X  we  first  consider  the  function  given  by 


, Kin  w y j (20) 

| U Hill  l/i 


.ml  , . 

in  which  :■){  has  been  replaced  by  r and  ( j8my  ) by  </  . The 
curve  given  by  (20)  in  the  polar  coordinates  (r,  9 ) is  now  easily 
understood.  The  principal  maximum  is  attained  at  f 0 and  v>  « n 
(and  multiples  of  ? ).  In  between  lie  n - 2 secondary  maxima 
which  are  separated  by  the  zero  positions  2, .. . «-i). 
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If  a straight  line  is  drawn  parallel  to  the  ordinate  axis  at 
the  distance  l/n  > it  touches  the  curve  at  the  points  which  cor- 
respond to  the  angles 


3 

2n  71 ' 


2n-_l 
2 71  ' 


(20a) 


It  is  seen  that  the  positions  of  the  secondary  maxima  are,  to  a good 
approximation,  defined  by  the  angles 


3 

^ ~ 2 n 71  ’ 


2 n - 3 
2 n : 


(20b) 


(due  to  the  synmetry  with  respect  to  the  ordinate  axis  we  can  limit 
ourselves  to  the  values  of  <p  & n ).  Fig.  9 shows  the  behavior 
for  n = 6 « here  the  broken  line  curve  is  an  ellipse  with  the 

semi-axes  1 and  1/6  on  which 
lie  all  the  maximum  values  of 


z 


sin  Gy 
G sin  y 


(20c) 


The  secondary  maxima  (in  the 
first  quadrant)  coincide,  to  a 
good  approximation,  with  the 
contact  points  between  the  el- 
lipse and  the  straight  line 
drawn  parallel  to  the  ordinate 
axis  at  the  distance  l/6. 


Fig.  8.  For  the  calculation  of 
the  straight  line  group. 


For  a more  ",cc urate  deter- 
mination, the  equation 


n tg  = tg  w 99  (21) 

resulting  from  must  be  solved,  ^his  may  be  done  simply 

if  equation  (21)  is  written  in  the  form 


95  = n arc  tg(ntg9>). 


(22) 


The  approximation  values  just  found  are  substituted  in  the 
right  hand  side  of  (22)  and  a better  approximation  is  obtained.  Then 
this  better  approximation  is  again  substituted  in  the  right  hand  side 

of  (22)  eto.  If  fo  ~ 2 ~2 » ' 71  i-s  substituted  as  a first  approximation 
in  the  right  hand  side  of  (22),  0 then  results  as  a second 


Pig.  9.  Marking  the  positions  end.  magnitudes  of 
the  secondary  maxima  of  the  directional  character- 
istic for  the  straight  line  group  of  six  radiators. 


approximation  where  e0  is  given  by 


tg  n e0  = 


»tg?>o 


(22a) 


If,  e.g.,  n = 6 and  m = 1 then  %-45°  and  from 


tg6«0==J  (22b) 

f„=*  i°  35'  so  that  <Pi  — 43°  25' , The  third  approximation  yields  e,  = 0°5' 
and  <p„  = ?>i  - «j  = 43°  20'  , Thus  6 tan  n-  5.66  agrees  with 
tan  6 <p»  = 5.67  to  within  an  error  of  less  than  1$. 

If  we  denote  the  coordinates  for  the  maximum  value  by  rm,  fm 
then  (20)  and  (21)  must  be  fulfilled  simultaneously.  Therefore 

JV^sm'Yn  ~ aina«9?m,  (23) 


\2*n<pm>  (24) 


and  dividing  (23)  by  (24): 

r!„  «>h>ih  cmin<Fm  (24a) 
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Adding  (23)  to  (24a), 


rm  [ji2  sin2  fm  + cos29?m]  = 1 . 


(24b) 


That  is,  the  maximum  values  all  lie  on  an  ellipse  with  the  semi- 
axes 1 and  1/n.  In  order  to  perceive  the  dependency  of  the  directional 
characteristic  on  d/X  and  y , we  have  only  to  consider  how  the 
discussed  r,  <p  curve  is  affected  by  the  transformation  SR  = r and 
cp  = sin  v into  the  r,  y plane.*  Due  to  symmetry,  we  can  confine 

ourselves  to  the  first  quadrant  |o  <;  y £ 0 ) and  from  the  corres- 
ponding range  can  state  the  positions  and  number  of 

the  principal  and  secondary  maxima.  The  size  is  manifestly  not 
changed  by  the  transformation. 


If  we  compare  two  arbitrary  straight  line  groups  I and  II  where 
I is  characterized  by  Mj.rf,,/,,,  and  II  by  »2,rf2.  h then 


ain 

?!•!  sin  ’ 
sin 

7i28iri9?2  ’ 


jiiit  . 

8lny» 

1 1 'L  . 

7 a ^ «ny. 


(24c) 


Since  one  is  concerned  with  only  small  values  of  y in  estimating 
the  directivity,  sin  <p  can  be  replaced  by  <?  . And  since,  for 
small  values  of  cp  , the  inequality 


sin  »,</!,  ^ sin  M,  Vi* 
«i9’i  "iVs 


(24d) 


is  fulfilled  when  and  only  when  »,  <px  < ,«2  , it  then  follows  that 

the  directivity  of  II  is  greater  than,  equal  to  or  less  than  that 
of  I according  as  »,pi2//2  is  greater  than,  equal  to  or  less  than 
«,r /,/z-i  . Viith  the  same  frequency  and  total  length  the  more  densely 
covered  group  therefore  has  a smaller  directivity  (broader  beam). 
And  two  groups  with  a different  length  and  a different  number  of 
radiators  can  possess  the  same  directivity. ' For  example,  since 

' ’X, ^ the  groups  w,  ^ Ik,  </,  =-  A/(>  and  n,  — 3,  rf2  - i have 

equal  directivity  while  their'  base  lengths  have  the  ratio  of  17 
to  12. 


For  a number  of  directional  radiators  (in  a rigid  wall),  the 
calculation  is  substantially  simplified  when  the  radiators  are 
uniform  (so  that  each  radiator  by  itself  would  produce  tne  same 
directional  characteristic) . If  we  imagine  a nondirectional  (point- 
source)  radiator  in  place  of  each  Individual  directional  radiator 


Apparently  the  , y plane  is  meant  here. 


-16 


„ ‘ 


and  call  the  directional  characteristic  of  the  system  thus  composed 
of  nondirectional  individual  radiators  SRX  , and  the  directional 
characteristic  of  each  individual  radiator  «tt2  , then  the  direct- 
ional characteristic  S)i  of  the  whole  arrangement  is  simply  given 
by  the  product,  i.e.,  9i  = &', -iH2. 

For  example,  the  directional  characteristic  of  an  arrangement 
in  a rigid  wall  consisting  of  two  identical  circular  piston  membranes 
of  radius  r and  with  a center- to-c enter  distance  a in  a rigid  wall 
turns  out  to  be 


SR  = cos 


na 

/. 


cos  /? 


2 J,  (krsmy) 
k rainy 


(25) 


»ie  prove  this  theorem,  first  expressed  by  Bridge"7,  as  follows: 

Let  there  be  given  in  the.XY-plane  a system  I of  n nondirectional 
(point-source)  radiators  with  the  coordinates 


2/1) - (^.2/i) 

• • • (*»,  2/;,); 

(25a) 

the  velocity  amplitudes 

, «>i  . 

. . «>;, 

(25b) 

and  the  radiating  surfaces 

F[ , Fi  . 

• • n • 

(25c) 

and  for  this  system,  set 

w\F\  -|-  w'3F’s+ 

■ + u''„F'„  ---  ,-f 

(25d) 

Furthermore  suppose  a system  II  of  m nondirectional  radiators 
correspondingly  characterized  by  the  quantities 


(•*’|  , U |)>  {•*■'■}  >//*)<  • • • (•>',.! Urn)  • 


(25©) 


, H’j  . . . II’,,, , 


(251) 


W'/  J.U/  t,:  t 

* I » * J»  * ‘ ^ HI  > 


«■?*?+  •••  -i-  <«',  .4' 


(25g) 

(25b) 


^ Sec  H.  Poincare,  Thoorie  de  la  lumiero,  p.  158. 
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Now  in  place  of  each  radiator  of  system  II,  the  system  I is  to 
be  set  so  that  the  new  system  III  is  correspondingly  characterized 
by  the  quantities 


(x'' + x\ , y; +1/,),  (.i-;; + 4 , ?/; + >/,)...  (x" + , y;  -1-  y'„ ) 
w"-w[,  w"  w'i , ...  w"  w'„ 

f;f\  f;  f:  ...  f;f:, 


p— 1,2,3  . . .m. 


(25i) 


Th6  directional  characteristic  3?"'  of  the  system  III  results 
then  from  the  general  formula  (6)  where  the  integral  must  be  replaced 
by  the  summation  sign.  Thus 


9T« 


W II 

1 X 1 V,,,"  F" F'  r' k I W + *<i>  '’m ''  ii'i’  + W)  ™» A] . 
A' A"  jFj  i’  >’  1 1 

11=1  4^1 


(25j) 


while  the  directional  characteristic  »{'  of  I is  given  by 

>){'  (25k) 

4 1 

and  the  directional  characteristic  'J{"  of  II  is  given  by 


r 


Hi 

1 N1 ,1," ‘tV c,,R  ■'  1 iv'-c 
4"  “r  ‘ 


'P'] 


(25m) 


p-  t 


From  this  it  follows  immediately  that  9i"‘  — iK'-iH"  . This 
theorem  can  serve  to  derive  another  group  arrangement  whose 
directional  characteristic  possesses  a simple  form  from  a simple 
straight  line  group  arrangement.  If  we  let  the  directional 
characteristic  of  a straight  line  group  which  consists  of  n 
nondirectional  radiators  with  equal  surface  placed  at  equal  in- 
tervals on  the  Y-axis  with  the  velocity  amplitudes  «’t.  «’*.  . . . «■„ 
be 


(25n) 


and  denote  in  particular  the  directional  characteristic  of  n unit 
radiators  by 


then 


|l,  1,  ...  i| 


xiu  M If 
n silk «/•  ’ 


(26) 


3. 


II.  I] 


Mill '•!(/' 
2 Kill  </' 


(26a) 
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If  we  now  replace  each  radiator  of  this  group  consisting  of  two 
elements  by  the  same  group,  then  from  the  original  group 

(S2=  x x (26b) 


there  manifestly  comes  the  following  group 

(SS  =-  x xx  x (26c) 

That  is:  ue  obtain  a group  of  3 elements  (121).  The  corresponding 
directional  characteristic  Ql,  2,  l]  is  then,  according  to  the 
theorem  of  Bridge  equal  to  cos2  9 0. 

Correspondingly,  it  follows  that 

is  = U,  1|[1,  nil.  l]  — U,  3, 3 l]  — co839>,  (26d) 


<S  = U.  (I').  (S)  ...  (2)]  = cos»9>.  (26e) 


In  order  to  find  the  symbolic  product 

• 'Hs  ^ . . . <»,.,]  [i,  , V • • bn]  ( 26f ) 

i.e.,  in  order  to  find  the  directional  characteristic  of  the  straight 
line  group  which  results  if  each  individual  radiator  of  the  one 
gx'oup  is  replaced  by  an  equal  and  similarly  directed  group  (so  that 
the  original  radiator  and  all  radiators  appearing  in,  its  place  lie 
on  the  same  straight  line),  we  form  the  following  rectangular  scheme: 
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It  is  then  easily  perceived  that  for  the  determination  of  3?„  ■ sJv6  f 
the  sums  must  be  formed  of  the  products  contained  on  the  diagonals. 

It  then  turns  out  that 


• %,  [flj  V (i.£ a^b j L c.'j6.i  -f-  Oj63,  (26g) 

l _ i “t  (tn , _ i bn , an  6„] , 

This  result  may  also  be  explained  in  a simple  manner  by  con- 
sidering the  product 


(°i  •*  -i-  «a  **  H- ' ■ H a„,  x"‘)  <6i  * + b»  *a  + ■ • • + br  a:")  ( 26h) 

- 6,  a:2  -|-  (a.,  6,  + bs)  X'3  -j |-  am  bn  a;”1  ■'  n 

Of  interest  is  the  converse  of  the  question:  With  a hypothetical 
straight  line  group  with  a fixed  number  and  distance  for  the  receivers, 
how  can  the  directional  characteristic  be  changed  by  changing  the  in- 
dividual velocity  amplitudes,  i.e.,  by  changing  the  sensitivity. 

Thus  with  six  radiators  which  are  arranged  at  a uniform  distance  d 
on  a straight  line  the  following  seven  directional  characteristics 
are  obtained: 
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The  corresponding  directional  characteristics  are  represented 
in  ?ig,  lo.  According  to  the  moment  of  inertia  theorem  it  is  immed- 
iately clear  that  the  directivity  of  9?,  is  greatest.  Since  the 
sensitivity  of  the  receivers  becomes  more  and  more  concentrated  near 


Fig.  10.  Directivity  of  the 
straight  line  group  of  six 
receivers  with  different 
sensitivity. 


Fig.  11.  Secondary  maxims 
of  the  straight  line  group 
of  six  receivers  with 
different  sensitivity. 


the  center,  the  directivity  must  then  of  necessity  decrease  more  and 
more.  Fig.  10  shows  this.  It  appears  of  considerable  importance 
that  an  opposite  trend  exists  for  the  sizes  of  the  secondary  maxima. 
The  part  of  the  curve  corresponding  to  <p > i is  greatly  enlarged 
in  Fig.  11  in  oixior  to  show  the  size  of  the  secondary  maxima  and 


it  is  evident  that,  with  regard  to  the  size  of  the  secondary  maxim,  9^ 
is  very  unfavorable.  Thus,  e.g.,  for  «,  sizes  of  the  secondary 
ma-srfma  remain  below  1$  as  compared  with  22$  for  9^  . ^ith  a 

fixed  arrangement  of  radiators  we  can  thus  influence  the  directivity 
or  the  sizes  of  the  secondary  maxima  in  the  most  favorable  sense  by 
adjusting  the  sensitivity  (i.e.,  the  amplification)  as  is  necessary. 

In  order  to  derive  the  directional  characteristic  for  the  rec- 
tangular group,  we  proceed  from  the  general  directional  characteristic 
of  the  straight  line  group.  Here  we  substitute  cos  £ for  sin  y 
in  Formula  (19)  in  order  to  show  that  the  relation  hold3  independently 
of  the  measuring  plane.  It  was  necessary  for  Formula  (19)  that  the 
measuring  plane  coincide  with  the  ZY-plane.  We  then  obtain 


sm  — — cos  j i 

to  = L A 

M . I nd 


jism  ^ ' C0B/ 


If  we  now  Imagine  each  individual  radiator  replaced  by  a group 
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Fig,  12.  Rectangular  group. 


Fig.  13.  Rectangular  piston 
membrane. 


of  m radiators  parallel  to  the  X-axis  (Fig.  12),  according  to 
the  rule  given  above  the  directional  characteristic  of  this  rectan- 
gular group  then  turns  out  to  be: 


. IniTidi  . Im/k/j  o 

HIU  l - ^ CO Brt  Bin  J ^ OOB/J 

. |n(/i  I . nd a . 
m »in  ^ 1 eon*  * » Bin  ^ co h(S 


If  we  allow  (2,  and  <ia  to  become  smaller  and  smaller  and  m 
and  n to  increase  arbitrarily  so  that,  in  the  limit,  wrf,  -*a  and 
nd^b  (Fig.  13),  we  then  obtain  the  general  directional  charac- 
teristic for  the  rectangular  piston  membrane: 
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Here,  it  must  be  supposed  that  the  membrane  oscillates  cophas- 
ally  on  both  sides  of  the  ^Y-plane  or  that  it  oscillates  in  au  in- 
finite rigid  wall  as  a simple  piston  membrane. 


In  place  of  the  densely  filled  circular  line,  a circular  group 
is  frequently  used  in  practice  which  consists  of  a definite  number 
n of  radiators  which  are  arranged  equidistantly  on  a circle 
of  diameter  d ^ If  we  choose  n as  an  even  number  (n  — 2 m), 
then  by  summing  up  the  response  of  each  pair  of  diametrically 
situated  radiators,  the  following  relation  for  the  directional 
characteristic  may  be  easily  derived: 
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= 1 V, 
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(31) 


Here  the  field  point  line  is  determined  by  <p  and  y (see 
Fig.  14) • For  the  larger  values  of  m , in  particular,  calcu- 
lation by  (31)  is  very  time  consuming.  By  the  use  of  Bessel's 
functions,  the  sum  may  be  transformed  into  an  infinite  series 
which  is  considerably  more  convenient  for  the  calculation.  It 
then  turns  out  that 


00 

— J0  ( siny)  -|-  2 ^ ^ sinyj  • coaipmtp.  (32) 

p 1 

The  practical  significance  of  thi3  relation  which,  at  first, 
appears  complicated,  is  immediately  recognized  if  it  is  considered 
that  the  first  terms  in  the  sum  Jt  , Jim  etc.,  very  rapidly 
assume  a practically  negligible  value  so  that  it  is  generally  suf- 
ficient to  consider  only  the  first  term  Jim  . 

v 

Moreover,  with  the  aid  of  (32)  the  frequently  important  ques- 
tion can  be  decided  as  to  how  densely  a circular  line  must  be  cov- 
ered by  the  radiators  in  order  that  the  directional  characteristic 
may  be  indistinguishable  from  that  of  the  continuously  covered 
circular  lino.  If  we  choose,  say,  four  receivers,  then  the  decisive 


Rlectr.  .‘hchr,  I'och.  Vol.  6 (1929),  p.  170.  (or  the  NRL 
translation  (,/-ll<i)  of  this  paper). 
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correction  term  is  smaller  In  absolute  value  than  1 2J4 (—--sin y)j . 

If  we  denote  the  part  from  y = o(J0=l)  to  y = yl  Where  y,  corre- 
sponds to  the  first  zero  position)  as  the  principal  part  of  the  charac- 


teristic of  = jJ— sinyj 


Fig.  14.  Calculation  of  the 
directional  characteristic 
of  the  circular  group. 


then  yx  is  determined  by  ~ 3in  ^ _ 2)4 

and  since  j4(x)  < 0,06  for  x < 2,4, 

the  principal  part  of  the  characteristic 
is  then  not  changed.  That  is,  the  direc- 
tivity of  a circular  group  of  four  radi- 
ators is  completely  equivalent  to  the 
directivity  of  an  arrangement  which  is 
densely  covered  with  radiators  regardless 
of  how  great  a radius  is  chosen  for  the 
arrangement.  It  can  likewise  be  concluded 
that  the  complete  behavior  of  the  direc- 
tional characteristic  of  the  circular 
group  can  be  represented  with  sufficient 

accuracy  by  ,/0  ^ sm  yj  if  the  number  n 

is  so  large  that  the  condition 


71  ^ 


(33) 


is  fulfilled. 

Since  the  circumference  of  the  circle  is  * d --- » •«  where  « is 
the  distance  measured  on  the  arc,  we  can  then  so  formulate  the  con- 
dition (33)  that  the  distance  measured  on  the  circle  between  two 
neighboring  radiators  must  be  somewhat  smaller  than  /.  (more  pre- 
cisely (a/1  < l - 2/m)  ) in  order  that  the  circular  group  directional 

characteristic  may  be  given  by  »"'y) . 


Closely  connected  with  the  beam  sharpness  of  an  arrangement  is 
the  matter  of  the  resolving  power  (Trennscharfe) . Then  it  is  a 
question  of  when  a receiving  arrangement  is  in  the  position  to  per- 
ceive separately  two  sound  sources  placed  at  a great  distance  away 
and  at  a small  distance  from  each  other.  If  we  imagine  a sound 
source  placed  synmetrically  on  each  side  of  the  perpendicular  bisec- 
tor of  a receiver  arrangement  (sources  1 and  2)  and  at  a great  dis- 
tance r and  if  we  imagine  furthermore  a circular  surface  as  the 
receiver  arrangement  (Fig.  1$),  being  rotated  about  the  X-axis,  then 
the  principal  maxima  produced  by  each  radiator  will  be  perceived 
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separately  if  the  distance  a between  the  radiators  is  sufficiently 

great.  If  one  assumes  with  Rayleigh  that 
the  maxima  can  then  be  separated  if  the 
maximum  of  one  coincides  with  the  minimum 
of  the  other,  the  necessary  condition 


fig.  15.  Resolving  power 
of  the  circular  pie  ton 
membrane. 


7ld 

-j-  sin  <p'^  3,83 


(3» 


then  follows  from  the  fact  that  the  equa- 
tion 


2J1 


nd 


=■■■  0 


- sin  <p 


(33b) 


i3  fulfilled  for  ^-sin<p-=3,83  Since 


sin  (p<x><p  — ajr 


(33c) 


it  follows  that 


djX  > 1,22  r/o 


(34) 


In  order,  therefore,  to  be  able  to  separate  with  an  acoustical 
objective  radiators  placed  at  a great  distance  r from  the  objective 
and  at  a small  distance  a from  each  other,  the  diameter  of  the  ob- 
jective measured  in  wave  lengths  must  at  least  be  equal  to  r/o. 

The  preceding  considerations  referred  exclusively  to  the  case 
where  all  parts  of  the  radiating  membrane  or  radiating  system  moved 
with  the  same  amplitude  and  phase.  In  practice  this  is  not  by  any 
means  the  case.  In  general,  the  force  exerted  on  a membrane  is  not 
uniformly  distributed  over  the  whole  surface  but  is  exerted  at  the 
midpoint  or  along  a line.  Since  there  is  no  absolutely  rigid  mem- 
brane (above  all  not  if,  in  addition  to  rigidity,  as  small  a weight 
as  possible  is  necessary  to  produce  a good  efficiency  and  frequency 
response),  the  force  acting,  say,  at  the  midpoint  will,  due  to  the 
(damped)  propagation  with  finite  velocity  of  the  tranverse  elastic 
waves,  be  able  to  act  neither  in  phase  nor  with  uniform  amplitude 
on,  say,  the  outlying  elements  as  soon  as  the  dimensions  of  the 
membrane  are  no  longer  small  compared  to  the  wave  length  of  this 
transverse  wave.  Furthermore,  the  membrane  is  generally  retarded 
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in  its  motion  by  its  support  at  the  boundary.  In  principle  making 

allowance  for  a different  phase,  in  our 
formulae  will  cause  no  difficulty,  we 
would  then  need  to  use  the  given  velocity 
amplitude  w(x,  y)  as  a complex  quantity. 

For  simplicity  we  will  confine  ourselves 
to  a variable  amplitude  and  accordingly 
assume  w(x,  y)  as  a real  function. 

As  the  simplest  example,  we  consider 
the  directional  characteristic  in  the 
Zl-plane  of  the  rectangular  membrane 
represented  in  Fig.  16  where  the  membrane 
is  shown  at  the  moment  of  greatest  de- 
flection. 

Here  w is  to  be  given  by  w(y)  - 1 — y2/b2 
We  then  find  that 


fig.  16.  Vibration  form 
of  the  rectangular 

membrane  »<*,»>- 1 - «*/»*. 
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(34a) 


so  that 


=-•  5 • 


(34b) 


and 


-ft 

05  3 /sinu  \ 

^ ^ — C08“)’ 


(35) 


where 


2 nb  . 

u = -y-  sin  y 


(35a) 


In  order  to  find  the  directivity  coefficient,  we  have  to  calcu- 
late the  value  of  u for  which  -=  0 .707 . It  then  results  that 
u~  1.8.  The  directivity  coefficient  is  then  equal  to  1,8  In  = 33° 
and  the  half- value  beam  width  y is  given  by  q?  = 33^  a/2  b . Using 
the  approximation  formula,  we  would  find  that 


(35b) 
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-w  •-  “■-  •-  ■““*  > 1 *" 


and  hence  that 


?>  = 33,5° -A/2b. 


(35c) 


As  a second  example  for  the  rectangular  membrane  we  choose 
(Fig.  17): 


w(x,  y)  — 1 — 2 f/8/&2  (35d) 

and  find  K,  ~ i- 

For  the  directional  characteristic  in  the  ZX-plano,  it  turns 
out  that 

9i  = .I4(^_C08u)_3.?! (35e) 

u*  [ u j U 


V.'e  here  determined  the  direc- 
tivity coefficient  from 

_ - cosu)  - 3 — - - 0,707 , t36) 

and  it  hence  follows  that  u = 4,4 . 
The  directivity  coefficient  is 
thus  4,4/71  = 80°  while  the  half- 
value beam  width  is  <P  = 80M/25 


Fig,  17.  Vibration  form  of  the 
rectangular  membrane  w<*,w=i-w. 


nere  something  new  occurs  since  SR  does  not  possess  a maximum 
but  a minimum  on  the  Z-axls.  a series  development  of  a?  yields: 

*-i+5- 

Fig.  18  shows  a comparison  of  the  directional  characteristics 
of  the  rectangular  piston  membrane  of  length  2b  and  the  two  latter 
membranes.  One  recognizes  how  substantial  is  the  change  if  a mem- 
brane with  nodal  lines  is  operative  in  place  of  the  desired  piston 
membrane.  On  the  other  hand,  this  knowledge  can  be  turned  to  good 
account  if,  with  a large  membrane  (whioh  is  necessary  in  order  to 
provide  the  required  energy),  it  is  a question  of  avoiding  a sharp 
convergence. 

As  the  third  example  we  choose  a tightly  stretched  circular  mem- 
brane. Here  we  can  generally  carry  out  the  calculation  for  the  curves 


^ 3«e  Note  7. 
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given  by  Mq)  — (l  - o*>*)*  . (Fig,  19  represents  the  curve 

w(e)  = (l  — S*l1ri)n  for  n = 0,  2,  4,  and  8.)  Ihus  we  find^ 


w-  = ,ri“i  und:  »» = 2"  • ‘ • (» -;■  l) ! %4l) , 


(37) 


wherein  JM+1(u)  is  the  Bessel’s  function  of  the  order  (»-fl)  and 

2.-r r . 

m = ^ sm  y • 


Fig.  18.  1.  Directional  characteristic  of  the  rec- 
tangular piston  membrane  in  the  symnetry  plane. 

2.  Directional  characteristic  of  the  rectangular 
stretched  membrane  in  the  symmetry  plane.  Di- 
rectional characteristic  of  the  rectangular  mem- 
brane with  nodal  lines  in  the  symmetry  plane. 


From  the  approximation  formula  calculated  with  the  moment  of 
inertia: 


r£  7ifc  siny 
n 4-  2 


%,  1 

we  find  the  half -value  beam  width  to  be 

(p  15°  • )'«  -|- 1 • A/r . 


(37a) 

(37b) 


The  complete  behavior  of  directional  characteristics  corresponding 
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Fig.  19.  Vibration  form#  «■(<■)  u-  -•*. /■*)** 
for  >1  - » . t , 4 . « . 
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to  the  vibration  forms  of  Fig.  19  is  represented  in  Fig.  20  for 
n = 0,  2,  ki  8.  It  is  to  be  mentioned,  however,  that  by  a linear 


Fig.  20.  Directional  characteristic  of  the 
circular  membrane  with  the  vibration 
form  rtji-a-iwfor  n~o, 2,4,8. 

combination  of  the  vibration  forms  (w  = (1  -e2/r2)"  3X1  arbitrary  ro- 
tat tonally  symmetric  vibration  form  of  the  membrane  (even  with  nodal 
lines)  can  be  represented  to  an  arbitrary  degree  of  approximation. 


„e  can  state  the  directional  characteristic  corresponding  to  a 
general  vibration  of  the  membrane 


u>(q)  :=  <*o  + U - QVr*)  + aa(!  — etlr*)i  + ■ + 0,1(1  — (38) 


It  is 


*-—i— 

a°  1’  2°J  + 3°s  + 


+ ] | °« 
n + 1 


2a, 


J,(u) 


+ 22 • 1 ! Ox",-  4-  • - 

+ + 


(39) 


An  exceptional  case  arises  if  the  expression  in  the  denominator 


®i  + • 


+ »+  la" 


(39a) 


vanishes.  Then  formula  (39)  for  9i  becomes  unusable.  However,  the 
difficulty  may  be  immediately  reeved  if  the  pressure  amplitude  is 
calculated  directly  by  equation  (11).  It  then  follows  tiiat 

v -u  Yd},'  (39b) 

r 
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and  for  y — o 


v = TTi{jw^>  y)dF  = 0 

F 


(39c) 


This  means  therefore  that  since  wm  ~ 0 the  pressure  amplitude 
at  a sufficient  distance  on  the  normal  axis  vanishes.  In  addition,  it 


Fig.  21.  Circular  membrane  with 
the  deformation  volume  zero. 


Fig.  22.  Relative  Bound  pressure 
amplitude  of  a circular  membrane 
with  the  deformation  volume  zero. 


is  clear  that  the  directional  characteristic  (which  represents  the 
ratio  of  the  pressure  amplitude  in  an  arbitrary  direction  to  the 
pressure  amplitude  on  the  nonral  axi,s)  loses  its  meaning.  In  other 
respects,  the  calculation  offers  no  difficulty. 

A suitable  example  is  given  by  w — l - 2 ei/ri  . Therefore  we 
have  a0  •.=  - l,  = + 2,  «#  = a3 . . . = o (see  Fig.  21) . 

Here  wm  - o and  p is  found  from 

*W«)  %M)  (40) 


where 


n<i  , 
x may, 


'H0(u)  - - 


2 J>(n) 


02  . O 


J8(n) 


(40a) 


For  small  values  of  u , it  is  recognized  from  tho  series  devel- 
opment 


9M«0  ■-«,(#) 


2*'.  0!  31 


• u‘  w‘ 

2*-  1141  + 2"  t 21 5! 


(41) 


that  even  for  directions  in  the  neighborhood  of  the  Z-axis  the  effect 
of  the  membrane  is  extraordinarily  small.  The  complete  behavior  of 
-JKo(it)  is  represented  in  Fig.  22. 
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(b)  The  directional  characteristic  for  a noise. 


In  practice,  of  the  sound  sources  which  are  to  be  measured 
(e,g.,  underwater),  two  essentially  different  types  are  to  be  dis- 
tinguished. Either  it  is  a question  of  artificial  sound  sources 
which  are  to  send  out  pure  tones  (light  ship  transmitters,  bells, 
signal  senders  on  ships)  or  of  more  or  less  natural  sound  sources 
which  possess  a noise-like  character  and  vdiich  are  generally  un- 
intended and  undesired  (propeller  noise,  machine  noise). 

Therefore  an  investigation  to  determine  how  the  preceding  con- 
siderations may  be  applied  in  such  cases  suggests  itself.  Here  we 
will  assume  the  case  occurring  most  frequently  in  practice  where  a 
noise  source  is  present  from  whose  continuous  spectrum  a definite 
frequency  range  (bounded  above  and  below)  is  received.  As  is  well 
known,  this  may  be  accomplished  without  difficulty  by  an  electrical 
filter.  It  is  additionally  assumed  that  the  receiver  receives  all 
frequencies  of  the  range  in  question  with  equal  intensity  in  the 
principal  direction  - i.e.,  that  the  transmitter  favors  no  frequency 
and  also  that  no  frequency  dependency  exists  on  the  path  from  the 
transmitter  to  the  receiver  in  the  medium  (absorption,  reflection) . 
We  will,  with  advantage,  define  the  quadratic  mean  value 


(42) 


as  the  directional  characteristic  where  and  na  are  the  limits  of 
the  range  and  9?  is  the  previously  defined  directional  characteristic 
dependent  on  the  frequency  n. 


If  we  choose  as  the  simplest  measuring  apparatus  two  nondirec- 
tional  receivers  separated  by  the  distance  d (which  is  small  com- 
pared to  the  wave  lengths  in  question) , then 


ffl,- 


(42a) 


Carrying  out  the  integration; 


sin  (*,  — *,)  , , , 


(43) 


wherein  for  abbreviation 


(43a) 


7Ul 

nd  . 

7i,  - sin  y — 
1 c 

j-  am  y = Xj 

7id  . 

nd  . 

n2  — sm  y 

siny  = x„ 
a2 

One  sees  that  for  x2  = xx  formula  (43)  goe3  over  into  the  original 
directional  characteristic  3^  in  (14).  It  is 

= *,  = * = Vi  + l COS 2 X - coax . (43b) 


We  will  initially  investigate  the  behavior  of  the  directional 
characteristic  when  the  pass  band  is  precisely  one  octave.  Here  we 
set  , 

x2  — 2 x,  xx  = x = -jj  sin  y 


and  obtain: 

a,  , + 


(44) 


From  Fig,  23  we  see  that  the  directional  characteristic  is  now 
substantially  different  from  the  previous  one  (for  a single  tone,  Fig. 
7,  curve  l).  Instead  of  the  zero  and  unity  value 3 appearing  period- 
ically with  increasing  x , there  appears  here  only  one  principal 


Fig.  23.  Directional  characteristic  of  a 
straight  line  group  (Length  a)  for  octavo 
reception: 

1.  Two  receivers,  2.  Three  receivers, 
3.  Four  receivers,  4.  Six  receivers, 

5.  Densely  covered. 
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maximum  with  the  value  unity,  rind  for  the  larger  values  of  x , the 
curve,  oscillating,  always  approaches  the  value  fi  . In  order 
to  investigate  the  influence  of  the  number  of  receivers  for  octave 
reception,  the  cases  are  represented  in  Fig.  23  when  the  receiver 
array  consists  of  2,  3,  4,  6 and  very  many  receivers  - the  total 
length  of  the  array  being  kept  constant  and  equal  to  a . One 
sees  how  the  directivity  decreases  here  also  with  an  increasing 
number  of  receivers.  Simultaneously  a steady  decrease  occurs  in 
the  value  which  the  curve  approaches  with  increasing  “f! 8kt  y . 

The  values  corresponding  to  the  receiver  ni’mbers  2,  3,  43  and  6 
are  respectively  and  )/j 

while  the  value  zero  corresponds  to  a large  number  of  receivers. 

In  general,  for  a straight  line  group  of  n receivers  at  the 
same  distance  d and  for  a frequency  range  from  v to  v + pv  it 
is  found  that 

®>  =]/  -n  + li  (»  - n)  cos [(7*  + 2 )»»*].  (45) 


In  order  to  learn  the  influence  of  the  size  of  the  pass  band, 
the  cases  n » 6, p = Oj  0.2j  0.5>  1}  3 are  represented  in  Fig.  24. 


Fir.  24,  Influence  of  the  bend  width  on  the 
directional  characteristic  of  a straight  line 
group  of  six  receivers. 
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It  is  seen  that  the  size  of  the  pass  band  determines  the  amount  of 
the  variation  above  and  that  the  greater  the  pass  band, 

the  smaller  this  variation.  In  Fig.  24, 

Curve  1 corresponds  to  the  case  n = 6 p = 3 

Curve  2 corresponds  to  the  case  n = 6 p = 1 

Curve  3 corresponds  to  the  case  n = 6 p = 0.5 

Curve  4 corresponds  to  the  case  n - 6 p = 0,2 

Curve  5 corresponds  to  the  case  n = 6 p = 0 . 

d 

The  values  of  j sin y are  read  on  the  abscissa  axis,  a = 5d 
is  the  length  of  the  base. 


As  a further  erK&mple,  we  consider  the  noise  reception  for  a 
dense  receiver  array  on  the  circular  line  and  on  the  circular  sur- 
face. The  corresponding  directional  characteristics  are 


r, 


(46) 


where  again 


*1 


(i 

81U  )' 


and 


7t<l  . 

SHI  y . 
A, 


For  octave  reception  (rs 


: 


(47) 


(47a) 


(46) 


This  is  represented  in  Fig.  25.  For  comparison,  the  ordinary 
directional  characteristics  corresponding  to  the  upper  and  lower 
limits  of  the  octave  are  also  inserted  (curves  1 and  3)* 
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Pig.  25.  Directional  characteristic  of  the  densely  covered  circular  line. 

1.  For  the  fundamental  (wave  length  \ ).  2.  For  the  octave  (wave 

length  ^ to  \/z).  3.  For  the  harmonic  (wave  length  x/z). 

(c)  With  artificial  compensation 

When  the  position  of  a sound  source  is  to  be  determined  with  a 
given  receiver  array,  this  can  be  done  by  rotating  the  array  until 
the  maximum  intensity  is  attained  on  the  indicating  instrument  (or 
for  the  ear) . When  the  elements  of  the  array  lie  in  one  plane,  this 
maximum  response  is  obtained  when  the  plane  of  the  array  is  perpen- 
dicular to  th8  sound  direction* 

However,  the  same  result,  which  depends  on  a coincidence  of 
phase,  can  also  be  attained  by  artificial  compensation  without 
rotating  the  array.  This  is  done  by  variable  electrical  delay 
circuits  which  are  inserted  between  the  fixed  receiver  array  and 
the  indicating  instrument,  ^ach  position  of  the  compensator,  which 
controls  the  delay  time,  corresponds  to  a definite  direction  in 
space  (in  the  measuring  plane)  for  which  the  receiver  array  is  in 
phase.  In  order  to  find  the  direction  of  the  sound  source,  leaving 
the  array  fixed,  one  now  lias  only  to  turn  the  compensator  and  read 
off  the  corresponding  beam  angle. 

Artificial  compensation  by  electrical  delay  circuits  is  quite 
important  for  directive  reception.  Not  in  the  least  contributive 
to  its  usefulness  is  the  simplicity  and  accuracy  of  the  operation 
of  the  electrical  circuits,  ’l'hese  generally  consist  of  a number  of 
uniformly  constructed  sections  which  contain  solf-inductanoe  and 
capacitance  and  are  so  assembled  that  the  self-inductance  coils  are 


Using  the  inciioated  argument  in  this  figure,  the  curves  for 
the  fundamental  and  the  harmonic  will  be  identical.  If,  however, 
sin  Y used  the  argument,  curves  similar  to  curves  (l)  and 
(3)  will  ce  obtained  but  curve  (l)  will  correspond  to  the  harmonic 
while  curve  (3)  will  correspond  to  the  1‘undamontal  froquenoy. 


inserted  in  series  while  the  condensers  are  in  parallel.  If  o is 
the  value  of  the  self -inductance  and  (s  the  capacitance  of  the  con- 
denser, the  delay  furnished  per  section  of  the  circuit  is  yg  • <$. 


Furthermore,  this  delay  can  be  made  quite  independent  of  the  frequency 
if  one  arranges  it  so  that  the  limiting  frequency  defined  by 

i i 

n = — . - 

;i  i a -is 

is  sufficiently  far  above  the  frequency  range  to  be  passed.  In  this 
way  one  succeeds  in  imparting  a pure  time  delay  per  section  to  a 
noise  signal  without  changing  its  character. 

If  we  have  an  arbitrary  receiver  array  arranged  in  space  then 
the  compensation  apparatus  can  be  obtained  in  a very  simple  manner 
for  a given  measuring  plane  by  a single  delay  circuit.  If  we  imagine 
the  receivers  projected  on  the  measuring  plane  (XX-plane)  and  that 
their  coordinates  are  given  by  (*1,2/1)  (**,!/*)...(*„ yn)  it  is  then  clear 
that  for  a compensation  in  the  measuring  plane  only  these  projections 
are  of  concern,  otherwise  expressed,  the  individual  receivers  can 
be  arbitrarily  displaced  perpendicularly  to  the  measuring  plane  with- 
out any  change  in  the  compensation  action.  If,  for  simplicity,  we 
confine  ourselves  to  three  receivers,  the  natural  directional  char- 
acteristic is  then  given  by 

ft 

J V*  cik{jrn*n*  a + j/aCoh/0^  (49) 

M “ 1 

where  a and  p are  the  direction  angles  cf  the  field  point  line.  If 
we  let  the  compensation  direction  in  the  measuring  plane  have  the 
direction  angles  <i‘  and  y , then,  upon  applying  the  delays,  the  arti- 
ficial directional  characteristic  is  given  by 


ft 

- - J %'  - co*y)  ■ ->  co#v*))t 


(50) 


If  we  allow  the  sound  source,  i.e.,  the  field  point  line  (a,/)), 
to  shift  so  that  it  mates  a complete  circuit,  W*  assumes  its  greatest 
value  when  the  direction  of  the  sound  source  coincides  with  the  com- 
pensation direction.  Exactly  the  saim  circumstances  occur,  however, 
If,  with  a fixed  sound  source,  we  allow  the  compensation  line  (y,  y) 
to  shift. 

If  we  regard  the  configuration  of  Fig.  26  as  a scale  drawing  of 
the  receiver  array  which  can  be  rotated  about  the  origin  0,  then  in 
the  rotation  the  correct  retardations  *„co8y-t -yni  *\p  are  attained 
simply  as  the  projections  on  a fixed  line  which  we  will  assume  to 
be  the  X-axis.  If  we  furthermore  imagine  perpendicular  contact  bars 
which  are  connected  with  the  individual  sections  of  the  delay  circuit 
(as  is  shown  in  Fig.  26)  and  place  sliding  contacts  at  (*ii/t)  (*ayj).  !■*»&) 
which  are  constantly  connected  with  the  one  terminal  of  the  corres- 
ponding receiver  while  the  other  terminals  are  connected  with  a 
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Fig.  26.  Compensation  apparatus  for  a group  of 
radiators  (l,  2,  3)  (small  in  comparison 
to  the  wave  length)  with  an  arbitrary 
arrangement  in  space. 


common  return  conductor  at  the  circuit  input,  then  each  impulse  from 
a receiver  reaches  the  input  of  the  amplifier  with  the  desired  delay. 
At  the  same  time,  by  an  appropriate  design  of  the  delay  sections  the 
time  delays  of  the  circuit  are  to  be  made  to  correspond  to  the  ve- 
locity of  sound. 


For  the  uncompensated  group  the  directional  characteristic  is 
independent  of  the  position  of  the  object  to  be  located.  This  is 
not  the  case,  however,  with  the  compensated  group.  Here  the  di- 
rectivity will  be  so  much  the  more  dependent  upon  the  angle  the 
mo i3  the  pattern  of  the  projections  of  the  receivers  on  the  measur- 
ing plane  departs  from  that  of  a circular  array.  If  we  consider 
the  straight  line  group  of  Fig,  g}  the  directional  characteristic  SR* 
for  the  compensation  direction  characterized  by  y0  is  given  by 


nnd  , . . 

- (»m  Y ~ »m  Vo) 


- Hill  y0) 


(51) 


In  Ffg.  27  the  directional  characteristics  are  represented  in 
polar  coordinates  for  a straight  line  group  («  = o,  d - A/2)  for  the  com- 
pensation directions  0°,  45° > 60 °,  and  90° j here  the  secondary  max- 
ima are  left  out.  It  is  noteworthy  here  that  the  directivity  is 
changed  only  slightly  at  ~ 45°  . However,  for  larger  angles  a 
decided  broadening  of  the  principal  maximum  arises.  Furthermore, 
the  principal  maximum  no  longer  lio3  symmetric  to  the  compensation 
direction.  This  las  the  consequence  that  too  great  a value  for  y0 
is  found  if  the  measurement  is  undertaken  in  the  usual  manner  wherein 
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the  compensator  positions  for  equal  loudness  on  both  sides  of  the 
maximum  are  determined  by  ear  and  the  mean  of  these  positions  is 
taken. 


fig.  27.  The  directional  characteristic  of  a straight 
line  group  of  six  radiators  with  compensation. 

)',-(>•,  2,  y,~4f>\  3.  y,  - 00* , 4.  »-00*. 

These  disadvantages  of  the  straight  line  group  will  be  avoided 
If  the  densely  covered  circular  group  is  used  so  that  the  measuring 
plane  coincides  with  the  receiver  plane . If  we  next  calculate  the 
directional  characteristic  of  the  circular  group  compensated  for  an 
arbitrary  angle  »0.  /?„,  y0  we  find  that 


y{.  . . I fill'lco*  * - - 1/ (COH //  - 

*2 nr)  ' 

If  we  let  x = r cow  f and  y —■  r Kin  f and  set 


(2 


vmj\  ■ ■ ei)«/J0 
COH  V ~ C08  \0  ’ 


(5] 


it  then  follows  that 


and 


__  JL_  j gtirtco89>  (co8a  — cosot)  + 3in?j(oo8^- cosft)],^^ 
0 

2.T 

gj^  __  JL_  ^gi*r\(co8a--cosa9),+  {co8/f  — cos/fo)*^9^  “ ? 

6 

% — J0(¥  V(cosot  — cosa0)2  + (cos/S  — cos/J0)2) 


(52) 


For  ap  = 90,  jS0  — 90 c,  SR*  must  become  the  uncompensated  directional  char- 
acteristic as  given  by  (15);  if  we  let  cos  «0  - cos  0O  = o in  (52)  it 
then  follows  that 

9t*  - Jo(~r  sill  /),  (52a) 


since 


]'  (cos  a — cosa0)2  ' - (cos/?  — c(wji0)2  = ]/  cos2  jc  -P  OuSs/S 


= y 1 — cossy  = cinyt 


(52b) 


If  the  sound  source  is  on  the  X-axis  and  the 
XZ-plane  is  the  measuring  plane,  then,  since 
«o  = 0,  A,=  !)0o,  /?  = 90°  (see  *ig.  2S): 


94  = J0  (”/  -2a>n2|-j 


(53) 


If  the  souM  source  is  on  the  X-axis  and  the 
XI -plane  (receiver  plane)  is  the  measuring 
plane,  then,  since  «0  — 0,  — 90°,  <x + 0-90°: 


34  » 1 


0 (x  2 8itl'2  ) ' 


(54) 


In  practice,  this  last  case,  in  which  the 
measuring  plane  and  the  receiver  plane  coincide, 

on  TU^  „„„4.  • is  of  -.'.articular  significance.  From  formula 

Fig.  28.  Different  po.i-  ( } r0co^izo  that  the  directivity  of  tte 

t oe.  of  the  maaBurlog  ^mpematod  group  agrees  «dth  the  directivity 
plane  for  the  comp.n-  of  ^ 

uncompensated  group  (see  formula  (15)) 

•a  ed  circular  group.  go  m aga^  have  to  choose  the  angle  20° 

as  the  directivity  coefficient.  If  a definite  directivity  is  required, 
the  ratio  of  the  circle  diameter  to  the  wave  length  is  thus  determined. 
If  wo  choose  a half  value  bean  width  of  13°  as  an  example  then  from  (18) 

" W =1,5. 


13  ' 


The  remaining  question  is:  how  many  radiators  are  necossaxy  in  order 
that  no  practical  difference  may  exist  between  the  resulting  direc- 
tional characteristic  and  the  directional  characteristic  of  the 
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densely  covered  circular  group.  If,  for  simplification,  we  allow  the 
compensation  line  to  go  through  a radiator,  we  then  find  an  analogous 
representation  by  BQ8sel‘s  functions^  '• 


% = Ja(~sm^j  + Jp  n sin  j cos  ~ 


np 


= 1 


(55) 


From  this  it  also  follows  that  the  directivity  of  the  compensated 
circular  group  agrees  practically  completely  with  the  directivity  of 
the  group  consisting  of  four  radiators.  .If  the  complete  characteristic 
is  to  agree,  the  inequality 


n — ~J~  i"  2 


(55ft) 


must  be  fulfilled.  Or,  the  distance  a of  two  neighboring  radiators 
must  be  somewhat  smaller  than  X/s  . lore  precisely 


i 

n 


(55b) 


In  other  respects,  the  formula  provides  a substantially  more  simple 
calculation  of  the  directional  characteristic. 


As  an  example  m calculate  the  directional  characteristic  for 
w = 6,  d = iX  and  find  the  following  tables  by  the  approximation  formula 

% = (4,7 1 sin  a/2)  + 2 J6  (4,7 1 sjj  a/2)  cos  3 « . ( 55c ) 


<x 

4,71  sin  «/2  ! 

J,{  4,71  »>n  a/2) 

2 «/# (4,71  shi oc/2)  cos3«  ! 

91* 

0 

0,00 

1 

i 

0 

1 

10  , 

0,42 

0,06 

0 

0,96 

20 

0,82 

0,84 

0 

0,84 

30  1 

1,22 

0,06 

1 0 

0,66 

40 

i 1,00 

0,45 

0 

0,45 

50 

1,09 

0,23 

0 

0,23 

60 

1 2.36 

0,02 

o 

0,02 

70 

! 2,70 

-0,14 

0 

-0,14 

80 

i 3,03 

-0,27 

-0,01 

-0,28 

00 

3,33 

-0,35 

0 

— 0,35 

100 

3,61 

0,30 

4 0,03 

-0,30 

no 

3,86 

- 0,40 

. 1-0,07 

-0,33 

120 

4,08 

•0,30 

4-0,10 

-0,20 

130 

4,27 

-0.37 

1 +0,12 

-0,25 

140 

4,43 

-0,34 

j 4-0,08 

- 0,26 

150 

4,55 

0.31 

+0,00 

-0,31 

100 

4.01 

-0,20 

-0,10 

-0,30 

170 

4,60 

- 0,27 

0,18  | 

- 0,45 

180 

4,71 

- 0.27 

- 0.21 

-11,48 

**  dlektr,  Nochr  ,-iechn. , Vol.  6 (1929),  p.  17G,  (or  the  NRL 
translation  (/;il4 ) of  this  paper). 

* i'he  series  called  for  by  the  summation  sign  should  be 
multiplied  by  the  faotor  2. 
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while  the  direct  calculation  by  the  formula 

9J,.  = J{cos[£ji(1  — cosa)]  + cos[j7t(cos(i3c  + 60°)  — cos60°)] 
+ cos[f  jz(cos(a  + 120°)  — cosl20°)] 


(56) 


yields  the  following  table 


a. 

•Pi  = 

1 — COS  <X 

7'j  “ 

COS(a-r60'5) 
— COS 60 3 

Ts  “ 

COS  (a  + 1‘203) 
— cos  120° 

TOSf 

3 

cos  --  -T  V, 

3 

cos  - .7  <! , ; 

i 

ill* 

o; 

0,00 

0,06 

| 0,00 

1 

1 

i 

1 

10 

0,02 

-0,10 

: 0,14 

1 

0,93 

0,95 

0,96 

20 

0,06 

-0,33 

: 0,26 

0,99 

0,73 

0,81 

0,84 

30 

0,13 

-0.50 

i 0,37 

0,95 

0,38 

0,64 

0,66 

•10 

0,23 

-0.67 

0,44 

0,86 

0.00 

0,48 

0,45 

50 

0,37 

-0,84 

0,48 

0,67 

-0,41 

0.41 

0,22 

60 

0,50 

-1,00 

0,50 

0,39 

-0,69 

0,37 

0,02 

70 

0,66 

— 1,14 

0,48 

0,00 

-0,90 

0,43 

-0,14 

80 

0,83 

-1,27 

0,44 

-0,37 

-0,99 

0,48 

-0.29 

90 

1,00 

-1,37 

0,36 

-0,71 

-1,00 

0,66 

-0,35 

100 

1,17 

-1,44 

0,26 

-0,93 

-0,97 

0,81 

-0,36 

110 

1,34 

-1,48 

0,14 

: -1,00 

-0,95 

0,95 

- 0,33 

120 

1,5 

-1,50 

0,00 

! -0,92 

-0,924 

i 1 

-0,28 

130 

1,64 

-1,48 

-0,10 

i -0,76 

-0,95 

1 0,93 

-0,25 

140 

1,77 

-1,44 

i -0,33 

: -0,59 

-0,97 

i 0,73 

-0,24 

160 

1,87 

-1,37 

, -0,50 

-0,31 

-1 

i 0,36 

-0,32 

160 

1,94 

-1,27 

-0,67 

-0,12 

-0,99 

i o 

-0,37 

170 

1,98 

— 1,14 

\ -0,84 

| -0,06 

-0,90 

! -0,37 

-0,44 

180 

2,00 

-1,00 

| -1,00 

i 0 

-0,69 

: -0,69 

-0,46 

In  agreement 
correction  tern 
nont 


\dth  the  general  procedure,  it  is  seen  that  the 
2 JQ  has  no  influence  to  <*  = 60°  . If  the  require- 


n sh 


2nd 

A 


+ 2 


is  not  fulfilled,  the  secondary  maximum  can  assume  considerably  greater 
values  than  the  extreme  values  given  by  J0(x) , 


Fig.  29.  Directional  characteristic 
of  the  compensated  circular  group 
with  a sufficient  number  of 
radiators. 


Fig.  30.  Directional  characteristic 
of  tho  compensated  circular  group 
vith  an  insufficient  number  of 
radiators. 


As  an  example,  we  calculate  the  case  n = 6,  djX  = 1,5  by  the 
approximation  fomula  (Fig.  30).  For  comparison  the  directional  char- 
acteristic for  n = 14,  d/x  — 1,5  which  is  given  simply  by  ./0(3jtsin«/2) 
is  also  plotted  (Fig.  29).  Both  directional  characteristics  agree 
completely  as  to  their  principal  parts  (i.e.,  to  the  first  minimum) 
but  then  depart  considerably  in  the  magnitudes  of  their  secondary 
maxima. 


B.  The  radiation  factor 
(a)  At  a fixed  frequency 


The  directional  characteristic  has  its  practical  significance 
when  it  is  a question  as  to  the  accuracy  with  which  a beamed  re- 
ceiver system  (which,  e.g.,  can  be  rotated)  can  locate  a distant 
sound  source.  Also  the  behavior  of  the  directional  characteristic 
vdll  determine  the  freedom  from  disturbance  in  certain  directions. 
From  the  transmission  viewpoint,  however,  beyond  the  question  as 
to  the  total  power  radiated,  it  is  generally  a question  of  con- 
centrating the  sound  transmission  in  a certain  direction  or  plane 
in  order  to  increase  the  efficiency  of  the  array.  Here  the  ques- 
tion is  how  great  is  the  sound  concentration  in  the  given  direc- 
tion as  compared  with  that  for  nondirectional  sound  radiation. 

In  order  to  calculate  the  total  power  of  a beamed  radiator, 
we  proceed  in  a manner  similar  to  that  used  in  introducing  the 
directional  characteristic  in  which  we  start  with  fomula  (3) 
for  the  nondirectional  system  and  then  calculate  the  influence 
of  the  beaming  by  the  use  of  a factor  6 which  we  designate  as 
the  radiation  factor . Thus  for  the  total  radiated  power  L 
we  obtain  the  relation 


1 F2 

L = L03  = -2-cattu&y  S.  (57) 

Here  <3  is  defined  by  the  integral  over  a sufficiently  large 
sphere  K with  the  radius  R: 


© = < (58) 

A' 

where  $ is  the  directional  characteristic  defined  by  (6)  and  dK  is 
the  surface  element  of  the  sphere  with  the  radius  R . Now, 
with  directed  sound  radiation,  the  sound  energy  passing  through 
the  unit  surface  of  the  sphere  with  the  radius  R in  the  direc- 
tion defined  by  the  directional  characteristic  D?  is  given  by 


2c  a 


/>!«* 
2 co 


1 

= 2 o • an  • 


«•*./’»  fit* 

A*  'in#' 
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(59) 


On  the  other  tiand  if  the  total  sound  power  given  by  (57) 
radiated  spherically,  the  sound  energy  passing  through  a unit 
face  would  then  be  given  by 

L 1 wlF1  <3 

471-8“  2C'a'n‘  A*  ‘ inW 


is 

sur- 


(60) 


For  each  direction  characterized  by  the  directional  character- 
istic the  ratio  of  from  (59)  to  L/inR*  from  (60): 


1, 


J>‘ 

2ca  ‘ 47I-R2 


91“ 

<S 


(61) 

states  how  many  times  as  great  the  sound  energy  is  in  the  considered 
direction  as  compared  with  that  for  spherical  propagation.  Gener- 
ally, it  is  referred  to  the  principal  direction  for  which  yt  = 1 , 

When  this  is  so  we  will  denote  the  thus 
standard  quantity  ]/£  as  the  condensation 
factor  f . 


4- 


U 


4 


ator 


We  will  illustrate  the  significance  of  the 
radiation  factor  in  a simple  example.  We 
consider  two  equal  radiators,  small  compared 
to  the  wave  length,  and  assume  that  the  total 
rotated  power  of  the  system  consisting  of 
the  two  radiators  is  to  remain  constant  while 
we  change  the  distance  between  the  two  radi- 
ators. On  the  basis  of  symmetry,  the  power 
radiated  by  each  of  the  two  radiators  then 
naturally  remains  constant.  «e  will  then 
find  the  sound  pressure  amplitude  at  a great 
distance  R on  the  normal  axis  of  the  system 
3l).  One  could  then  be  led  to  the  following  fallacy:  Radi- 
(1)  produces  the  sound  pressure  amplitude: 


fig.  31. For  illus- 
trating the  radi- 
ation factor. 


(Fig. 
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(61a) 


Likewise  radiator  (2)  yields  tho  pressure  amplitude 


v; 


(S  __  1 


2 AW 


' • F . 


(61b) 


Since  tho  signals  from  the  two  radiators  are  in  phase  on  the 
normal  axis,  the  resultant  sound  pressure  amplitude  at  P must  be 
equal  to  twice  that  due  to  an  individual  radiator  (independently  of 
the  mutual  distance  of  the  radiators).  But  we  know  that  a direc- 
tional effect  dependent  on  tho  distance  of  the  two  radiators  - 
i.e.,  a sound  condensation  variable  with  the  distance  of  the  radi- 
ators - exists  on  the  normal  axis.  The  error  is  due  to  the  fact 
that, with  constant  sound  radiation  of  the  individual  radiator,  its 
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velocity  amplitude  is  by  no  means  independent  of  the  sound  pressure 
produced  by  the  second  radiator  on  its  membrane.  In  addition  to  the 

•vork  which  the  individual  radiator  must 
do  with  an  undistorted  sound  field,  it 
must,  in  the  presence  of  a second  radi- 
ator, overcome  the  pressure  exerted  by 
this  on  its  membrane.  This  additional 
‘work  will  be  different  with  the  phase 
difference,  i.e.,  with  the  distance 
between  the  two  radiators,  (^uite  sim- 
ilar to  this  is  the  case  where,  in  a 
half  space  cut  off  by  a rigid  wall,  we 
bring  a radiator  of  constant  sound  power 
nearer  and  nearer  to  the  wall.)  In  order 
to  investigate  this  quantitatively,  we 

Fig.  32.  For  calculating  the  **£  radnated  pot, or  of 

' radiation  factor.  the  system  consiati^  of  equal  radi- 

auors  (1)  and  (2)  at  the  mstance  d 
(Fig.  32). 


By  (57) 


L = n c - a-  n\ol 


2 .1  .T 

O = fiftt  fwdK  = iJs,/d?/Facos*(ni,'<*aa)mn«,  * 


®=i(' 


sin  2 nd/X\ 


Big*  33,  curve  1,  shows  the  dependence  of  the  radiation  factor 
on  •./a  . If  v.e  were  able  to  keep  «•„  for  the  two  raciiators  constant, 
according  to  (62)  the  radiated  power  ‘would  vary  in  the  same  manner 
as  3 in  Fig.  33.  Conversely,  if  we  keep  L constant  and  vary  the 

distance,  >r„,  mu3t  vary  with  1.3 
since,  with  constant  L , 
must  remain  constant.  And,  with 
constant  total  radiation,  the  sound 
intensity  retains  its  maximum  in 
the  symmetry  plane  if  the  condensa- 
tion factor  / -1/3  is  a maximum. 

We  easily  compute  that  this  is  the 
case  ford//.  =0.715  ana  that  the 
maximum  value  of  f here  becomes 

Fig.  33.  Radiation  factor  of  a sys-  ?pua^“  2*554*  This  moans  that 
^ WQ  radiate  a definite  sound 


tern  consisting  of  two  radiators 
(distance  d) 


power  by  two  individual  radiators  - 
each  nondirectional  in  itself  - 


1.  In  phase.  2.  In  phase  opposition.  nondirectional  in  itself  - 

/ 2 instead  of  using  one  nondirectional 

The  factor  R should  also  appear  in  the  numerator  of  the 
fraction  t-reoedini*  the  double  integral  in  the  right  member 
of  this  equation. 


radiator,  then  in  the  maximum  case  2.55-fold  sound  energy  per  unit 
surface  can  be  attained  in  the  symmetry  plane  as  compared  with  the 
nondirectional  sound  radiation.  In  the  same  manner,  we  can  inquire 
as  to  the  maximum  condensation  factor  with  arrays  of  3,  4,  5,  etc. 
vibrators  arranged  on  a straight  line  equidistantly  and  always  find 
one  quite  definite  value  of  d/X  for  which  this  is  the  case.  This 
yields  the  table  given  below. 


n 

/max  = j/Srnat 

ffiT  d/X 

2 

= 2,55 

= 0,715 

3 

= 4,25 

= 0,77 

4 

= 5,9 

= 0,825 

;> 

= 7,7 

= 0,865 

0 

= 9,5 

= 0,90 

Fig. 

For  n = 
34. 

2,  3,  4 

If  we  denote  the  radiation  factor  for  n 
similarly  arranged  radiators  by  3, „ there 
thus  results  the  general  formula''^’ 


n-  i 


m “ 1 


. sinm2wrf/A\ 

mT2 Sd/AJ- 


(64) 


is  represented  in 


Fig.  34.  Condensation  factor 
of  the  straight  line  group 
(d  = distance  between  two 
radiators). 

1.  Two  radiators,  2.  Three 
radiators.  3.  Four  radia- 
tors. 4.  Five  radiators. 

5.  Six  radiators. 


If  d is  equal  to  X/2  or  an  inte- 
gral multiple  thereof,  then  from  (64): 


S„  = 1 jn, 


(64a) 


i.e.,  the  condensation  factor  is  then 
equal  to  n . This  has  been  pointed 
out  by  lord  ^ayleigh  The  calcula- 
tion of  the  circular  piston  membrane  is 
also  due  to  hijir4?i  For  this  it  turns 
out  that 


G — 
— 


(TldlW 


2nd jX  \} 


(65) 


wherein  d is  the  diameter  and  J t 
denotes  the  Bessel  function  of  the  first 
order.  From  this  it  results  that  for 
large  d/  X > the  condensation  factor  f 
of  the  circular  piston  membrane  is 

I — 2n  FI?.*,  (66) 

where  F denotes  the  surface  of  the 
membrane. 


^ Ann.  d.  Phys. , Vol.  7 (1930),  p.  964. 

^ On  the  nroduotion  and  distribution  of  sound.  Phil.  Mag,,  1903, 
pp.  289-305. 

^ ’The  theory  of  sound,  sec.  302. 
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Also,  the  radiation  factors  of  the  circular  membrane  whose  ve- 
locity amplitude  is 

= (66 

may  be  calculated  with  the  aid  of  the  Vessel's  functions.  For 
u?i  — 1 — s2lr 2 there  results  the  radiation  factor 


(5  = ‘>2  o'|  . 91  (* * i _ 5 (2ar)  2i/3(2x)\ 

* --  ^•(2,  + I!3!  ■ j; 


and  for  w2  = (1  — n2/r2)2  .'the  radiation  factor  is 


= 23-3!  3'f— ~ + + 4</6(2x)\  (68 

\ 3!  ~ 1 ! 4!  ' 2!  3!  a;*"  i"7  / ' 


In  general,  by  the  series  development  for 


it  turns  out  that 


« _ , _ (*/2)*  ^ (x,2)* W2)6  , 

~ n + 2 ' 21  (»  + 2)(»  + 3)  3 ! (n  + 2)  (h  + 3)  (n  + 4)  ' ’ 

^ ,__x*  , (2n+_5)**  _ (6< 

'is"  ~ 1 3 (#  + 2 j ' 2 ! 3 • 5 ■ (n"+  2)»  (»  + 3') 

(2n  + 7)x«  __  _ 

3!  3 • ft  • 7(n  4-  2>l  • (w  + 3)'(»  + 4)  + 

The  radiation  factors  S,,  S„  S,  are  represented  in  Fig,  (35). 
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Fig.  35.  1.  Radiation  factor  of  the 
circular  piston  membrane  («-». 

2,  Radiation  factor  of  the  stretched 
membrane  (»» l-e’w. 

3.  Radiation  factor  of  the  stretched 
membrane  i«  - u - sWi. 


Fig.  36.  For  the  calcu- 
lation of  the  radiation 
factor  of  the  rectangular 
piston  membrane. 


For  a rectangular  piston  membrane  whose  one  side  (b)  is  small 
compered  to  tne  wave  length  while  the  other  side  (c)  is  arbitrary 
(Fig.  36),  the  radiation  factor 


(cos  y) 


& ■ 


sinyiy  : 


© 


sin*cJi/A 


2 c 7/1 


(cn/A)2 


, 1 fain  t Jl amsc?r/A  1 c,2  c.t 

cn/XJ  t (i c«/A j*  ‘ cnjX  1 A 


(70) 


results . 


If  one  replaces 

the  sine  integral 

% 

~ , faint  . 

Six  = / — dt 

(70a) 

by  the  approximation 

0 

Si*-!-  C-“(l-2/x»)-^, 

(70b) 

one  obtains 

© 

If  . , sm2a;  . cos2a:l  * cn 

= 21  » ~ */*  “ "27-  + TF-  > x~-x- 

(71) 

For  larger  values  of  z(xs2)  , 
raat.ion 

© = £ or 


one  will  oe  able  to  use  the  approxi- 


l 

2x* 


as  is  seen  in  the  following  table 


X 

,S 2 1 »ln*  x 
x *• 

X 

A*<2/  hUi1/ 

X X* 

.7  l 

2 ~x  ~ 2 !rJ 

0 

1 

2,5 

0,563 

0,548 

0,2 

0,1)04 

3 

0,473 

0,467 

0,4 

0,983 

3,5 

0,406 

0,408 

0,6 

0,949 

4 

0,358 

0,361 

0,8 

0,931 

4.5 

0,323 

0,324 

1 

0,897 

5 

0,295 

0,294 

1,2 

0,858 

5,5 

0,27 1 

0,269 

1,4 

0,813 

0 

0.249 

0,248 

1,8 

0,76(1 

0,5 

0,230 

0,230 

1.8 

0,710 

7 

0,214 

0,214 

o 

0,073 

7,5 

0,200 

0,200 

The  third  term  in  the  brockets  shou  ld  be  - ~ . The 
last  term  in  the  brackets  should  be  * 
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The  maximum  condensation  factor  f then  turns  out  to  be  2 cjX 
for  large  values  of  cnj).  . This  therefore  means  that  a narrow 
rectangular  piston  membrane  (in  a rigid,  wall)  with  one  side  equal 
to  10  wave  lengths  sends  out  in  the  symmetry  plane  20  times  as 
much  sound  energy  per  unit  surface  as  compared  with  nondirectional 
radiation. 


If  b is  not  small  compared  to  the  wave  length  © can  generally 
be  represented  by  the  following  series'^. 


~ / cti \ 1 23  / 6 n \ 2 icn\ 


bny 

T! 


9 2 


)- 


(72) 


where 


Fig.  3?.  Function#  for  the  calculation  of  the  radiation 
factor  of  the  rectangle. 


Those  functions  aru  represented  in  Fig.  37*  Here  it  is  seen 
that  U3  soon  as  cj).  bocomos  greater  than  5»  w v’8.  •••  can  be  re- 
placed by  v’o  • •'•t  then  follows  from  (72)  that 


A.rm.  j.  Phys. , Vol.  7 (lyoO),  pp.  ‘J53-y57. 
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©=  0>o 


(74) 


The  integral 


is  found  by  the  series  development: 

. , ft) +$±  i£L+ 

v 3 * 1 ! -2 ! 't'  5 • 2!  3!  7-314! 


(74a) 


(75) 


or, for  the  larger  values  of  bnj).  , better  by  the  Bessel  function 
series 


* = in/f 


(76) 


For  the  larger  values  of  6 nil  this  series  is  considerably  more 
convenient  to  evaluate.  It  turns  out  that  for  6n/A>3  the  value 
of  the  bracketed  sum  differs  from  unity  only  by  (a  few)  percent^ 

Mio  value  of  the  function 


x • y (z) 


is  represented  in  *'ig.  38. 


j»/a 

2 /' sin* (ae cosy)  , 
nj  *eos8y  ™ 
o 


(77) 


Fig,  ?8.  Auxiliary  function  for  the 
calculation  of  the  radiation  factor. 


This  i«  aajudged  to  ' e the  Penning  of  the  following  sentence 
in  vhich  one  or  rrore  vorae  have  rp.jprently  been  left  out:  *Unu  zwar 
i r^ibt  rich,  dene  filr  bir/%  > 3 oer  Wert  urr  in  der  eckigen  Kleirmer 
atisiieriden  Suffire  nur  urn  Prosente  von  Bine  sbweicht." 
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If  we  therefore  assume  that  c/A>5  and  bjl  > l 
densation  factor  then  becomes 

, _ cn/X  • bjt/X  2 nF 

1 ~ n(2  ~~V~- 


the  con- 


(78) 


The  oondensation  factor  is  represented  for  values  of  bnjX 
and  cji/A  between  0 and  10  in  Fig,  39.  Here  are  drawn  the  curves 


Fig,  39.  Oondensation  factor  of  tue  rectangular 
piston  membrane.  (The  numbers  on  the  curves  are 
the  condensation  factor.) 


for  which  the  condensation  factor  possesses  a constant  value.  We 
find,  e.g.,  the  same  condensation  factor  f = 7 Tor  a square  for 
which  bnH  = c ;r//l=s3!&as  for  the  rectangle  6^  = 0, 5:  cjt/A  — 10  even 
though  the  surface  of  the  rectangle  is  only  lialf  as  large.  For 
a < 1 and  6 jt/A  < 1 (quarter)  circles  result.  For  I > 10 * 
(in  the  a»re  central  part)  the  curves  are  hyperbolas.  This  last 

means  that  with  sufficient  accuracy,  one  can  write  f — 2nFlkl  . 

In  the  calculation  of  directional  charaoteristics,  we  have 
remarked  that  whenever  individual  parts  of  the  membrane  vibrate 
in  opposite  phase,  the  direction  effect  is  substantially  affected. 
This  influence  must  also  be  effective  for  radiation  factors.  As 
a very  simple  example,  we  calculate  the  radiation  factor  of  two 
radiators,  small  compared  to  the  wave  length,  at  the  distance  d 
which  vibrate  in  opposite  phase.  A computation  analogous  to  that 
on  page  43  yields 
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1 /i  _ sin2«rfM\ 
® 2 \ 2nd/X  I' 


(79) 


We  obtain  the  corresponding  curve  by  reflection  of  the  curve 

(63) 

on  the  horizontal  line  y = ^ (see  Fig.  33)* 

The  maximum  of  f occurs  here  for  -^  = 7,725,  i.e.,  for 
d/X  = 1,23  and  it  amounts  to  2.294.  This  is  the  condensation  factor 
at  the  position  where  9J  = l and  corresponds  to  y = 24°  . While, 
with  two  in-phase  radiators,  we  could  obtain  a 2.55-fold  condensa- 
tion (for  y = 0°  ),  vdth  two  oppositely  phased  radiators  there  re- 
sulted at  most  a 2.29-fold  condensation  (for  y = 24°  ).  For  the  ex- 
ample on  page  30  where  wm  vanished,  the  radiation  factor  must 
naturally  lose  its  meaning  and  we  must  calculate  the  total  power  Lx 
by  the  formula  (59) • Then  from 

it  follows  that 


L^2¥n*,Fi -SJi»(z)]2dA'  * 

K 

. F*  1 J,{2*)  . 4 J}(2*)  8Jt(2i) 

= inca -jr  rp-  ^ + “?  "**  . 


(80) 


(81) 


If  the  Bessel's  functions  are  replaced  by  the  corresponding 
power  series,  there  results 


F* 

£,  " 4,tcct  ^ 


1^.  2 • i.-4 

3!  0!  " 


2-3**  3 • 4 ■ J'B 

4!  7‘f  + "'5!  81 


(82) 


If  the  power  radiated  by  a piston  membrane  of  the  same  size 
with  the  velocity  amplitude  w = 1 


2 n-c-a  • 


F* 

X*' 


M**)y 


is  compared  with  1,  for  small  values  of  x , one  finds 

t,  __ 

L,  ~ 3!  0!' 


(83) 

(83a) 


* X8  should  replaoe  X^  in  the  denominator  of  the  fraotion 

preceding  the  integral  sign. 
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this  becomes 


For 


r = 


nd 

T~ 


LJL2  = 1/17280. 


(83b) 


This  means  that 
power  as  that  for 
must  be 


if  the  membrane  for  is  to  radiate  the  same 
L2  , then  the  amplitude  of  the  first  membrane 

V 17  280  = 131 


times  that  of  the  latter. 


For  very  large  values  of 
(81)  and  (83): 


d/l  , on  the  other  hand,  we  have  by 


(83c) 


(b)  With  artificial  compensation 


From  Fig.  27,  it  is  to  be  seen  that  the  radiation  distribution 
is  also  substantially  changed  by  artificial  compensation.  In  general, 
the  radiation  factor  for  the  group  consisting  of  two  radiators  is 

n H'^/2 
dip  COS2 

0 -fi/2 

If  we  introduce  a new  variable  of  integration  x by 


7zd 


(siny  — 8iny±)  oos ydy. 


(*) 


it  is  found  that 


a:  — -r-  (ainy  — Biny*) 


(84a) 


nd 

~ Qo^xdx, 

n d 

r T(1  - W»r.) 


2nd 


(85) 
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sin 


In  Fig.  40,  <5j.  is  represented  for  sin  y*  = o,  4,  4,  f , 1 . Since 

yt  = 0 corresponds  to  the  uncompensated  case,  rt  is  seen  that  for 


Fig-  40.  B&diation  factor  for  two  compensated 
radiators  {Distance  d»  compensation  angle  yj). 

1.  Rln  Yk  ~ 0;  2.  sin  ;•*  --=»  J ; 3.  sin  y<  -•  1 ; 4.  sin  vk  --  2 ; C>.  sin  y*  *=--  I 


dfi  < i the  uncompensated  case  always  yields  larger  values,  £.g., 
for  dp.  = |j(s?  = 0,83;  <Bk  = o,5. 

From  this  it  follows  that  the  condensation  factor  of  the  com- 
pensated group  in  the  direction  of  the  line  of  the  radiators  (i.e., 

Yk  = 90° ) is  1.66  times  as  great  as  the  condensation  factor  of  the 
uncompensated  group  in  the  maximum  direction.  If,  therefore,  with 
two  nondirectional  radiators  at  the  distance  d < A/2  , the  greatest 
sound  intensity  for  one  definite  direction  is  desired,  the  best 
results  will  be  obtained  by  working  with  a compensated  array. 

In  general,  for  a compensated  straight  line  group  which  con- 
sists of  n equidistant  radiators,  the  radiation  factor  turns 
out  to  be 

0 1 I “ VV,  (85a) 

n h* “ tit  ■ 4;r<//A 

m 1 


If  the  direction  of  the  compensation  line  coincides  with  the 
direction  of  the  array  (i.e.,  >■,.  -•.=  wr  ),  one  finds 


l>  v 


/i  --  m 


niit  (in  • 4 ml'). 
hi  • 4 7i  i//A 


(86) 
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Comparing  this  with  the  radiation  factor  for  the  uncompensated 
case  (i.e.,  >^.  = 0) 


n-  i 


n*  2 {n  ~ m) 


8mm 


w = i 


r-r-j 

' indj). 


(87) 


we  conclude  that  the  straight  line  uncompensated  group  with  the 
receiver  distance  d has  the  same  radiation  factor  as  the  group 
compensated  in  the  array  direction  with  the  receiver  distance  d/2. 
Or,  otherwise  expressed,  the  sound  concentration  of  this  straight 
line  uncompensated  group  on  the  normal  axis  is  equal  to  the  sound 
concentration  of  the  group  compensated  in  the  array  direction  with 
half  the  receiver  distance. 

In  order  to  find  the  radiation  factor  for  the  densely  covered 
circular  group  with  a compensation  direction  we  have  to  cal- 

culate the  integral 

2/i  .7,2 

©*  = I d<p  I dyW*  sin y 
6 6 

where,  according  to  (52),  JR  is  given  by 

$ = JQ  )'  (cob«  — cosa0)2  + (cos/S  — cos/?0)2j 
oince  cos  « = sin  V COS  <p , cos  (I  = sin  y sin  <p  it  follows  that 
JR  = J0{kr ) sin2)/  — 2sinysiny0cos(9?  — <p0)  + sin8y0). 


(88) 


(88a) 


(8bb) 


If,  for  abbreviation,  we  set  « -■  k rainy,  v - lrsiuy0,  it  then  follows 
by  the  addition  theorem  for  Bessel  functions  that 


JR  = J0  (}'  M2  -j-  l’2  — 2 K V C08  {(f  — (f ;„)) 

'XJ 

- J0WJ0W  + 2 2l,V„(!i)JN(i')cos7i(^  - f0) 

N«1 


(89) 


If  we  now  form  the  integral 

2/i 

u 


2.i 

f cos  m ((/■  — vu)  cos  71  (9-  — If0)  (If 
0 


| 0 fiir  771  4=  w 
I fiir  m - 77 


(89a) 


(89b) 


-54- 


then  since 


the  terms  with  unequal  indices  fall  out  and  there  results 


2n 


^-jWdcp  = Jl(u)Jl{v)  4-  2J\(u)J\(v)  + • • 


It  than  turns  out  that 

?!  rz 

©*  = Jl(kramy0)  • f Jl(krainy)  sinyiy 


■•>/ 2 


+ 2 J](kramy0)  f J\(kr  siny)  8inydy+  ••• 
o 


or  since 


tt/2 


f JJ,  (rainy)  ainydy  = ^jj2n(Z)d£ 

0 0 

it  follovre  that 

. kr 

®t=JUkrainy0)  • ■Jj0(2S)d$  + 2J*(krsmy0)-±fja(2£)d$ 
o o 

k r 

+2J|(Af8iny0)-j-Jjt(2f)if+ 


The  calculations  of  the  integrals 

X 

y„(x)  = f Jin(2x)dx 


may  be  easily  carried  out  by  the  relation 

I J,('2r)dx  = 2^J,r+s„+1(2x) 


with  the  aid  of  the  aessel  function  tables. 


(90; 


(91] 


(92) 


(93) 


(93a 


(94) 


(92)  is  obtained  if,  in  the  well-known  equation 

"i  a 

9l(u:Biny)  = ® j J-t.{2xainycoaq>)(l<p 
u 


(a 


both  sides  aro  multiplied  by  »inyrf><  and  are  integrated  from  0 to  */2 
and  the  relation  piven  by  N'lelsen  (Kandb.  der  Zylincerfunktionen, 
p.  280,  formula  l,l»>  =--(>]) 


*.3  .t  2 


I j x sin  10/ {x  siivm  siny)  dw  dy  ■=  * /(*) 


(b) 


u o 


Is  used. 
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The  corresponding  functions 


■Mi 


are  represented  in  Fig.  41. 


In  Fig.  42  the  radiation  factor  (for  the  circular  group) 


is  drawn  as  a function  of 
wave  length).  Here: 


2 nr 


( r = radius  of  the  circle,  X = 


Curve  1 represents  the  uncompensated  case  (i.e.,  yk  = y0  ==  o°), 
Curve  2 represents  the  case  yk  — 30  , 

Curve  3 represents  the  case  yk » 90°. 


Besides  this,  the  curves  4,  5>  6: 

®.(-D 


resulting  from  this  are  drawn. 


2nr  , 
"A 


From  the  latter,  it  is  recognized  that  with  increasing 

2nr  2nr„  1 2nr\ 

A > A H A ) 

approaches  nearer  and  nearer  to  the  value  1/2. 


This  means  that,  for  larger  values  of  r/X  and  with  great  di- 
rectivity, the  condensation  factor  / — l/<3 1 is  given  by  . 2 . 
'Otherwise  expressed  this  says;  The  length  of  the  circumference  of 
the  circle  measured  in  X , when  multiplied  by  2,  gives  the  size 
of  the  condensation  factor.  In  connection  with  the  earlier  con- 
siderations on  the  condensation  factor  for  uncompensated  radiator 
arrays,  the  straight  line,  a circular  line,  a circular  surface 
or  a rectangular  surface,  we  can  formulate  the  following  general 
theorem; 


With  great  directivity  the  condensation  factor  is; 


1 . Two  times  the  length  of  the  radiating  line  measured 
in  wave  lengths  for  line-shaped  radiator  arrays . 


2.  2 tt  times  the  surface  of  the  radiating  arrangement 

measured  in  for  surface  radiator  arrays . 
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Part  Two 


THE  SOUND  FIELD  IN  THE  NEIGHBORHOOD  OF  THE  RADIATOR 
3.  The  group  of  two  radiators 


For  the  previous  considerations,  the  assumption  was  made  that 
the  field  point  was  at  a sufficient  distance  from  the  radiating 
surface.  For  this  part  of  the  sound  field,  which  is  generally 
taken  into  consideration  first  in  the  usual  practical  problems, 
there  results  a simple  representation  when,  to  the  expression 
which  characterizes  the  nondirectional  radiation,  a factor  is 
affixed  which  is  independent  of  the  distance  of  the  field  point 
and  depends  only  on  the  direction  of  the  field  point  line. 
Furthermore,  this  characteristic  function  (the  directional  char- 
acteristic) was  only  dependent  upon  one  quantity  (e.g.,  the  func- 


The  sound  field  was  therefore  substantially  determined  by  this 
one  function  - the  determination  being  quite  general  for  ary- 
given  frequency  (wave  length)  and  for  any  given  dimension  of  the 
radiating  system  (r/A).  ^'he  conditions  for  the  calculation 
and  representation  of  the  nearby  field  are  considerably  more 
difficult.  First,  we  are  forced  to  calculate  the  adjacent  field 
at  so  great  a number  of  points  that  the  complete  field  can  be 
obtained  by  interpolation,  and  secondly  it  is  necessary  to  carry 
out  this  representation  for  each  particular  case  which  i3  char- 
acterized by  the  ratio  of  the  geometric  dimension  to  the  wave 
length.  The  diversity  of  the  problem  has  now  thus  become  sub- 
stantially greater.  will,  therefore,  have  to  confine  the 
representation  of  the  adjacent  field  to  special  cases.  The  sound 
field  will  be  represented  when  we  draw  curves  in  the  neighborhood 
of  the  radiating  system  which  correspond  to  a constant  pressure 
amplitude . 

It  must  first  be  made  clear  when  a field  point  is  to  be  re- 
garded as  belonging  to  the  nearby  field  and  when  this  is  not  the 
case.  The  term  "nearby  field"  could  lead  to  the  mistaken  idea 
that  it  is  only  a question  of  the  geometry  of  the  radiators  so 
that  one  could  perhaps  say  that  with  a radiating  circular  piston 
membrane  of  the  radius  , all  field  points  are  no  longer  to  be 
regarded  as  belonging  to  the  nearby  field  which,  e.g.,  are  at  a 
greater  distance  from  the  membrane  center  than  j.0  times  the 
radius.  Actually,  this  definition  is  not  sufficient  since, 
besides  this,  it  is  also  a question  of  the  wave  length.  In  order 
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to  investigate  this  more  closely  we  must  start  from  the  formula 
generally  valid  for  distant  and  nearby  fields  and  determine  under 
what  conditions  this  becomes  formula  (6)  given  for  the  distant 
field. 


We  assume  that  the  radiating  surface  of  the  membrane  is  in 
a rigid  infinite  wall  or  acts  as  a double  membrane,  where  the  one 
part  vibrates  symmetrically  to  the  other  part  at  each  instant  as 
is  represented  in  Figs.  1-e  and  1-f . Here  the  radiating  surface 
must  consist  of  sevsral  component  parts  situated  in  the  same  rigid 
wall,  or  of  several  individual  double  membranes  all  of  which  have 
a common  symmetry  plane.  We  assume  the  radiating  surface  to  be 
at  the  zero  position  in  the  XX-plane  and  the  velocity  amplitude 
given  by 


rn  - w(x,  y)eUot  , (95) 

^ere  we  will,  in  general,  assume  w(x,y)  to  be  a real  function. 
Physically,  this  means  that  all  oscillating  membrane  elements  pass 
through  the  zero  position  simultaneously,  and  reach  their  extreme 
positions  simultaneously  so  that  (besides  nodal  lines)  only  nations 
of  the  membrane  elements  which  are  in  phase  or  in  phase  opposition 
are  possible.  Basically  no  difficulties  exist  in  prescribing  the 
velocity  amplitude  at  every  point  of  the  membrane  with  respect  to 
amplitude  and  phase  when 


w(x,  y)  = u(x,  y)  + iv(x,  y)  (9$a) 

is  regarded  as  a complex  function. 


'Fhen  for  an  arbitrary  field  point  P in  the  upper  half  space 
(because  of  the  rigid  wall,  we  can  confine  ourselves  to  the  half- 
space  2 S 0 ) the  behavior  of  the  sound  pressure  according  to 
Rayleigh^  is  given  by 


Fig.  43.  For  the 
definition  of  the 
nearby  field. 


p = «*<-»  + *'».-? J w(x,y) e y-dF. 


(96) 


Here  the  integration  is  to  extend  over  the  sur- 
face F radiating  into  the  half  space,  r is 
the  distance  of  the  element  of  integration  dF 
from  the  field  point  P and  R is  the  distance 
of  the  field  point  from  the  coordinate  origin 
(see  Fig.  43). 


If  the  field  point  P has  the  coordinates, 
•ro.  ?/».  ~n » and  the  mid-point  of  the  surface 
element  dF  has  the  coordinates  x,  y , the 
relation 


rl  R"~  4 xa  4-  >/  — 2u\fu  - 2 </,//„ 


(96a) 


The  theory  of  sound,  Sec.  278. 


I 


m 


I 


w ' 


« i 


then  follows  from 


./?“  ' - .?n  j-  f/u  " t - 2,. , 

r2  -•=  (.t:  — ,r0)2  + (//  — y/„)2  4 22 


(96b) 


Therefore 


r 

R 


i _ 2(*cos*j-t-  ^cos^)  , z2  + y 2 1'1 

1 " it  ' ir-  I ’ 


(97) 


xcos  \ -f  ycos  /?  ] xJ -f  .t/2  ] /xcos  a -r  y cos  jl\z 

' " ~R"  " • 2 R3  ~ ■>  \ j • 


In  this  development,  terms  of  higher  order  than  the  second  are 
neglected.  In  place  of  the  earlier  formula  (which  was  derived  on 
assuming  a sufficient  distance  of  the  field  point): 


p = ci("‘t  + x.i-iR)  j V'(x,  l/)eik<-*°«  *■'-  W'vMdF 
}■’ 

there  now  follows  from  (96)  and  (97): 

P r_  ^ . J2-  i'«)  j yj  , f,iKi  1-.W.,  ■ Jjnw/fl 


(98) 


I A’  f1  v*  ( j* <•« is  \ i/nistt)1 
•J  « « 


xcoha-  i/ros/f  1 

1 H ' ? .>  l x- y - ixnw\-v  yvmflf] 


dt\ 


(99) 


(99)  can  therefore  be  replaced  by  (98)  when 


/Jrj/*  • i/1  (jt»k  \ 1/ 1*<  is*  /)’)* 

^ I /<  /< 


(99a) 


i - 


XCOK\  ■ l/COH //  I J-  • I/2  I /XCOS  \ t/COS/b2 


/f 


• 2 «» 


1 /j'l-u 

2 1 


7 


can  be  replaced  by  1.  Since  this  is  to  hold  for  all  values  of  \ 
and  /}  , it  is  easy  to  see  that  because  of  the  denominator,  one 
must  have 


R3 


1 


and  because  of  the  numerator  it  must  be  that 


(100) 

(101) 


here  the  exponent  of  (viv‘i<*  + lf  cos  ) in  the  second  exponential 
of  the  integrand  should  be  changed  from  3 to  2, 
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Since  in  the  most  unfavorable  case 


]x2  + y7  = qm 


(101a) 


is  the  greatest  distance  of  a membrane  point  from  the  coordinate 
origin.,  (100)  asserts  that  the  field  point  distance  R must  be  large 
compared  to  the  greatest  linear  dimension  of  the  radiating  surface 
(in  the  XY-plane ) . From  (101)  it  furthermore  follows  that  one  must 
have 


n ■ Cm 
A - R 


<1 


(102) 


With  snail  values  of  X (more  precisely,  if  /.  ,i  H ) condition 
(102)  is  therefore  more  discriminating.  If  we  assume  a piston  membrane 
of  5 cm  radius  and  first  use  a wave  length  of  15  cm  and  then  a wave 
length  of  1 cm,  then,  according  to  (100),  the  sufficient  distance  R 
would  be  given  by  R v f>  cm  while  in  the  second  case  it  would,  accord- 
ing to  (102),  be  given  by  R ...  7.3  cm . 


For  two  radiator's  which  are  small  compared  to  the  wave  length, 
a simple  addition  replaces  the  integration  of  formula  (10).  »Ve  ob- 
tain the  pressure  of  the  resultant  field  in  the  following  form: 


P =• 


Hal  i -i  j) ' 


2 -A 


"'l  A’,  ’ I If.,  F.,  ' 

r,  - - r, 


Id  « 


(103) 


Here, 

a are  the  mean  velocity  amplitudes  of  t\  and  t\, 
t\,Fs  are  the  radiating  surfaces 

rl,ri  are  the  distances  of  the  radiators  from  the  field 
point. 


If,  furthermore,  we  introduce  the  abbreviations 


«’i  f 1 
2AS 


11  , 


tr.  y. 
:!/• 


r 

) 


A 


then  there  results  from  (103)  the  relation 

;»  hr 

c ■ a ! x u 


r: 

/ 


u > 


(103a) 

(104) 


for  the  pressure  amplitude  p. 


Since  we  leave  the  phase  out  of  consideration,  we  have  to  in- 
vestigate the  expression  on  the  right  only  as  regards  to  its  magni- 
tude. By  a simple  calculation,  we  obtain 


/» 
C (T 


I",  ' I'f^rrU  - „) 


("  1 1'-"*  »<'■  "I 


(105) 
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or 


p 

c • a 


2 ah 
*‘J 


cos  2n(x  — ij)  . 


(105a) 


For  certain  values  of  (x  - y),  cos  2 w (x  - y)  (and  therefore 
also  ) assumes  simple  values  which  are  shown  in  the  following 
table: 


x -v 

1 2 cos  2 n(x  ~ y ) 

i 

0 

+ 2 

■i  x -r  hly 

2*  -/r 

-r-  1 

l«*  **  + Pif  + <tb:xy 

..  i 
...  i 

0 

| n'l.x'2  + b-,y- 

. i 

- 1 

■Ja'.'x'1  4-  Ii*'!/-  — abjxy 

r \ 

-2 

a:.r  - b y 

: Y 

- 1 

\<i*;x1  + h*  ;if  — ablxy 

i*'  i‘ 

0 

1 <i' lx"  4 lr,y- 

. 

T 1 

| a*  x2  -h  Ir.y2  ab/xy 

1 

O 

n x ■ hly 

: (l  t >) 

m 1 

pi\>'  i-  lt2jir  4-  ab/x  y 

etc. 


The  values  o'*  p/(cc)  corresponding  to  the  points  lying  on 
the  hyperbolas  x y = constant  can  thus  be  simply  calculated  and, 
by  interpolation,  curves  of  constant  pressure  can  be  constructed 
for  which  x and  y are  the  distances  measured  in  wave  lengths  of 
the  field  point  from  the  two  radiators. 

Of  particular  importance  is  the  determination  of  the  positions 
where  the  pressure  amplitude  vanishes.  ( For  this  it  is  manifestly 
necos3ary  that  'f  " and  that  ,■  .... ,,  1 (m  u,  i,  2,  . . .)  . 

(«itixmt  limiting  the  generality  a can  be  assumed  greater  thanb.) 


From  this  it  follows  that 


X 

y 


</  2 m • l 

ii  h 2 


Ji  2 m • I 
a b ‘ 


(105b) 


iloreover,  in  order  that,  the  circles  described  by  x and  y yield 
a xHial  intersection  point,  the  condition 


.r  • »/  i / /.  r ii. 

mu3t  be  fulfilled  for  a given  distance  d of  the  two  radiators. 
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(105c) 


That  is,  the  condition 


« + i 2»  + I^i^2j«  + I 

a — fc " 2 2 (106) 


must  be  satisfied. 

To  each  value  of  m for  which  the  inequality  (106)  is  fulfilled 
corresponds  a zero  position  of  the  pressure.  If  dll  < J there  are 
therefore  no  zero  positions  at  all.  If  a = 2,  6 = 1,  and  dll  — 4, 
then  from 

3(2  -f  l)  Ss  8 s 2m  + l (106a) 

there  result  the  solutions  m = l.  m — 2,  ™ = 3 ♦ If  a = 3.  b = 2, 
and  d{l  = 1 then  from 

6 (2m  + 1)  fe  2 2m  + 1 (i06b) 

only  the  solution  m = 0 results . 


Tig.  44.  The  zero  positions  of  the  eound 

pressure  of  two  point  radiator*  at 

the  distance  d ■ 6Xfor  different 
deformation  volumes. 

Fig.  44  represents  the  conditions  for  d/X**5,  a«3and  b « 2, 
Here  the  hyperbola  branches 

r,  - ra  “ A/2,  3/2  A,  5/2  A,  7/2  A,  9/2  A (106c) 

and  the  zero  positions  (1),  (2),  (3),  (4)  lying  on  them  are  drawn. 
These  zero  positions  are  the  points  of  intersection. of  the  circles 
which  are  described  about  the  point  radiators  Ax  and  A%  with  the 
radii 

r,  3 (2  m n <i 

>:  To  x 

2 (2  m 4- 1 ) <1 

"""  id  ~ A 
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r, 

X 


(m  — 1,2,  3,  4) 


(107) 


Moreover,  all  the  zero  positions  lie  on  one  circle  which  divides 
the  distance  AiAz  harmonically  with  the  ratio  a:b  and  whose  center 
lies  on  AtA2  (see  Fig.  44).  The  radius  of  this  circle  is 

r _ ab  « (107a) 

— b2  T‘ 


A representation  of  the  curves  of  equal  pressure  may  be  attained 
by  graphical  means.  For  this  purpose  wa  draw  the  locus  of  the  termini 
of  the  vector 


(107b) 


and  of  the  vector 


(107c; 


This  was  done  by  allowing  x or  y to  increase  by  0.05  so  that  each 
vector  followed  from  the  preceding  by  a rotation  of  2 ^ *0.05  = 1£P. 
One  then  needs  only  to  draw  the  straight  lines  through  the  zero  point 
intersecting  each  other  at  an  angle  of  18?  and  to  lay  off  on  them  the 
lengths  a/x  or  b/y  as  the  case  may  be.  Two  spirals  thus  result  which 
wind  around  the  null-point  with  ever  closer  windings.  If  these  spirals 
have  been  numbered  with  the  corresponding  x and  y values,  all  the  so- 
lutions x,  y of  the  equation 


a 


(10  B) 


can  be  given  when  a straight  line  of  the  length  o moves  so  that  the 
initial  point  slides  on  the  one  (x)  spiral  and  the  end  point  slides 
on  the  other  (y)  spiral.  4ich  position  of  the  line  o defines  by 
its  initial  and  end  points  on  the  spirals  a system  of  values  x,  y 
which  satisfies  r;uation  (14)'^'* 

As  an  example,  we  choose  a~b--l  t i.e.,  two  radiators  of 
equal  intensity.  In  Fig.  45>  the  spiral  r, e'-'*/*-  is  represented 
for  all  values  from  x ®0.5  to  x = 5.  In  order  to  clarify  tho  draw- 
ing, only  the  spiral  points  corresponding  to  the  individual  values 
y =0,5,  y =0.6,  y =0.7,  etc.,  to  y =2.5  are  represented  for 
tho  second  spiral  u < . Those  are  obtained  very  simply 
by  reflecting  the  corresponding  points  of  tho  x spiral  vdth  respect 
to  the  conter.  'Iheso  latter  points  are  marked  by  small  round 
circles  near  which  tho  corresponding  figure  is  given  in  a square. 

If  we  describe  a circle  with  the  radius  unity  about  such  a point 
of  the  y-spiral  (o.g.,  the  point  l),  we  obtain  by  the  scale  given 
on  the  x-spiral  all  tho  values  of  x (lying  between  0.5  and  5)  which 
satisfy  the  equation 


^Apparently,  eq.  (104)  is  meant  here. 
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We  thus  read  from  Fig.  45  the  solutions  x = 0.50;  0.58;  1.32; 
1.70;  2.29;  2.71;  3*28;  3-72;  4.29;  4-72.  The  relation  becomes  still 
clearer  if  we  draw  the  curve  defined  by 

K*2""  ^ , ()  (108b) 


in  rectangular  coordinates . One  half  of  this  curve  is  given  in  Fig. 
46.  The  other  half  is  secured  by  reflecting  the  first  half  on  the 
straight  line  y - x = 0.  In  addition  to  this  the  curves  for 

iZ.y  J 3 A (109) 

y - 2 -°*  ! -----  + — p-j-  2 -0 


and 


+ 


y 


•2  =--  0 


are  drawn  in  Fig.  4 6'^'. 

In  order  to  obtain  from  these  curves,  which  are  independent  of 
the  radiator  distance,  the  corresponding  curves  of  constant  pressure 
for  a given  radiator  distance  (e.g.,  d/X  = 3),  we  have  to  draw  the 
two  points  Ax  and  at  the  distance  d/X  =3  and  describe  a.bout 
these  points  circles  whoso  radii  are  given  by  the  coordinates  x and 
y of  the  desired  point.  Here,  however,  only  the  coordinates  x and 
y are  to  be  considered  which  lead  to  real  intersection  points  of 
the  two  circles.  Manifestly,  this  depends  essentially  on  the  radi— 
itor  distance.  From  the  condition  for  the  intersection  of  the  two 
circles 

'■>  1 r*  ''  ' ~ '>  (lC9a) 

it  follows  that 

1/  i-  X (l,/,  1/  - . (109b) 


O Those  ou  rves  can  be  identified  as  follows:  For  their  inter- 
section with  y-x  = 0,  we  have,  by  the  piven  conditions  and  eq.  108: 

J- 1 2 - c 

X f 1 

or  . a 1 = ££ 

] civq  fcx  “ 1 ^ **'  2- 

The  value  of  the  express  * -n  on  the  left-hand  side  is  1.  hence  the 
abscissa  of  the  inter  .on  of  the  curve  with  y-x  - 0 is  x = £ 

Tii i s assumes  the  values  1,  4/3,  2 and  4 for  c — 2,  3/2,  1 and 
l/2  respectively. 
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(109c) 


If,  in  Fig,  46,  we  draw  the  three  straight  lines 

y + x — (!■/.,  y — x — rf/A,  y — x — 0 


(where,  because  of  symmetry,  we  assume  r.,  > r,  and  can  confine  our- 


selves to  one  quadrant),  a rectan- 
gular stri}'  is  then  bounded  by 
these  ; Ine:  which  contains  those 
and  only  tnose  points  x,  y for 
which  the  condition 


V -f  * (HI  y - x --I  0 (109d) 


is  fulfilled. 

The  corresponding  strip  for  d/X  = 3 
has  been  shaded  in  Fig.  46.  One 
recognizes  that,  e.g.,  of  the  curve 

■ it.tx  e i 2 1 v 

I , '•  , - 1 0 (109e) 

only  the  broken  part  enters  into 
consideration.  If  we  transfer  the 
four  curves,  insofar  as  they  are 
contained  in  the  shaded  strip,  we 
thus  obtain  the  corresponding 
curves  of  constant  sound  pressure 
(for  one  quadrant).  Here  each 
point  of  the  curves  in  Fig.  46 
with  the  coordinates  (x,  y)  yields 
a corresponding  point  in  Fig.  47 
as  the  intersection  point  of  the 
circle  described  about  -4,  with  the 
radius  x and  the  circle  described 
about  At  with  the  radius  y . It  is  important  that  one  recognizes 
from  Fig.  46  whether  the  constant  pressure  curve  consists  of  one  con- 
tinuous curve  or  how  many  separate  curve  sections  result.  Thus  for 
dA  0 3 the  constant  pressure  curvos  corresponding  to 


Fig.  46.  • The  functions 

i~!-~  ' ..r  Irir  c ■ j.  1. 2. 

/ V 


I _ -i  x 
x 


I . » u 
;/ 


(109f) 


yield  two  separate  curve  sections  for  c « 0.5  and  three  for  0 “ 1 
while  c “ 1.5  and  0 = 2 each  produce  one  continuous  curve  section. 

In  Fig.  47  the  corresponding  curves  of  constant  pressure  are  repre- 
sented. a he  complete  spatial  distribution  is  obtained  when  one  allows 
the  whole  structure  to  rotate  about  . Three  separate  surfaces 

then  result  for  C = 0.5,  five  for  c = 1 and  two  each  for  c = 1.5  and 
c ~2.  If  we  allow  d/X  to  assume  smaller  values,  the  corresponding 
siiaded  area  in  Fig.  46  becomes  steadily  narrower  and  the  number  of 


-11* 
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Fig.  47.  Curves  of  constant  pressure  amplitude 
for  two  point  radiators  at  the  distance  d = 3X, 


extreme  values  always  decreases  until  finally  for  d/X  c £ only  one 
maximum  (on  the  normal  axis)  remains  for  all  curves. 

In  Fig3.  48  and  49  the  sound  fields  are  drawn  for  d/X  * 1 and 
for  Fig.  48 

«’|F,  _ _ , 

2 A*  2 A2 

while  in  Fig.  49 


According  to  previous  considerations,  a null  point  (in  the 
spatial  sound  field,  a null  circle)  must  appear^  . One  sees  from 
Fig.  49  that  particularly  in  the  neighborhood  of  this  null  position, 
a fair  dissymmetry  of  the  sound  field  is  produced. 


v For  these  values  of  a and  & (i.e.,  a = 1,2,  ^-=0.9)  condition 
(106)  Is  satisfied  by  the  single  value  in=sO.  Consequently,  by  eq. 
(,105b)  Z=  1.  This  point  lies  within  the  0,05  curve  of 

Fig,  49.  In  the  previous  example  with  & = the  condition  for  a 
null  point  can  only  be  satisfied  at  infinity. 
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4 


The  circular  piston  membrane 


The  calculation  of  the  sound  field  of  a piston  membrane  for  an 
arbitrary  field  point  in  the  neighborhood  of  the  membrane  raises 

considerable  difficulties.  Therefore  we  first 
Investigate  the  case  where  the  field  point  lies 
on  the  normal  axis  of  the  membrane.  Then  the 
calculation  may  be  very  simply  carried  out'^jk 
Since  w(x,  y)  is  to  be  constant  ( = «%  ),  the 
calculation  of  the  integral 


- !-r 


— dF 


(109g) 


is  necessary. 


If  we  introduce  the  polar  coordinates  g,  <p 
for  the  surface  element,  then  dF  = gdgdf 

. and  since  r2  = £2  + «2  , gdo  = rdr.  (Fig.  50). 

Fig.  50.  For  the  calcu-  ’ ^ 

lation  of  the  circular  then  fon0W3  that 
piston  membrane. 

2t  « 


Jr= 


v/r  • .>• 

•tr  >'  2n 


2 n/er  ikrdr  — - ■.*  [«- + !•  - e- <**  . 


0 0 


^3ing  the  easily  to  be  derived  relation 


i ' (• " < /< 


. . \ /i  - ' 1 ' . , 
- — 2 1 sm  „ c 


there  results; 


| k | 1]  i *{l«'  e'  4 

./  2/.  sin  I j|  /f2  -f- 12  - - ;!!'• 


Upon  substitution  of  this  result  in  (11),  there  follows; 

fir,  I I Uf  ' - * (I  li<  l'  r hi 

v>  2<‘-a-  H-0-sm|  2 \\tf-  !•  1 - - 

The  relative  sound  amplitude  is  therefore, 

fg  - sin  { o ll  I"  'S  ~ ~lj  • 

From  this  it  follows  that  /<>-o  has  the  value  zero  for 

k | A’2  + :l  — kz  ■ 2.t  , 4.t,  ... 


(no) 

(110a) 

(110b) 

(111) 

(111a) 

(111b) 


Baokhauo.  H. , and  Trendelenburg,  F.;  Ober  die  Riohtwirkung  von 
Kolbenmembranen,  Z.  teohn.  Phys.,  Vol.  7,  p.  630  (1926). 
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and  the  value  2 for 


k]'W  + z2~kz  = 71,  3n,  ... 


(111c 


If  one  denotes  the  values  of  % corresponding  to  these  extreme 
values  by  ?0  and  z„,  the  null  positions  z0  are  given  by 


A2 


0 2 »/A 

and  the  "positions  of  the  maximum"  by 


(»=  l,  2, ...) 


*2  / i \* 

* (2»"+  1)/A  " • 


(112) 

(113) 


(»  — 0,  1,  2,  ...) 


With  an  Increasing  radius,  the  number  of  the  zero  and  maximum 
positions  increases.  No  extreme  values  can  occur  on  the  normal  axis 

for  membranes  whose  radius  is  smaller 
than  ^ . The  complete  behavior  of  p/c<7 
for  field  points  on  the  normal  axis  of 
the  piston  membrane  is  represented  in 
Fig.  51,  The  ordinate  is  the  ratio  z/R. 
The  four  curves  correspond  to  the  values 
kR  ■ 6,  10,  20,  and  40  respectively. 

For  points  with  a sufficiently  large 
value  of  z , formula  (111)  transforms 
into  (l)  according  to  a previous  con- 
clusion (p.4  )•  By  the  considerations 
on  page  60 , z is  sufficiently  large  if 


*/?2 

As" 


and 


Then,  however, 

tl'jp  + s*  kz  kz^ i -t  !>  * 

and 


~ k; 


n!{‘‘  (H3j 

As 


2 sin 


i ,T«2 


As 


st  A’2 
As 


F 

A •s’ 
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Fig.  61. 
tude  <w»«> 


Preeeure  ampli- 
-n  the  normal  axia 
of  the  circular  piston  of 
the  radius  B. 

3.  * ft  - 20, 


1,  kit  - 0,  2,  kit  10, 
4.  ktl  - 40. 


so  that  (111)  does  in  fact  transform  into 
(1)  (see  footnote  1).  If  the  field  point 
distance  z is  chosen  as  six  times  the 
radius,  there  then  results  for 

n A2  . 

As  (or  for  L ' ) in  the  cases 


JcR»6,  10,  20,  40,  the  values  £,  ^ so  that  a field  point 

distance  greater  than  three  membrane  diameters  (at  most)  can  be  re- 
garded as  sufficiently  great  in  the  case  kR  * 6.  This  does  not 
hold  by  any  means  for  the  case  kR  * 20  or  for  kR  =40. 


It  may  be  mentioned  that  formula  (111)  can  be  generalized  to 
the  case  where  the  radiating  surface  consists  of  a sector  of  a 
circular  ring  with  the  limiting  radii  and  n2  and  the  central 
angle  % instead  of  a complete  circular  surface.  If  the  field 
point  i3  then  located  above  the  center  of  the  circular  ring,  the 
formula  analogous  to  (111) 


•w. 


,•  c • cr  sin  |-|  [y*2 + Vz4+ i2?j}  p*  ‘ + : 


-i{y?r^+v««+fif}) 


(114) 


is  then  valid. 


If  the  pressure  at  one  special  point  for  a piston  membrane  with 
an  arbitrary  boundary  is  desired,  one  can,  after  resolving  the  radi- 
ating surface  into  the  corresponding  partial  domains,  apply  this 

formula  and  sum  up  the  effects  of  the 
individual  component  surfaces.  Here, 
one  has  only  to  see  that  the  neglected 
parts  of  the  surface  are  so  small  that, 
taken  together,  they  make  practically 
no  contribution.  The  subdivision  there- 
fore depends  essentially  on  the  magni- 
tude of  the  wave  length. 


Jig.  62.  for  the  calcu- 
lation of  the  sound 
pressure  on  the  edge 
of  the  circular 
piston  membrane. 

(Fig.  52)* 


Tha  calculation  of  the  integral  may 
be  simply  carried  out,  if  the  field 
point  lies  0£  the  ecigo  of  the  membrane. 
If  we  choose  the  field  point  £ at  the 
origin  of  the  coordinate  system  with 
the  polar  coordinates  ei  and  f then 


.■»/’>  2Jleob'*» 


ft^ldF  = 2 [d<p  f e~ik"'dQl  = fk  fd<p[  1 - 
;?  u o’  o 


(lUa) 


Now 


.t/2 

A le-it2Hcotlfd(p  = J0(2M)  - iU0(2kli).  (114b) 

0 
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Hero  J„  and  Ha  indicate  respectively  the  Bessel  and  Struve 
functions  of  zero  order.  It  then  follows  that^ 

P = C-  oe*'w'  [— f-~+  + i HQ(2k.R)\ . (115) 


With  the  aid  of  the  available  tables'^' for  J0  and  Hn  we  csja 
represertt  the  behavior  of  the  pressure  amplitude  on  the  edge  of  the 
membrane  as  a function  of  2nRj\ . In  Fig.  53,  in  addition  to  this 
representation,  the  behavior  of  the  pressure  amplitude  at  the  center 
is  indicated  which  is  given  by  (111)  with  z <=  0 by 


v 

ac 


— 2 sin 


nit 

T 


(116) 


The  two  curves  exhibit  a notably  different  behavior,  while,  at 

the  center,  the  relative  amplitude  p/ccr 
varies  periodically  between  the  values 
zero  and  two  with  increasing  membrane 
radius,  at  the  boundary  of  the  membrane 
p/cd  approaches  closer  and  closer  to 
the  value  £ with  increasing  radius. 

From  surface  considerations,  one  is  prone 
to  assume  that  with  a piston  membrane  a 
sound  field  is  developed  immediately  in 
front  of  the  membrane  which,  with  in- 
creasing membrane  radius  R (for  R » / ) 
corresponds  more  and  more  closely  to 
the  sound  field  of  a plane  wave.  Ac- 
cording to  the  foregoing,  this  is  by 
no  means  the  case.  As  is  well  known, 
there  is  a similar  fallacy  when  one 
allows  a plane  wave  to  fall  perpendic- 
fi«.  83.  Pressure  amplitude  ularly  on  a shield  with  a round  open- 
(p/oo)  at  the  center  (l)  and  and  believes  that  by  decreasing 
at  the  edge  (2)  of  the  the  size  of  the  opening,  one  can  di- 
circular  piston  membrane  minish  to  the  point  of  extinction  a 
of  the  radius  E.  steadily  contracting  acoustic  "ray" 
while  in  actuality  with  a decreasing  radius  a (e  ■?.'!)  of  the  opening 
a more  and  more  nearly  hemispherical  divergence  takes  place. 

We  have  already  seen  that  the  pressure  at  a great  distance  can 
be  simply  calculated  if  the  velocity  amplitude  of  the  circular  mem- 
brane w(q)  is  given  in  the  form: 

w(p)  - «0  f a,  ( i - i-  djj  -| t «»^i  — j . (Il6a) 


McLachlan,  W. t On  the  Aoousfcio  and  Inortia  Pressure  on  a 
Vibrating  Circular  Disk.  Phil.  Mag.  Ser.  7,  (1932)  p.  1022. 

^ Theory  of  Bessel  Functions,  Cambridge:  G.  N.  Watson  1922. 


Here  e is  the  distance  from  the  center  and  R is  the  membrane 
radius.  >*e  will  now  stow  that,  in  the  simple  case  n = 1 and  for 
field  points  on  the  normal  axis  and  on  the  edge  of  the  membrane, 
a simple  calculation  is  also  possible.  »»q  set 

«>(<?)  = 1 --/•£•  (117) 


Then  the  corresponding  pressure  amplitude  pf  is  given  by 


P/  = 


(117a) 


It  is  sufficient  to  calculate  the  pressure  amplitudes  p0  (for 
f = 0)  and  pL  (for  f »l).  Then  from 


Po 


(117b) 


it  follows  that 


P/  = I1  -/)Po  + /Pi- 


(117c) 


If  the  field  point  lies  on  the  normal  axis  (fig.  50)  at  the  dis- 
tance a than  gi  ~ r*  - z4,  and,  3ince  dF  =>  rdrdtp  , it  follows  that 


Vr'  + 1< 


Vr'  + 


(1174) 


t l 

We  therefore  find  for  the  field  point  on  the  normal  axis 

Pi  - | - SW"'kH?  i1  + (U8) 


and  by  (110) 


pQ  = e1*'  _ e- <*>'«• Vt«_ 


(118a) 


If  the  field  point  is  on  the  edge  of  the  membrane  then  sine® 
(Fig.  52) 


£a  = ifs  -f.  (>j  — 2p,iJcos(/', 

l 2 ci  c08o <?> 

/(3  H 1 P3  ’ 

(118b) 

•V«i 

pi"  - ->'/  «*?• 

, 2 e,  4*0 »q>  £>{ 

a w ■ 

(118c) 

» 0 
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If  we  integrate  with  respect  to  ^ 


> 


there  results 


V'C  = jg  • ~J  ji  -icos<p  - (o  + i cos<p)j  dcp  , 

U 

.t/2  »/2 


(118d) 


Using  the  relation^' 

.t/2 

.?-[e-'*™»dO  = J0(x)-iHoi  x), 

6 

nJ2 

~ j cosOe~  txCMUd9  — — //,(£)  — » Jj(x) 

6 


(U8e) 


it  then  follows  that  for  field  points  on  the  boundary 


1 - J0(2kR)  Jt[2kR)  , . U0(2kR) 

ft  = — FB* HOT  ‘+*pb,'-»,‘ 


J0(2kR) 

2 


tf,  (2  A- A) 
' 1-ii 


4 

71  ’ 


(119; 


Here  «/0  and  J!  are  Bessel  functions  of  the  zero  and  first 

order  and  H0  and  Hi  are  the  Struve  functions  of  the  zero  and  first 

order.  Witn  the  aid  of  the  available  tables  for  •/„,  J\>  H0,  Hx  the 
calculation  offers  no  difficulties . 

The  variation  of  the  pressure  on  the  normal  axis  for  the  case 
tR  * 10  is  represented  in  Fig.  55  and  Fig.  56  for 

tv  = J — / • (?a/A4,  (/  “ 0 » d:  ± i > ± 1) 

The  corresponding  velocity  amplitudes 

W=»l ~/-<>3/tfa  (1=0,  ± j . ±h  ±1) 

are  drawn  in  Fig.  54.  furthermore  the  variations  of  the  pressure  at 

the  center  and  on  the  boundary  are  represented  as  a function  of  kR 

for  w « l - {}»/*•  and  «■ » 1 - *#*/**  in  Figs . 57  and  5a. 


'\//  Molachl&a,  W.  j 


Bessel  funotiona  for  Engineers,  p.  167 


-75- 


Pig.  54.  Tho  velocity  amplitude* 

«•..  1 l/»0,  ±J,  i J,  ±1). 


Pig.  56.  Preeeure  amplitude  on,  toe 
normal  axle  for  the  velocity  amplitude 

>i  i it.  i . ~ 1 ii  : ;~i. 


Pig.  55,  Preeeure  amplitude  ■!><•  •>  on  the  normal 
axle  for  the  v*i""tty  amplitude 


These  tw  cases,  where  the  field  point  is  (a)  on  the  normal 
axis  and  (b)  on  the  boundary  of  the  circular  piston  membrane,  are 
the  only  ones  for  which  the  evaluation  of  the  integral  is  generally 
possible,  i.e.,  expressed  in  terms  of  well-known  functions  which  are 
available  in  tables. 


■nann 

HUHaUH 


tig.  57.  Pressure  amplitude  at 
the  center  (l)  and  at  the  bound- 
fry  point  (2)  for  the  vibration 
fora  v - 1 - ew  and  the  corre- 
sponding components  ta. . u»,,  j,). 


Pig,  58.  Pressure  amplitude  at 
the  center  (l)  and  at  the  bound 
ary  point  (2)  for  the  vibration 
form  with  the  corre- 

sponding components 


If  we  now  turn  to  the  general  calculation  of  the  field  of  a 
circular  piston  meobmne^V  wil^  resolve  the  integral 


into  its  real  and  imaginary  parts. 


We  therefore  obtain 


p = c • a 


^ Elolctr.  Naohr.  Teohn.,  vol*  12  (1936)  pp.  16-30. 


wherein 


Va 


Vm 


\t 


sin  hr 


dF, 


1 f cos  hr  , -r, 

jJ-y-dF 

F 


(121) 


This  is  advantageous  since  the  relative  pressure  components  va 
and  vm  exhibit  an  essentially  different  behavior.  If,  for  ex- 
ample, the  wave  length  is  so  great  that  2nrl?.<l  , pa 
becomes  very  small,  i.e.,  the  component  pm.  (which,  with  the 
velocity  amplitude  to  = wae'"n  possesses  a phase  displ '.cement  of 
90°)  predominates.  Physically,  this  means  that  the  membrane  works 
on  the  field  almost  wattlessly  so  that  it  moves  a dead  mass  back 
and  forth  without  radiating  practically  any  sound. 


lit  order  to  be  able  to  evaluate  the  integrals  (121),  we  must 
first  relate  the  variable  r to  the  corresponding  surface  element 

(determined  by  its  polar  coordinates  o,cp  ) 
If  we  denote  the  angle  which  the  field 
point  line  OP  forms  with  the  Z-axis  by 
y-  and  let  the  distance  of  the  field 
point  P from  0 be  rt  while  we  let  r be 
the  distance  of  the  field  point  from 
the  surface  element  dF,  then  one  has 


t — /rj  + Q * — Sr^coaijpsiny. 


(122) 


Fig.  69.  For  the  calcu- 
lation of  the  circular 
pieton  mesfcraaa. 


Here  the  field  point  P is  assumed  to 
be  in  the  X-Z  plane  in  Fig.  59  which,  due 
to  symmetry  conditions,  manifestly  does 
not  limit  the  generality.  The  integrals 
(121)  then  assume  the  fom 


p» 


I d<t  j QdQ 


Hill  hr 


2.. 


Pm 


I d<r /(,<!(> 


cos  It  r 


(123) 


where  the  value  of  r given  by  (122)  is  to  be  substituted.  Since 
an  integration  in  closed  form  does  not  appear  generally  possible, 
one  will  attempt  to  so  transform  the  integrands  by  a series  devel- 
opment that  the  variables  of  integration  appear  separated.  This 
is  done  by  a series  development  whose  terms  are  formed  from 
spherical  harmonics  and  Bessel  functions.  ,Jm  lias  the  following 
relations  ^ i 


1 V(2,i  -t  \)s^x) sn(tj) /*.(«•#),  (124) 

ji*+y,-2jryco»a 


Watson,  0.  N. i Theory  of  Bessel  Functions,  p.  366,  Cambridge 

1922. 
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(125) 


cos  | a;2 + y® — 2 xy  cos  § 
\'x2+y*— 2xycoa  & 


x,  (-«  + l)<5,(a;)Cn(y)P)1(co8^),  (igt/) 

n ~ 0 

CO 

^T(  2 « + 1 ) 0„  (x)  8n  (y)  p„  (cos  0) . (a:  2?  y) 

H**0 


(126) 


Here  P„(cos??)  is  the  Legendre  spherical  harmonic  and 

+ »(*).  <?„(*)  = (127) 


where  the  J's  indicate  Bessel's  functions.  These  functions  may 
be  relatively  simply  represented  for  small  values  of  a as 
rational  functions  of  x , sinx  and  cos  x . Thus 


S'o(x)  = sinx, 

„ . , sinx 
Sl[x)  — ■■  ■ — C08X, 


St(x)  = ljsinx-  ~ cos*, 

o , . /1«  0\  ■ 


C0(;r)  =•  cos.r, 

Cj(;r)  = sinx  -t-  ^ . 

C2(x)  = ^ sinx  + — ljcosx,' 

C3(x)  -(‘j-ljsin.r 


(128) 


wherein  between  the  Sjx)  and  CK(x)  the  important  relation 

S.W(lllt.W-S„iWe,(,r)=l  (129) 


exists. 

Using  the  series  in  (124),  (125),  and  (126)  we  can  now  calculate 
the  integrals  (123)*  We  first  carry  out  the  integration  tern-wise 
with  respect  to  <p  . Here  we  use  the  well-known  relation  from  the 
theory  of  spherical  harmonics: 


2;1  2. -I 

] 1\h* i (oosysinyJi/f^O,  / 2.-? /'^(eoHy).  t.  1 ,»  J 


0 

0 

From  this  it 

follows  that  all  odd  terms  drop  out  and  there 

\\  - 

results: 

<• ' 

V 

' J. 

V* 

i VI  *- 

, > (4h  -t-  I)/',.  (l»)  t'j,  (<<»';  )S.H{krl)  1 S.KU)ilt, 

* • i *>— < 

U31. 

•X 

« 0 " 

Vm 

J VI 

...  ' (4h  . 1 1 /'i K (t w y ) r ; , ( i ) / .S'; , (c ) •/ ..  . 

(132) 
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These  term-wise  integrations  are  permissible  as  long  as  the 
series  (124)  and  (125)  are  uniformly  convergent.  While  this  is 
always  the  case  for  (124),  it  holds  for  (125)  only  if 
That  is,  while  the  development  (131)  is  valid  without  any  limita- 
tion (i.e.,  for  arbitrary  positions  of  the  field  point),  (I32)  is 
only  valid  if  the  field  point  distance  is  greater  than  the  membrane 
radius.  We  must,  therefore,  seek  another  development  for  pm  for 
field  points  in  the  neighborhood  of  the  membrane.  In  order  to 
carry  out  the  calculation  when  ri  < Qi  , we  imagine  the  membrane 
of  radius  ei  divided  into  a smaller  circle  of  radius  ^ and  a 
circular  ring  of  the  width  ei  ~ri  . If  we  denote  the  two  re- 
sulting regions  by  f , and  Ft  , the  integral  over  the  surface  F 
breaks  down  into  two  summands.  The  development  (132)  is  valid 
for  the  first  summand  except  that  kri  replaces  the  upper  limit 
of  integration  kQx  . For  the  region  Ft , formula  (126)  is  now 
x,o  be  applied,  however,  so  that  the  following  development  is  valid 
for  Pm  at  a field  point  distance  smaller  than  the  membrane  radius: 


00 

pm  =•  + 1)  jP*»(0)  Pin  (cosy) 

1 n “0 

krt  kg, 

X + $>„(£>•,)  j <?*«(*)<*]. 

0 ’ kr, 


U33) 


For  a practical  application  of  the  formulae,  a tabular  com- 
putation of  the  functions  t x 

P*,  Sn,CH,  [SiH(x)dx,  {Cia(x)dx 

0 (133a 

is  necessary.  For  the  spherical  harmonics  Pn{ oosy)  , tables  for  the 
values  of  n from  0 to  20  and  for  y =0°, 5°.,.90°  are  given^; 

The  functions  5„(«).  Cn{x)  are  available  for  values  of  * smaller  than 
2 for  integral  values  x * 1,  2,  3.  ...10nP,  and  for  the  inter- 
mediate values  x ■ 2.2,  2.4,...»9.8v&  The  tables  for 

* X 

fsin(x)dx,  and  lctn{x)dx 
0 

are  given  for  n ■ 0,  1,  2,  ...  10  and  x * 1,  1.25,  1.5,  • ••  10  in 
the  appendix.  For  the  calculation  of  the  latter  and  for  small  values 


^ The  symbols  x and  y used  here  refer  to  eqs.  (124)  - (127), 
inoluaive,  and  are  to  be  distinguished  from  the  ooordlnatea  (x,  y) 
of  the  eurfaoe  element. 

^ Phil.  Trane.  Roy.  Soo,  Lond. , Vol.  203  (1904),  p.  100. 

Rep.  Brit.  Aaeoo.  Adv.  Sol.  1916,  pp.  97-l07j  1922,  po.  263-270. 
^ Hep.  Brit.  A6soc.  Adv,  Sci.  .1914,  pp.  87-102. 

& Elektr.  Naohr.-Techn. . Vol.  15.  (1938),  p.  73. 
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of  x > a series  development  with  respect  to  x is  advantageous. 
This  results  at  once  by  a ternwwise  integration  of  the  well-known 
series  for  S. ln{x)  and  Cu(x)  . Thus  it  follows  that 


x 

/-Sn  + * 

Sin(x)dx  = lVg7^  (4 


0 

X 

x 


{in+  1) 


+ 2 2(2»  + 4)(4»+3)  1 2.4.(2n  + 6)(4?i  + 3)(4»  + 5) 

[n  . , , 1.3-5...  (4 7i  — l) 

jCitt(x)dx=  — l 


)• 


(2^1  + 


+ i 


2(2n  — 3)  (4n  — 1)  r 2 • 4 (2n  - 5)  (4» 


+ ...'1 


(134) 

(135) 


For  the  larger  values  of  x it  is  more  advantageous  to  return 
to  the  simple  functions  Sn(x)  and  Cn(x)  . For  this  we  proceed  from 
the  equation"$7. 


V 

(x)  = /~(>.  f 1)  ( 1) 


r 


lv+i-  M-U 


(136) 


If  we  set  v =- n — l,  p — n,  a = m,  wo  then  obtain 


~2('m) (2n  + W2n  + 3)  • • • (4n  - 2m  - + » Jm_, ,_,(*) 

m -0 


(136a) 


From  this  we  find  by  integrating  term-wise  considering  the  re- 
lation 


Jxf+ljj)(x)dx  = xp  + xJv+  x{x) 


the  desired  relation 


(136b) 


(~1)”/C3n(x)<tx  = S0(x)  + ( J)  (2n  + 1)  Gf] 

-(J)(2*+  D(2«  +3)^  + ... 
-(-l)"(“)(2n  + l)(2>.  + 3)..  (4»  - . 


037) 


Nielsen,  N. » Ilandbuoh  der  Theorie  der  Zylinderfunlction,  p.  269, 
Loiptig,  1904. 
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-»  x 

tor  the  determination  of  fstn{x)dx  , we  proceed  from  the 
formula  ^ 0 


Jy(x) 


>_(x)  _ yij-i)*'-'  /' 

p!  (p  — *)!  \. 


( — l)**  ' I v * p «\  / 2 \« 


J v — ip  r 5 (*p) 


S " 0 


and  obtain  for  v — 2n  + i , p = n,  s — m 

, /i  , V>(-1)'  ’/a  + 5 + w-  l\/2\"'  T 
■/in  + ‘ (a:)  ~ n (n~  m) ! ( m ) (x  ) J~m  + » (x‘) ; 

Ml  - 0 

from  which  we  find  by  term-wise  integration  considering 


the  relation 


i 

x>>:-  i 


( - 1 rjs,  n (x)dx  = - C0(I)  + ( J ) (2  * + 1 ) S°y 


(137a) 


(137b) 


(137c) 


(2*  + l)(2«.+  3)6i,J)-  r 


(138) 


l),'("l)(2n  + l)(2»+3)... 


(4n-l)Sll.1(x) 


K. 


here  the  constant  K is  to  be  determined  so  that  the  right  hand 
side  of  (138)  vanishes  for  x * 0 . 


^ince 

Sm"  1.3. 5.’ -2.-1  ttw  *»0)  («&) 

we  then  obtain  the  value  of  K by  the  equation 

i+(;)e. -hi  .(;)«»■■ 


+ « >1I 


,,/  » ) (2 « + l)(2i<  -I  :t)  ..  (4»  - 1) 
i • ;t  • r> ...  [2«  - i) 


(139) 


For  the  expression  on  the  right-hand  side  we  find  the  value 


( l )M 


, 2 ■ 4 ■ li  ■ 2 ii 

l • 3 • :>"■  2 II  -■  I 


(139a) 


Therefore,  it  turns  out  that 


(-()“/. s\,  Jr)  (I/  { - 1)- 


2 ■ 4 • »l  • • • 2 h 
1 2 « - l 


~ O0(x)  -r(J)(2«  • l)^-  -•  l!<2*  i + 

ir( ” )(u*»f  + DC..  • 3).-.(4«  - 


(140) 


& Niolsen,  N, » Handbuoh  der  Thaorie  der  Cylinder funkt lira,  p«  269, 
Lelptlg,  1904* 
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We  will  first  apply  the  formulae  when  the  field  point  is  assumed 
on  (i.e.,  directly  in  front  of)  the  membrane.  Then  there  is  to  be 
substituted'.  ^ 

cosyr^  1,  p2„(i)--=  1 and  p2n(0)  = (_1)«L;3^-2«-i 

Six:  different  values  of  A-gj  are  used  {0 «5>  2,  4,  6,  8,  and  10)  while 
the  field  point  assumes  all  values  from  r1  = o to  r1  = p,  immediately 
in  front  of  the  membrane  for  each  value  of  ko1  . In  the  calculations 
we  make  use  of  the  fact  that  the  portion  due  to  the  interior  surface 
elements  (o  < rj  is  found  by  formula  (.115)  so  that  for  p„  and  pm 
we  only  have  to  add  the  parts 

lr~  2 (4»  + l)PtAO)PiMStulkrl)fsi^x)dx  (T40b 


CO  £ n 

(‘JH"  + i)PuAO)l\Al)S,n(krl)fb,M(lx 


(140c, 


due  to  the  ring-shaped  region. 


ft  ft  f»:  fc 


3.  Sound  pressure  amplitude 
- immediately  in  front  of 
circular  pioton  membrane 
(Bad tut  'h  ). 

3.  <(>■•-«,  3 kt>.  -4.  *.  ):j,k  _ 6 

a a.  *>.,  — la. 


From  tno  two  components  p‘ 
and  p,„  thus  obtained  p =-.  | -i-  pi 
is  then  calculated  and  a curve  is 
drawn.  %re  rjg,  is  the  ab- 
scissa so  that  the  abscissa  1 
corresponds  to  the  boundary  point 
of  the  membrane  (Fig.  60). 

One  sees  that  only  for  very 
small  values  of  to,  (*■/_>,  < A) 
is  the  pres, sure  amplitude  to  some 
degree  constant  ovor  the  whole 
surface  of  the  membrane  and  that 
for  larger  values  of  an 

ever  increasing  waviness  makes 
its  appearance  while  the  value 
at  the  center  continuously  os- 
cillates between  zero  and  two 
and  the  boundary  value,  always 
smaller,  fluctuates  about  the 
value  £ 


V For  the  condition  as  her*  s^tec , i.c,,  the  field  point  on  the 
membrane,  the  condition  would  soem  to  :*>  cosy  ~0.  fins  consider- 
ation applies  likewise  to  the  remainder  of  equations  (HOJ  «nd  to 
editions  { 1.10b)  nnc!  (140c), 


i 


As  an  example  of  the  computation  with  an  arbitrary  position  of 
the  field  point,  the  calculation  procedure  for  the  case  leg  = 10,  kr  — 6 
will  be  explained.,  First  the  quantities 


°2n  — (4n,  + 1)  -P2n(0)  Stn(kr)  I Sin(z)dx 

o 


(140d) 


are  easily  calculated  using  the  tables  in  the  appendix.  Thus  one 
obtains 


2 n = 

1 0 

2 

4 

6 

8 

•10 

= 

i -0,086 

+0.124 

i +1,785 

! -1,694  | 

+0,222 

-0,008 

From  this  we  find  the  value  on  the  normal  axis  (y  = 0): 


10 


+ 0,342. 

2n=*0 


(140e) 


According  to  earlier  considerations  (p.  70  ) this  value  is  given 
by  cos  kr  - 00a ^(kr)2  + {kg)2  . In  fact,  this  yields 


cos6  — cos  )/l36  = 0,960  — 0,618  = + 0,342  . 


(140f) 


This  gives  an  important  control  for  the  correctness  of  the  coefficients 
«sh  . <»e  also  calculate  the  coefficients  63h  for 

Pm  =2  b*«  (cosy) 

and  obtain  the  following  table: 


2n  — i 

0 

2 | 4 

I 8 

8 

^ kr  \ 

(4  it  + 1)  P,,(Q)  C,,^r)  fSt.Mdx  f 0,006 

1 

+ 1,317  i +0,209 

I 

-0,631 

+0,287 

(4m  + i)g,.(0)  1 

+0,012 

! 

1 

i 

1 

4*0,093 1 — 1,727 

J -0,421 

+0,640 

kr 

i 

! 

lh  % 

+0,018 

+ 1,410  -1,518 

-1,052 

+0,927 

2«  -7  "1 

10 

12  14 

10  i 

18 

1 i ! 

1 

Compared  with  the  coefficients  aln  , the  last  terms  of  the  two 
series  here  are,  say,  larger  by  a power  of  10  so  that  we  have  to 
expect  a greater  deviation  of  the  sum  2\>„  from  the  true  value  y,„ . 
In  order  to  obtain  a better  control  with  the  directly  to  be  found 
values: 


k r 

siny'2(fcrj2  — sin kr  = ^ (4;.-  + l)PJn(0)  j Sin(x)dx 


(140g 


kn 

sin  )l(krf+{k^~  sin  f2(krj*  = £ (in  + l)P.,fl(0)  fc,n(z)dz  (tyOh 


we  observe  that  the  summands  of  the  two  series  in  the  tables  for  the 

larger  values  of  n oscillate  less  and 
- less  about  a mean  value.  life  take  account 

of  this  by  dividing  the  last  term  in  the 
two  series  by  2 and  therefore  put  the 
value  “0.031  in  place  of  -0.061  and  the 
value  -O.O38  in  place  of  -O.O76.  If 
we  denote  the  thus  obtained  sums  by  2i 
and  V„  we  then  have  V,  = 4- 1,091 
and  2a  = -1.591,  while  the  true  values 
give 


sin  y;72  - sin  6 = 1,087 


(0401 


w so  so  x x x* 
/— 


sin]' 136  — sin  p2  — — 1,595  (,140j 

as  a result  so  that  a sufficient  agreement 
obtains  hero.  Of  courso,  the  device  can- 
not be  applied  generally  (y  p 0)  for  tne 

calculation  of  pm  - 2 4>„  y) 

0 to  10  x w so  so  so  x so* 

since,  due  to  the  factor  P,„(oosy) 

ig.  61,  Sound  pressure  am-  regular  oscillations  are  no  longer  present 
plitude  (t>/«o  of  the  circu-  But,  even  here  one  will  be  able  to  count 

lar  piston  membrane  with  on  a deviation  of  at  most  10/o  from  the 

the  radius  y.  (*». -10)  and  true  value.  For  acoustic  calculations, 

constant  field  point  however,  this  accuracy  is  in  general  to 

distance  r(*,“w-  be  regarded  as  quite  sufficient,  being 

tho  tables  for  the  spherical  harmonics,  the  quantities  • P,*  (oiw y) 
and  PlH( cosy)  are  then  found  for  y - 5",  io°  etc.  There  thus 

results  the  following  table  iHere,  the  values  given  for  p„  and  p„ 
are  derived  from  computations  originally  carried  out  to  four  decimal 
places,  therefore  the  values  of  p„  and  pm  deviate  at  tidies  in 
the  last  decimal  place  from  the  sum  of  the  numbers  standing  above 
these  values): 


plitude  (vito)  of  the  circu- 
lar piston  membrane  with 
the  radius  ?,  <*»,-io>  end 
constant  field  point 
distance 


-0,086  -0,086 
+0,124  +0,122 
+ 1,785  +1,717 
-1,604  -1,560 
+0,222  +0,193 
-0,008  -0,007 


-0,086  - 0,086  -0.086  -O.iKv,  - -i.iiSfi 

+0,1 18:  +0,1 1 1 +0.102  +(i,,!:":  -.-.-•77 

+ 1,639: +1^22  + 0.847  -0.!'- 
— 1,191  j— 0,675  -0.122  +0,8i.,.  n.t.34 

+0,116^0,021  -0,056  -0,000  - 0.076 
-0,003'  +0,001  +0,003  +0.003! -0.001 


- 0.086  -0,086  -0,086 
+0,063  +0,047  +0.031 
-0,305  -0,569  -0,725 
+0.895  +0,548  +0,251 
-0,026  +0,031  +0,066 
-0.002  -0,002  -0,001 


+0,342 . +0,379  +0,493  +0,597  +0,688  +0,705; +0,593  ,,-,^  -0,031  -0,463 


7 =•  50°  y — 55° 

y — (iO'  7 uf>°  J 7 

)'  --  80“  )■  r. 

0 —0,086  ; —0,086  ■ —0,086  - 0,088  - 0,086  -0,086  - 0,086  - 0,086  -0,086 

2 1 +0,015  j -0,001  —0,016  -0.029  - 0,040  -0,049  - 0,056  -0,064  - 0,062 

4 ( -0,762  l -0.686  -0.516  -0,277  -0,007  +0,256  +0,474  ; +0,618  +0,669 

6 • -0,095  j -0,380  -0,548  -0,531  -0,354  , -0,073  +0,223  +0,446  +0,529 

8 +0.085  1+0,032  : -0,016  - 0,053  -0,062  ■ -0,038  +0,005  + 0,045  +0,061 

10  . +0,001  I +0,002  : +0,002  ; — : -0,002  ! -0,002  -0,001  +0,001  +0,002 


+0,019  j + 0,019 
+ 1,410  —1,393 
-1,518  +1,460 
-1,052  - 0,968 
+0,927  +0,804 
-0,166  -0,373 


4 

— 1,518 

6 

-1,052 

8 

+0,927 

10 

— 0. 166 

12 

+ 0,265 

14 

-0.214 

16 

+ 0,107 

18 

-0,137 

+0,019. +0,019  +0,019 
+ 1.345  +1,265  +1,160 
-1,2931 -1,038' -0.721 
-0,740; -0,41 8| -0,075 
+ 0,484!  +0,089  -0,234 
-0,149  + 0,077' +0,187 
4-0,035:  -0,100j  — 0.104 
+0,014  +0,086  +0,034 
- 0,036  — 0,060'  +0,013 
+0.045  + 0,031 1-0,035 


+0,019  +0,019  +0,019+0,019:  +0,019 
+ 1,030  +0,880  +0.714  +0,535  +0,355 
-0,375  -0,035  + 0,270  + 0,484  + 0,616 
+0,214  +0,393  +0,432  +0,340  + 0,156 
-0,370  - 0,314  -0,107  + 0,129  +0,276 
+0,142  + 0,003  - 0,118  --0  138; - 0,054 
- 0,007  + 0,080  + 0,074  - 0,006 : -0,073 
0,050  - 0.055  + 0,012  +0,055,  4 0,022 
+0,056  +0.001  -0,043  -0,011  +0,036 
-0,021  -» 0,031  +0,016  0,028.-  0,012 


-0,479  - 0,276  -0,049  + 0,244  + 0,626  +1,003  +1,269  +1.370  +1,341 


+1.066  1+0,564  0,(107  -0.533  0.964  - 1,186  | - 1.032  - 0,503  1 
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Fig.  62.  Sound  field  (/>,<•  a) in  the  neighborhood  of  the  circular 

piston  membrane  '.«e, 
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1-.  ^ig.  61  the  two  components  pa  and  pm  are  first  repre- 
sented individually  since  a valuable  control  for  the  results  of  the 
calculation  is  furnished  by  the  continuous  sequence  of  the  calcu- 
lated values.  Finally,  the  complete  behavior  of  the  relative 

pressure  amplitude  v — \i>i  + pi,  is  obtained  by  geometric  addition 
of  the  corresponding  values  of  va  and  p,„  . (See  Fig.  61.) 

In  this  manner,  the  relative  amplitudes  fp*  +pl,  were  calculated 
for  three  piston  membranes  (whose  radii  were  given  by  * 4, 

6,  10)  for  such  a great  number  of  points  of  the  nearby  field  that 
the  curves  of  constant  pressure  amplitude  could  be  drawn  by  inter- 
polation. The  results  are  represented  in  Figs.  62  to  64. 


Fig.  63.  Sound  field  in  the  neighborhood  of  the 
circular  pieton  membrane 
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Here  the  numbers  on  the  horizontal  axis  aro  the  values  *\<» , , so 
that  the  points  which  are  the  length  of  the  membrane  radius  distant 
from  the  membrane  center  lie  on  tt.o  unit  semicircle,  dno  recognizes 
that  the  essential  character  of  the  field  is  determined  by  the  loca- 
tion of  the  zero  and  maximum  values  (which  are  found  only  on  the 
normal  axis).  If  we  first  oonsider  the  case  !•(»,  u,  we  see  that 
the  amplitude  has  approximately  the  value  zero  at  the  center  of  the 


—88— 


gSi  ■*: 

Hfr 
Sf  ‘ 


membrane  and  that  on  the  normal  axis  with  increasing  distance,  the 
amplitude  first  increases  rather  rapidly  and  then  rises  more  slov/ly 
to  the  value  two,  and  falls  off  very  gradually  from  there  on.  Off 
of  the  normal  axis  lie  two  peaks  at  the  height  of  the  center  which 
are  marked  by  the  high  pressure  lines  1.4.  Beyond  these  peaks  a 


Fig.  64.  Sound  field  ip.o  in  the  neighborhood  of  the 
circular  pie  ton  membrane  <**-io>. 
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uniform  drop  takes  place  vdth  increasing  radius.  If,  for  comparison, 
we  now  consider  the  field  for  *•<?,«  10  , a striking  similarity  of 
the  configuration  of  the  upper  part  of  this  sound  field  with  the 
foregoing  appears.  If  we  imagine  a horizontal  straight  line  drawn 
through  the  zero  value  on  the  normal  axis,  we  see  that  this  cuts 
the  field  into  two  parts  of  which  the  upper  part  is  extraordinarily 
similar  to  the  wide  field  for  hh  ^ « as  is  shown  in  Fig.  63.  The 
lower  part  of  the  field  kgx  <■-  U>  is  essentially  different  since  now 
a minimum  and  a maximum  have  made  thoir  appearance  both  to  the  right 
and  left  of  the  two -point . 

If  we  r»w  allow  the  radius  of  the  piston  membrane  to  increase 
continuously,  we  then  know  (from  the  simple  formulae  for  points  on 
the  normal  axis)  that  the  maximum  and  zero  values  become  increasingly 


I* . 


-89- 


more  numerous  and  progress  upwards  while  new  zero  and  maximum  values 
continually  make  their  appearance  at  the  center  of  the  membrane. 
According  to  the  above  considerations,  however,  beyond  this  we  have 
to  expect  that  the  constant  pressure  curves  lying  to  the  right  and 
left  (of  the  center)  also  simply  move  upwards  with  the  upwardly 
travelling  zero-  and  two-points  of  the  normal  axis  without  an  essen- 
tial change  in  their  character.  If  we  have  a membrane  with  an  ar- 
bitrarily great  radius  and  mark  the  two  last  two-values  of  the  ampli- 
tude on  the  normal  axis,  we  will  then  be  warranted  in  expecting  that 
the  accompanying  field,  which  lies  between  the  planes  parallel  to 
the  membrane  through  these  two  maximum  positions,  agrees  in  charac- 
ter with  the  field  represented  in  T'ig.  64.  Moreover,  according  to 
Fig.  64,  one  would  suppose  that  upon  travelling  from  any  maximum 
position  on  the  normal  axis  in  the  vertical  direction,  just  as 
many  extreme  values  will  be  encountered  as  if  one  travelled  from 
the  same  maximum  position  in  a horizontal  direction. 


Part  Three 


THE  SOUND  FIELD  OF  THE  SPHERICAL  RADIATOR 
5.  The  simple  spherical  radiator  of  definite  order 


An  essential  assumption  for  the  radiators  previously  discussed 
was  the  existence  of  a sound  reflecting  rigid  infinite  plane  wall. 

In  its  position  of  rest  the  radiating  surface  coincided  with  this 
wall.  The  calculation  of  the  radiation  process  was  then  accomplished 
by  the  calculation  of  an  integral  over  the  radiating  surface.  For 
the  problems  treated  in  the  following  chapters,  the  radiating  sur- 
face forns  a part  of  a sound  reflecting  rigid  sphere  of  fixed  radius. 
here  exists  a more  general  problem  since,  aside  from  the  amplitude 
and  extent  of  the  radiating  surface,  the  extent  of  the  rigid  wall 
can  also  be  changed,  "hile  formerly,  for  the  oircular  piston 
membrane  e.g.,  one  characteristic  function 


sufficed  to  completely  specify  the  sound  field  for  an  arbitrary 
field-point  position  at  a great  distance  and  an  arbitrary  membrane 
size  (in  relation  to  x ) this  is  now  no  longer  possible.  On  the 
other  hand  there  are  now  quite  definite  distributions  of  the  ve- 
locity amplitude  for  which  the  solution  becomes  conceivable.  Here 
it  is  also  particularly  important  that  the  calculation  of  the 
neighboring  field  offers  no  difficulty . This  simple  solution 
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exists  if  the  velocity  amplitude  is  given  by  one  spherical  harmonic 
of  definite  order.  This  assumption  of  the  velocity  distribution, 
which  at  first  appears  artificial,  immediately  finds  its  justifica- 
tion in  the  fact  that  an  arbitrarily  given  velocity  distribution 
can  generally  be  reduced  to  that  given  by  a spherical  namonic. 

In  the  computation  of  practical  problems  we  will  limit  ourselves 
to  those  processes  which  are  rotationally  symmetric  with  respect 
to  an  axis  through  the  center  of  the  sphere.  The  corresponding 
spherical  harmonics  are  then  given  by  the  well-known  Legendre 
functions 


PJft),  p = cosy  P0(n)  = l,  PM  ft,  PM  = -2:r [ft2  — yj etc , 


(140k) 


By  a spherical  radiator  of  zero  order,  we  will  understand  a 
pulsating  sphere  for  which  the  individual  elements  of  the  surface 
vibrate  cophasally  outwards  and  inwards  with  constant  velocity 
amplitudes.  (See  *'ig.  1-a.)  The  velocity  on  the  sphere  is  given  by 

ro  = w>0e'wi.  (14Cto) 


With  a sphere  radius  r0  and  for  a field  point  distance  r , the 
sound  pressure  is  completely  determined  by 

Po  = Pt0) «.  I “ C lw(+  "/2  - '•»,  (141) 


where,  for  abbreviation, 


• (T  • F 
Mr~ 


(241a) 


From  this  follows  the  pressure  amplitude 


V o 


P(°)w0  -pi  - 

y °|'l  4 k*r‘ 


cyn  o 
' W j/  iTi»r5  ’ 


042) 


For  kr0  < 1 wo  obtain  the  formula  given  previously: 

Po  — “oT.V  u°  • (142a) 


The  relative  pi'essure  amplitude,  VqIc<jwq 
given  by 


Po  = i **>;! 
oon  ir’j'i 


is  then  generally 


(142b) 


The  curves  of  constant  pressure  amplitude  are  therefore  very  simply 
represented  in  the  spatial  sound  field  by  concentric  spheres. 
long  as  kr0<l  , the  pressure  increases  quadratically  with  kra  . 

As  long  as  kr0  > 1 , it  increases  linearly  with  ;-r0  . 
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The  simple  spherical  radiator  of  higher  (nth)  order^ls  charac- 
terized by  the  fact  that  the  variation  of  the  velocity  amplitude  to, 
is  given  by 


(142c 


The  pressure  pn  at  the  field  point  determined  by  the  polar 
coordinates  r,  y is  then  generally  given  by 

p„  = + (143) 


Here 


Fn(ix)  = (1  + ix)fn(ix)  - ixj'n(ix) 


(144) 


and 


...  . 7 w(n+l)  . (»— l)(»)(n-H)(n+2)  1 • 2 • 3 • • • 2w 

/.F)  — ‘t- 2ix  +"  2 • 4 (»*)*  1 2 -4 *6  ■ • • 2 »(»*)"  ’ 


(145) 


Tt  is  better  to  introduce  the  Bessel  functions  Sn{x)  and  Cn(x) 
defined  by  (127)  page  79  in  place  of  Fn(ix)  and  /„(tx)  . Bue  to 

the  relations 


in+1e-ixfn(ix)  = Sn(x)  + tCB(x), 

in+1e~izFn(ix)  = xSn+1(x)  — nSH(x)  + i[xCn+l{x)  — n<4(*)]. 


(145a, 


we  than  obtain 


P„  = p101.  + 


S^kr)  + \C,(kr)  p . 
V.(Fr0)+  i V.W,)  • 


where,  for  abbreviation, 


Un(x)  = xSn+i(x)  - rt  S„(x), 
F„(x)  = a :CB+,(x)  - nCn(x). 


(146) 

(147) 


Then  one  has 


I/fl(x)  sc  ein  x — x cos  x , 

(7,(x)  = (2/x  — x)ainx  — 2cosx 


V0(x)  — cos  x + x sin  x , (147a! 

Fj(x)  = 2sinx  + (2/x  — x)  cos  x etc. 


Hor  the  spherical  radiator  of  the  first  order  the  velocity  ampli- 
tude h>,  is  given  by 

tn,  — u\Pl(n)eiul  ~ u',cosy  • ( 147b, 


ASh  ayleigh:  Theory  of  Soupd,  Seo.  323  and  following. 


This  corresponds  to  the  motion  of  a sphere  vibrating  to  and  fro 
as  a rigid  body.  Here  the  pressure  is  given  by: 

-5-r  -cosy  (147c) 

170) 

while  the  pressure  amplitude  px  is  given  by 


kr 


2+  t^ir0 


Pi 

caw1 


»rj 

|/4  -|-  &r 


•^>'1  + Pi  cosy. 


U47d) 


As  long  as  /fcr0«l  , the  pressure  amplitude  increases  propor- 
tionally to  k*rl  . As  long  as  kr„  » l , the  pressure  amplitude  in- 
creases proportionally  to  kr0  . If  we  denote  the  maximum  pressure 
amplitude  occurring  at  the  field  point  r — r0  and  y = 0 by  pm  and 
require  the  curves  of  constant  pressure  amplitude 

Pi  = “Pm  («  = 1 ; o.9 ; o,8;  etc . ) , (I47e) 


it  then  follows  that 


Pi 

Pm 


(147f) 


If  we  choose  kr0  = i as  an  example,  then  for  the  field  points 
on  the  symmetry  axis  (y  = 0)  there  results: 


a 


-!- . } 1 h + X 

\r2  kr  | T 


(147g) 


and  from  this: 

kr  - AyT +7T+ 8““ 


(147h) 


Accordingly  one  obtains  for  the  values  of  Pi/Pm  on  the  normal 

axis 


PllVm  “ 

0,1 

| 0,2 

0,3 

0,4 

0,5 

0,8 

0,7 

0,8 

0,9  | 1 

rjr0  » 

7,14 

I 3,87 

2,53 

1,98 

1,05 

1,44 

1,28 

1.17 

1,07  1 

If  we  now  imagine  the  circles 

r -=■  r9,  r = 1.07  r0,  r = 1,17  r0,  r = 1 ,28 r0  , etc,,  (1471) 

constructed,  by  calculating  the  y-  values  we  can  easily  state  the 
points  for  the  required  values  of  Pi lPm  . *or  example,  on  the  circle 

r * 1.28  r0  , the  values  of  cos  y for  the  corresponding  values  of 
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Pl/Pm  are: 


%■ 

) 

II 

0,7 

0,tf* 

0,5 

| 0,4 

0,3  | 0,2  | 0,1 

1 

6 

5 

1 4 ; 

3 .!  2 ; 1 

cosy  =s 

7 

7 

1 i 

7 j 7 j 7 

As  an  additional  example 

i we 

give 

the  corresponding  tables  for  kr0 

PllVm  = 

0,1  | 0,2 

| 0,3 

j 0,4 

! 0,5  1 

1 

\ 0,6  1 

0,7  | 0,8  ; 0,9  | 1 

*■/»■«  = 

9,81  j 4,91  | 3,28 

2,16 

1,97 

i 1,65  | 

1,42  | 1,24  i 1,11  | 1 

For  kra>i  , the  corresponding  values  of  r/r0  are  simply  given 
by  l ~ . therefore,  in  thin  case 


Pi/P»  = 0,1 

0,2 

0,3 

0,4 

| 0,5 

0,6 

1 °‘7 

| 0,8 

0.9 

1 

r/r0  ==!  10 

1 6 

3,33 

2,50 

| 2,00 

| 1,67 

1,43 

1 1,26 

! 1.11  1 

1 

In  Fig.  65  are  drawn  the  curves  of  constant  pressure  for  the 
spherical  radiator  of  the  first  order  (with  kr0  ~ 1 ) . Because 

of  the  given  antisymmetric  distribution  of  the  velocity  amplitude 
with  respect  to  the  equatorial  plane,  the  character  of  the  field 
is  essentially  different  a.3  compared  with  that  of  the  2ero  radiator 
since  the  pressure  is  nov.  zero  in  the  equatorial  plane. 


From  (143)  the  pressure  amplitude  p„  for  the  spherical  radi- 
ator of  the  nth  order  is  given  by 


Vn  = P°  • 


• «i !(Vr) + <’*(*r) 
I'llkr*)  +yUkr„) 


Pn  (cosy) , 


(147 j) 


where 


n v*  j 

L°  J 

L L...J 

1/x2  )-  1 
x 2 + 4/a-2  | 

2 

3 

#’.(*)  < ClM 
Ul(x)  + V\(x) 

1 

l + x* 

(t  - 3/.r2)2  4 (3/a-)2 
1(4  -W  + (i  • 9/af 

(1  - 15/a.1)2  + (6/x  + 15/a-2)2 
(7  - 60/x2)2  + (x  - 27 lx  60/a:5)2 

Tt  is  noteworthy  that  the  directional  influence  - i.e.,  the 
ratio  of  the  pressure  amplitude  at  an  arbitrary  field  point,  ( »’.  y) 
to  the  pressure  amplitude  at  the  field  point  at  the  same  distance  r 
on  the  symmetry  axis  ( r,  y --  0 ) - is  given  independently  of  r 
simply  by  /’.(wwy) so  tliat  the  behavior  of  the  sound  field  is  obtained 
without  difficulty  from  the  values  on  the  symmetry  axis.  In  partic- 
ular, the  rttdal  lines  /\,(o<>sy)  0,  at  which  the  velocity  amplitude 


1,0 


%S  3,0 


IS  1,0 
r 

n — 


Jig.  65.  Curves  of  constant  sound  pressure  for  a 
radiator  of  the  first  order. 


ft 


on  tho  radiating  sphere  is  equal  to  zero,  determine  quite  generally 
the  conical  surfaces  in  space  for  which  t he  pressure  amplitude  in 
space  vanishes. 


Calculations  for  spherical  radiators  of  a higher  order  will  not 
be  carried  out  since,  from  a practical  standpoint,  scarcely  any  cases 
occur  in  which  a radiator  behaves  as  a single  radiator  of  higher  order. 


6.  '■‘■'he  compound  spherical  radiator 

Important  for  the  following  is  the  investigation  of  the  sound 
field  which  is  produced  by  a combination  of  a series  of  spherical 
radiators  of  different  order.  For  a sphere  with  tne  radius  r0, 
whose  velocity  amplitude  to  is  represented  by  a superposition  of 
spherical  harmonies 


TO  = [.40P0(cos  y)  + .4,  P,  (cosy)  + •••  + AnPn{oos  (146) 

it  turns  out  that  the  resulting  pressure  is  given  simply  by  the  cor- 
responding superposition  of  the  pressure  due  to  tne  individual  spher- 
ical radiators.  Tno  pressure  p at  the  field  ;oint  (r,y)  is  then 


^ = p(u)  gi  i + 


III  rJ  0 


_S»jk_r)  + iCn(kr)  • 
Um(icr0)  + i Vm(kr0) 


AmPm(  cosy). 


(14'-;) 


Here  the  Am s are  constant  quantities  which,  in  a given  Case, 
may  be  complex  (in  order  to  take  account  of  different  phases).  The 
principal  significance  of  this  theorem  lies,  however,  in  its  revers- 
ibility. In  general  the  above  velocity  distribution  does  not  exist 
in  the  form  (146) • The  essential  thing  now  is  shat  by  the  develop- 
ment theorem  of  spherical  1 '.a monies  an  (to  within  general  continuity 
requirements)  arbitrary  function  in  the  form  (146)  can  be  represented 
on  the  spherical  surface. 


Thus  the  calculation  of  the  sound  field  for  an  arbitrarily  given 
velocity  amplitude  on  the  sphere  is  then  possib-e.  he  will  illustrate 
this  in  a special  example. 


On  the  part  of  a rigid  immovable  sphere  which  is  characterized 
by  0 >,  y h.  y0  (a  spherical  cap),  a constant  velocity  amplitude  «• <-  1 
is  given  while  for  the  remaining  part,  «■  -■  o (Fig.  66).  ho  first 
have  the  problem  to  determine  the  coefficients  .4,,,  of  the  function 

(150) 


«’  - A0 P„(eo» y)  + .4,  P,  (cosy)  + ... 


so  that 

w = lfor  0°gy^y0, 
w — Ofor  y0  < y SI  180° 

Here  the  orthogonality  relations  for  spherical  harmonics 


(150a) 


Pig.  66.  The  radiating 
spherical  cap. 


and  thus  generally, 


+ 1 

Jpm{/j)Pn(/i)d(i  = o,  for  w 4=  n (150b) 

+ 1 

jpU^dfi  = 2»Vr 

-i 

become  useful. 

When  we  multiply  both  sides  of  equation 
(150)  by  P«{p)  and  integrate  from  *1  to 
■fl , we  immediately  obtain 

+ 1 +1 

I wPn(n)dp  = A„J Pl(/i)d/j,  = An  • (150c) 

-l  -i 


An  = —^~fu'Pn(fi)dn.  (I50d) 

-1 

Now  in  our  example  w is  only  different  from  0 in  the  region 
008  y0  Si  ,u  £3 1 so  that 


1 

-4,  = --f-f  P,An)dfi  (150e) 

COB  Y l 

By  the  well-known  relation  of  tho  spherical  harmonics 

(2tt  + \)PH(x)  = -j-  [P#+l(^)  P d-lWl  (151) 


it  follows  that 


(2)i+lj  jPn(n)dn  =>  ?».,(*)  - P^(x),  (n^  1) 


jP^dfj,  - 1 — x. 


(152) 


Thus  there  results  the  required  representation  for  wt 

00 

W - [P«.i(oo8y0)  - P*  + 1 (008  y0)j+.(0O8y)  (153) 

H “ J 

and  therewith  also  the  required  pressure  by  (149)  upon  substitution  of 
the  resulting  valuos  for  ^4.. 
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If,  as  a special  case,  we  choose  a radiating  hemisphere  with  a 
constant  velocity  amplitude  w = 1 which  is  completed  by  a rigid 
hemisphere  on  which  *r  » Q , we  then  have  to  set  cosy0  = cos90o=0  • 
•^ince 


p8n(0)  = (—iy> Li--5  "l2”-!)  (153a) 

and  i>2n+i(0)  = 0 j we  lind: 

w = ^ + | A(c°sy)  — |-|p3(cosy)  + ~-j|p6(oosy)  + •••.  (153b) 

and  for  a radiating  spherical  cap  with  y0  = 60°: 


w = 0,260  + 0,663  Px{ cos  y)  + 0,469  P,(cos  y)  4-  0,082  P8( oos  y) 

- 0,264  P4(cos  y)  - 0,306  P5(cos  y)  - 0,067  PJcos  y) 

+ 0,198  P 7 (cos  y)  + 0,245  Pg(cos  y)  -f  0,057  P,( cos  y)  — • • • 


(153c) 


In  Figs.  67  and  68  the  corresponding  approximation'cuu  ves 

m — n 

U»(")  = ^ amP„  (cob  y) 
m-0 


(153d) 


are  represented.  The  deviation  of  the  approximation  function  from 
the  theoretical  valu«  unity  may  be  determined  for  y = 0°  , Since 
Pm(l)  = 1 in  wMt  all  terms  except  the  first  two  and  the  last  two  then 
cancel  out  so  that 

<_0.,  - 1 - (153e) 


In  the  present  case,  y0> 
equal  to 


00° 


, the  error  therefore  appears 


P,-.(Q)  + P.(Q)  t 
2 


Therefore  for  n » 0,  1,  2,  . . . etc.  it  equals 


+ i> 


. _L . i.  + — • — etc , 

2 2*2  2-4 


In  order  to  obtain  a bettor  agreement  with  the  first  formula, 
we  will  introduce  the  radiating  surface  F — 2 rcrj(l  — oosy0)  for  the 
representation  of  the  radiating  spherical  cap  defined  by  y0  and 
having  the  constant  velocity  amplitude  ir  , and  we  will  write  the 
pressure  p in  the  form 

p » Cofeu««+«m)Ll^ . fo)  /i>  p0).  (153f ) 

C ■ of  • IV 

Here  A>~  2>  is  the  pressure  amplitude  if  the  dimensions 

of  the  radiator  are  small  compared  to  the  wave  length  and 


i!(r*w<.P  o) 


S'(tr)  + '(\(kr) 


-I 


«x> 

‘po  V/y/dlPn  lW-P^.t/'o)] 

M * l 


A'.U-r)  4 iC,{kr) 
P.tJfol-HT.tir,)' 


(154) 
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Fig.  68.  Approximation  of  the 
function  ^ 

spherical  harmonicB. 


A casa  of  practical  importance  which  has  been  treated  by  Rayleigh 
exists  when  the  radiating  surface  is  very  small  so  that  N = cos  y0 
can  be  replaced  by  unity.  Then  formula  (154)  becomes  inapplicable. 
However,  we  can  easily  remove  this  disadvantage  when  we  consider  that 
by  (151)^>: 

- 1 t“o)  -Pn  + 1 (/*<))]  j = 2m  + 1 

He  then  obtain  for  a point  radiator  on  a rigid  sphere 


(154a) 

(155) 


where 


$(r,r0liu,  1)  = V(2n+ l)Pn(/i) 


«0 


Sn(kr)  + i C„  [k  r)_ 
U.(kr0)  + iV,(kr0) 


(155a) 


A substantial  simplification  of  the  formulae^  occurs  if  the  field 
point  is  at  so  great  a distance  r that  one  can  set  fn(ikr)  ~ l , 
bince 


SH(x)  + iCn(x)  = i»+‘e-‘'  /„(»*) 
Sn(kr)  + iCn(kr)  becomes  in  + le-‘*'. 

He  then  obtain 

where 


(155b) 


(155c) 


°o 

f(vv\= - - u 1 ~y.n  + 1 [fi.  - 1 (ft  o)  ~ P*  + 1 (/VI  /un 

nY’Yo)  f4(ir0)  + .K0(ir0)  1-^Z/  t/,(*r0)  + i F„( jfcr,)  (15b) 

ii“  l 

For  practical  calculation  it  is  advantageous  for  the  summation 
to  group  tl.  berms  with  the  even  n and  with  odd  n when  we  write 


Hy<  yo)  = 


+ T— r@i  + 


(157) 


Uo{kro)  *f  * y o(^*ro)  1 *”  /*o 

whore  <3j  contains  the  spherical  harmonics  /*,(/»), /*«(/<),.. . etc.  and  ©8 
contains  the  spherical  harmonics  l\ (ft),  P3(fi),  ...  etc  • 1'hen  since 
/’■•(-/<> -p„ut)  and  Pj.H(-rt- -p„n(h)‘- 


1 --  fid 


©a 


and 


/(y.  180°  — y0) 


4 + '4  [ Ho 


1 + 


(158) 

(159) 


Eq . (151)  is  to  bo  used  when  L’Hosni  tal's  rulo  is  a;- died  in 
evaluatin'1;  the  indetormi  ns  *.e  form  on  tho  loft  of  ( 1 5 'in  ) . 

^Formulae  (I53f)  and  (154)  are  referred  to  hero. 
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L&nifestly,  corresponding  relations  are  also  valid  for  field 
points  which  do  not  lie  at  a great  distance.  We  need  therefore 
oniy  to  oarry  out  the  calculation  for  0^ys;90o,  0^y0si90o  and 
the  values  for  90°  ssys;  180°  and  90°£y0£l80o  then  result  from  (158) 
and  (159). 


Physically  (156)  and  (157)  signify  that  one  can  replace  an 
arbitrary  velocity  distribution  on  the  spherical  surface  by  a dis- 
tribution symmetrical  to  the  equatorial  plane  and  one  antisymmetric 
to  the  equatorial  plane.  The  first  contains  the  spherical  tiannonics 
of  even  order  while  the  second  contains  the  spherical  harmonics  of 
odd  order.  By  (158)  one  finds  the  pressure  at  the  image  of  a field 
point  when  one  inserts  the  antisymmetric  part  with  the  opposite 
sign  in  the  formula  for  the  field  point.  Furthermore,  it  follows 
from  (158)  and  (159)  that 


(1  - /i0)/(180°  - y , y0)  + (1  + /*0)/(y,  180°  - y0)  = (l59a) 


This  says  that  the  radiations  due  to  (158)  and  (159)  together  yield 
the  radiation  of  a spherical  radiator  of  zero  order. 


as  an  example,  in  the  table  on  page  102,  the  expressions 


Un(krB)  + i F.(Jfcr0)  [p»(coayo)  Pn+i(cosyQ)]  = A„  + iB„ 


(159b) 


are  given  for  kr0  = io  and  y0  = 0°, 


10°,  20°,  ...  90 u and  the  correspond- 
ing values  /(y,y0)  for  y = 0°  ^nd 
y = 180° are  calculated. 


to  k u to  to  m m m 
x — 


Fig.  69.  Pressure  amplitude  &■ 
a function  of  the  central  angle 
of  thn  radiating  surface  with 
constant  deformation  volumes 
1.  for  the  principal  direction 
>r-<».  2.  for  the  opposite  direc- 
tion 


In  Fig.  69,  the  magnitudes  of 
/( 0°,  y0)  and  /(180°,  y0)  are  represented. 
These  curves  give  a conspectus  of  the 
pressure  amplitudes  in  the  principal 
■direction  (y  = 0°)  and  in  the  opposite 
direction  1 (y  = 180°)  as  a function  of 
the  size  of  the  radiating  surface 
(y«)  with  constant  deformation  volume 
(F  • w)  . ,.e  immediately  recognize 

that  for  small  values  of  Yo  in  the 
principal  direction,  the  double  amp- 
litude occurs  as  compared  with  the 
amplitude  for  nondiroctional  radi- 
ation / c*aF*«’\ 

~ 2 Ar  I' 


If  one  desires  the  largest  possible 
radiating  surface,  the  greatest 
possible  pressure  amplitude  in  the 
principal  direction  and  simultaneously 
the  greatest  possible  extinction  for 
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the  opposite  direction,  the  figure  shows  that  this  is  attained  for 
7o  m 20°  approximately. 

With  the  aid  of  the  tables  for  An  + iBn  , it  is  then  easy  to 
calculate  the  complete  sound  field  (for  a great  distance  of  the 
field  point)  for  every  value  of  y0  . One  needs  only  to  compute 
the  expressions 


; + * J — f— • 'E  (An  + i B„)  Pn  (cos  y) 


(159c) 


and  to  plot  }?A2  + Ji  as  a function  of  7 . * w components  h,  J 

corresponding  to  the  values  y0  — 10°,  30°,  60°,  90°,  120°,  150°,  and  170°  are 
plotted  in  Figs.  70  to  76  and  the  relative  pressure  amplitudes  ob- 
tained therefrom  are  represented  in  Figs.  77  to  83. 

As  long  as  the  radiating  surface  is  small,  one  may  expect  a 
directional  effect  similar  to  that  of  the  piston  membrane  with 
increasing  ya  « For  the  larger  values  of  y0  (y0  > 90 °)  the  radi- 
ation will  gradually  approach  that  of  the  nondirectional  radiator 
until  it  becomes  identical  with  that  of  the  spherical  radiator  of 
zero  order  for  y0  = 180°  . The  number  of  terras  which  must  be  cal- 

culated to  obtain  the  result  with  sufficient  accuracy  increases 
nearly  in  proportion  to  the  magnitude  of  kr0  . And,  indeed,  it 
turns  out  that  the  individual  terms  rapidly  decrease  as  soon  as  n 
has  become  greater  than  kr0  so  that  the  necessary  number  lies 
between  kr0  and  2 k r0  . This  is  connected  with  a general 
property  of  the  appearing  Bessel’s  functions,  which  exhibit  an 
essentially  different  behavior  if,  with  a fixed  argument,  the 
function  values  are  plotted  with  respect  to  the  index  n. 

As  long  as  n<x  f sn(x)  and  Cn(x)  vary  between  positive 
and  negative  values.  Ag  goon  as  n becomes  greater  than  x , 

SH(x)  decreases  sharply  and  monotonically  while  C„( x)  increases 
in  the  same  way.  The  3ame  thing  is  true  of  Vn(x)  and  VK(x)  . 

In  Figs,  % to  86 , UJx)  and  Vn(x)  are  represented  as  a function 
of  n for  x — 4, 10,  20,  and  40  . 

As  an  example  of  the  general  case  when  the  field  point  lies 
in  the  neighborhood  of  the  radiator  we  choose 

kr„  =•  5,  r « 2 r0  (i.e.i-r  = 10),  y0  = 30° 
and  calculate  the  table  found  on  page  110. 

For  the  pressure  amplitude  p at  the  field  point  f = 2r0,  y = 0 
wo  obtain  from  this 


e • a F • w i 


P - - ! 2,03  + 2,24  i \ » ■ 3,01 . 
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Fig.  70,  The  component®  of  the  Bound, 
pressure  for  ft-io0- 


Fig.  71,  The  components  of  the  sound 
pressure  for  c-so0. 
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Fig,  81.  Sound  pressure  amplitude  o»/p,>. for  a 
radiating  spherical  cap  (r.-ian. 


Fig.  82.  Sound  pressure  amplitude  (p-p.)  for  a 
radiating  spherical  cap  <r.- iwn. 


If  we  imagine  a small  spherical  radiator  with  the  same  F ■ w 
substituted  for  the  radiating  surface,  then,  in  the  absence  of  the 
rigid  sphere,  this  would  produce  at  the  same  field  point  (i.e., 

field  point  distance  = r/2) 
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the  pressure  amplitude 

c ■ a • F ■ w „ 
Po  ~ oTi  2 


(159©) 


The  complete  behavior  of  the 
pressure  amplitude  at  the  distance 
r — 2 r0  is  represented  as  a function 
of  y in  Fig.  88.  The  correspond- 
ing components  [R  4 ^J  = ^anPn(cosy)] 

are  given  in  Figure  87. 
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Fig.  87.  Sound  pressure  components 
of  the  radiating  spherical  cap 
()■.“ 3o°)  for  a small  field  point 
distance  6). 


Finally,  we  will  calculate  and 
illustrate  by  curves  of  constant 
pressure  amplitude  the  sound  field 
of  a point  radiator  located  on  a 
rigid  sphere  - the  field  being 
confined  in  the  neighborhood  of 
the  sphere.  A3  an  example  we 
choose  kr0  — 1 and  calculate  the 
pressure  for  the  field  points 
whose  distances  r from  the 
sphere  center  are  determined  by 

kr  - 1.5,  2,  3>  ky  and  5*  H«re,  we  first  calculate  the  expressions 

S,{kr)  + iC„(kr) 

(/,(i)  + »7H(i)-  (1591') 

(see  the  tables  on  page  112)  and  find  therefrom  ^ ~2Latn  and 

is  — + i . If  we  denote  the  undistorted  (i.e.,  existing  in 
the  absence  of  the  rigid  sphere)  pressure  amplitude  by  p0~-  . 

then  the  relative  pressure  amplitude  is  plp0  =>  iis  +i’»!  . In  the 

same  manner,  the  relative  amplitudes  for  y ~ 5°>  10°,  0tc.  would 
be  calculated  from  the  quantities  y.aHPK(w%y)  . 


an  = (2  n + 1) ' 


If  wo  collect  the  values  thus  found  and  plot  the  calculated 
field  points  together  with  the  corresponding  pressures,  we  can 
then  draw  the  curves  of  constant  pressure  by  interpolation.  One 
recognizes  from  Figure  89  that  the  sound  field  possesses  an 
essentially  different  character  because  of  the  presence  of  the 
rigid  sphez'e  although  the  wave  length  is  more  than  three  times 
the  diameter. 
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7.  The  disturbance  of  a sound  field  caused  by 
a rigid  sphere 

(a)  The  derivation  of  the  general  formula 


The  influence  which  a plane  rigid  infinitely  extended  wall  exerts 
on  the  sound  field  produced  by  a radiator  may,  as  is  well  known,  be 
easily  perceived  when  a second  radiator  is  introduced  which  is  located 
at  the  image  position  of  the  original  radiator  with  respect  to  the 
wall.  The  required  sound  field  (lying  on  the  side  of  the  original 
radiator)  then  results  simply  by  superposing  the  two  sound  fields 
produced  by  the  individual  radiators.  In  order  to  prove  this,  we 
have  to  show  that  the  boundary  condition  for  the  sound  field  on  the 
rigid  wall  (namely,  that  the  velocity  component  perpendicular  to 
the  wall  vanishes)  is  fulfilled.  Novr,  however,  it  immediately 
turns  out  from  the  symmetrical  positions  of  the  radiators  that  the 
velocity  vector  resulting  from  the  two  radiators  for  all  field 
points  on  the  wall  lies  in  the  symmetry  plane,  i.e.,  possesses 
no  component  perpendicular  to  the  wall.  Ljore  generally,  one  can 
state  that  if,  in  an  arbitrary  sound  field,  a surface  is  present 
such  that  the  velocity  component  perpendicular  to  this  surface 
vanishes  for  all  field  points  on  this  surface,  then  this  surface 
can  be  replaced  by  a rigid,  completely  reflecting,  surface  without 
any  change  in  the  sound  field.  Our  problem,  to  investigate  the 
influence  of  a rigi  -phcre  on  an  existing  sound  field,  is  then 
equivalent  to  the  problem  of  superposing  a second  field  on  the 
existing  sound  field  so  that  on  a prescribed  spherical  surface 
the  resulting  velocity  component  perpendicular  to  the  spherical 
surface  shall  /anish. 


..e  proceed  from  the  undistorted  sound  field  and  assume  that 
this  is  produced  by  a point  radiator  located  at  the  point  A with 
the  radiating  surface  F and  (constant)  velocity  amplitude  w . 
Then  the  sound  pressure  Vi  for  the  field  point  defined  by  the 
polar  coordinates  r , y is  given  by 


C 

Pi  - 


f T i>l  (<“  ( - i'  H i n i)  u*  k mi.  y 

UR  e 


(150) 


^ere  R is  the  distance  of  the  radiator  at  A from  the  coo:r- 
dinate  origin  0 . It  is  assumed  to  be  groat  compared  to  r . 

From  (160)  with  the  aid  of  the  relation 

c • ait  f r (l6l) 

there  results  the  radial  velocity  comi»nent 


(i6la) 


According  to  the  foregoing  considerations,  we  now  have  to  super- 
pose an  additional  sound  field  whose  radial  velocity  amplitude  t>* 
is  equal  in  magnitude  and  opposite  in  sign  to  q on  the  sphere 
r = r0.  The  problem  is  thus  reduced  in  principle  to  the  earlier 
calculation  where  the  velocity  distribution  on  a sphere  was  set 
forth.  % previous  considerations,  the  solution  can  be  stated 
immediately  if  the  velocity  distribution  on  the  sphere  is  repre- 
sented by  a series  of  spherical  harmonics.  In  order  to  bring  this 
about  we  proceed  from  the  relation 


oo 

tikrcm?  =2  (2n  + l)Pn(cosy) . 

n - 0 

By  (l6l)  it  then  follows  that 

o° 

or  - - ~ e““‘ (2ft  + L) Pn (cos y), 


n » 0 


(162) 

(162a) 


and  since 


|S.(*)V  _ _ U,(x) 


Clo2b) 


the  required  velocity  represented  in  the  desired  form  is: 


LiULfHut-  kR),  A... 

2 Air  A*r« 


oo 

y,(2»  + l)»"!71I(A;r0)i,B(cosy)  . 

n-=*0 


(163) 


By  (149)  we  can  immediately  specify  the  sound  field  produced  by 
. The  pressure  pa  at  the  field  point  r , y is 


A 

in  the 


_ eiU-’t  -kR  ! x!*) 

- UP 

00 

1 'srt(_2n+  l)t*f-/,(t»’(,)./».(c08y) 
' krjL  V.(&rt)  + iVnltrt)" 


[Sn(kr)  + iC\(kr)]. 


(164) 


n-0 

little  manipulation  thon  gives  the  total  sound  field  existing 
neigtiborhood  of  this  sphora  in  the  form 


A,  + « 


c • a • F • M> 
2 A A 


p 0>'t~kR  t rt.'S) 


[Sn(kr)  Vn(kra)-CAkr)Un(kr0) J. 


(165) 


(b)  The  sound  reflection  on  a rigid  sphere 

as  the  simplest  case,  we  first  investigate  the  pressure  ampli- 
tude produced  at  the  field  point  A (,i.e.,  at  ihe  position  of  the 
sound  radiator  itself)  by  reflection  on  the  sphere.  ,<e  will  compare 
tills  with  the  sound  pressure  amplitude  wtiich  is  reflected  perpendic- 
ularly by  a rigid  wall  at  the  same  distance  (as  from  the  sphere). 

If  wo'  introduco  the  reflection  factor  Z -■  1 3 1 


Attyleigh,  Theory  of  .Sound, 


ec. 


334. 
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by  the  relation 


8-&S™ 

n«0 


(2»  + l)Um(kr„) 
U.(kr0)  + iV,(lcr0)' 


(166) 


it  then  follows  from  (164)  that,  since  kr>l,  Sn{kr)  + iCn(kr)  becomes 

equal  to  and  cos  y becomes  unity,  one  has  for  the  reflec- 

tion amplitude  pr  due  to  the  sphere: 


coF  • w _ 
Vt  — ro  ^ > 


(167) 


while  the  reflection  amplitude  produced  by  the  rigid  wall  is  given  by 


Therefore 


c • a-  F • w 
~ iXR 


(168) 


Pr 

Pr 


(169) 


A very  simple  relation  for  the  reflection  factor  now  results, 
with  increasing  Jfcr0  , 2 approaches  the  value  £ so  that  for  the  larger 
values  of  kr0  one  can  with  good  approximation  set  the  reflection 
factor  equal  to  (169)  therefore  3tates: 


The  pressure  amplitude  reflected  by  a rigid  sphere  with  the 
radius  r0  at  a great  distance  R is  (for  the  larger  values  of  kr0 ) 
given  by  rJR  if  the  pressure  amplitude  reflected  by  an  in- 
finite wall  placed  at  the  same  distance  R is  set  equal  to  one, 

T0  represent  the  reflection  factor,  the  two  components  R and 
J are  computed  separately  and  plotted.  The  functions  Un(kr) 
and  necessary  for  this  are  easily  found  by  (147)  with  the 

aid  of  the  tables  for  S„(x)  and  C(x)  . The  result  of  the  cal- 
culation ia  represented  in  Fig.  90. 


If,  furthermore,  R + U is  drawn  as  a vector  in  the  complex 
plane,  one  then  has  the  phase  y(tgx  .HR)  a,a  & function  of  hQ. 

It  turns  out  that  the  phase  remains  practically  constant  while 
the  reflection  factor  increas-33  decidedly  (kr0-  1,2)  and  later 
increases  uniformly  when  the  reflection  factor  remains  constant 
(Fig.  91).  For  the  general  esse  ( y t 0 ),  we  define  the  dis- 

turbance factor  |3  (Fr0,  y | by 


x> 

>')  = L}ra  V (-ir 

« 0 


'(2>i  -1-  5 )V,{kr0) 
C,U'r0)  )•  1 


P.(cosy) 


(170) 


and  calculate  the  corresponding  components  for  kr0  ~ 1, 2, . . . 10  and 
can  then  represent  for  every  value  of  kru  the  complete  behavior  of 
the  disturbance  factor  as  a function  of  y (Figs.  92  to  99)*  It 
then  turns  out  that  for  y-  180°  (i.e.,  the  direction  opposite  to 

that  of  the  sound  source)  j^;  increases  more  and  more  with  kr0 

and  for  the  larger  values  of  kru  is  very  nearly  determined  by  4 kr0  . 
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Fig.  90.  Reflection  factor  Z of  a rigid  sphere 
and  its  components  (R,  J). 
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However,  even  for  a direction  opposite  to  that  of  the  radiator,  the 
disturbance  amplitude  is  small  compared  to  that  of  the  undistorted 
field. 


(c)  The  sound  field  in  the  neighborhood  of  the  sphere 


In  the  immedi&te  neighborhood  of  the  sphere  we  will  have  to  ex- 
pect a substantial  distortion  of  the  undistorted  field.  If  we  first 
inquire  as  to  the  behavior  of  the  pressure  amplitude  on  the  line 
connecting  the  radiator  with  the  center  of  the  sphere  on  the  side 
towards  the  radiator,  then  maxima  and  minima  will  appear  here  as 
they  did  in  front  of  a rigid  plane  wall.  As  an  example,  we  will 
investigate  the  case  kr0  = 2 when  we  calculate  the  pressure  for 
the  field  points  defined  by  kr  = 2, 3, . . . .10  . If  we  set  hr  = kr0 
in  (165),  then  due  to  the  relation  (from  129  and  147) : 


Sn{kr0)  Vn(kr„)  ~ Cn(kr0)Un(kr0)  = kr0, 


it  follows  that 


P = 


00 

|iaA  F ■ »«*<"»  - * ■ « + (2  n + 1) 

«-0 


i*  + lPH(co&  y) 

V*{kr<,)  + i ‘ 


(170a) 

(171) 


In  accordance  with  the  reciprocity  law,  this  is  in  agreement  with 
the  previously  discussed  case  where  the  radiator  was  assumed  to  be  on 
the  spherical  surface  and  the  field  point  was  assumed  to  be  at  a 
gi'eat  distance  (JSq.  156). 

For  convergence  reasons,  it  is  better  not  to  calculate  the  field 
composed  of  the  original  (undistorted)  and  the  distorted  field 
(p  = p. -t- p, ) by  (165)  but  to  first  calculate  the  distorted  field  P*  by 
(164)  and  to  add  the  undistorted  field  from  (160). 

Therefore,  we  first  calculate  the  disturbance  factor  by  (164) 
when  we  set  ^(cosy)  — 1 . Then 


3(kr0,  kr,  0)  =» 


1 (2»+ l)».*{/,(fcro) 

kr  U.(krB)'+iV,(krv) 

n—U 


\Sn(kr)  + iCK(kr)  J 


An(kr)  + *'4„(fcr)  ■ 


(171a) 


The  quantities  AH  + %Bn  for  kr0  = 2 and  kr  =>  2,  3,  ...  10  are  given 
in  the  following  table: 
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To  facilitate  later  calculations,  the  sums  _2i 
for  odd  n were  calculated  separately* 


for  even  n and 


The  values 


?>H2i+2W'*'| 


(171b) 


then  give  the  pressure  amplitudes  in  front  of  the  rigid  sphere  in  the 
direction  of  the  radiator.  were  the  undistorted  pressure  amplitude 
of  the  radiator  is  set  equal  to  unity.  We  see  that  the  pressure  ampli- 
tudes in  front  of  the  rigid  sphere  vary  in  a manner  similar  to  that  of 


the  reflection  on  a plane  wall  except  that  the  variations  with  in- 
creasing distance  from  the  sphere  decrease  very  rapidly.  (See  Fig, 
100.) 

The  complete  calculation  of  the  distorted  sound  field  is  then 
accomplished  with  the  aid  of  spherical  harmonic  tables  from 

Q{kr0,  kr,  y)  - &~<^y  ,=  ^(AH  + iBn)  Pn(oosy)  = R + %J.  (171c) 


The  components  R and  J are  first  calculated  as  a function  of  y 
for  kr  = 2,  3,  4, ...  10  and  are  plotted  in  Figures  101  to  104*  The 
quantities  fff+T* 

are  found  from  this  and  their  variations  about  the  undistorted  value 
one  are  represented  in  Figure  105*  Here  it  is  seen  that  the  magni- 
tude of  the  variations  decreases  with  increasing  kr  while  the 
number  of  the  variations  increases.  If  the  field  points  given  by 
kr  — 2,  3,  4, . . . io  are  numbered  corresponding  to  the  values  in  Figure 
105,  the  curves  of  constant  pressure  in  the  neighborhood  of  the 
sphere  can  be  drawn  by  interpolation.  The  result  is  given  in 
Figure  106.  The  greatest  value  (1.66)  lies  on  the  surface  of  the 
sphere  opposite  the  sound  source  ( ),  The  smallest  value 

(0.66)  likewise  lies  on  the  sphere  ( Y ~~  136°  ).  For  field  points 


at  a distance 


r>5  rn 


undistorted  amplitudes  remains  below  10$. 


the  difference  between  the  distorted  and 


r:  ! 
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Fig.  105.  The  deviation  of  the  distorted 
pressure  amplitude  in  the  neighborhood 
of  the  sphere  (radius  r ) from  the 
undistorted  amplitude  for  the  field 

point  (r.  r)  hr,  - 2;  hr  ~ 2,  3,  4 . . , 10, 


Fig.  106.  The  curves  of  constant  sound  pressure  in 
the  neighborhood  of  a rigid  sphere  with  the  radius  rc 
when  the  <>ound  source  is  at  a greet  distance  and 

irr0  . 2. 
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